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Resume 


Les nanofluides ont récemment été considérés comme l'un des fluides de travail les plus 
populaires en transfert de chaleur et en mécanique des fluides. Une estimation précise des 
propriétés thermophysiques des nanofluides est nécessaire pour l'étude de leurs performances 
en transfert de chaleur. Le coefficient de conductivité thermique et la viscosité sont quelques- 
unes des propriétés thermophysiques les plus importantes qui influencent directement 


l'application des nanofluides. 


Une telle augmentation de la conductivité thermique du nanofluide ne peut pas 
seulement être imputée à la plus grande conductivité thermique des nanoparticules en 
suspension. Mais, l'amélioration vient d'une combinaison de plusieurs facteurs physiques 
simultanés d'importance variable. À cette échelle nanométrique, le comportement thermique 
ne partage pas forcément les modèles des structures plus grandes, ce qui impose de faire plus 


de recherche pour en concevoir d'autres mieux adaptés afin d'expliquer ce que l'on observe. 


Dans ce contexte, l’objectif principal de ce travail de recherche est de développer un 
modèle numérique simple et précis capable de modéliser la conductivité thermique des 
nanofluides d’Oxydes métalliques/Eau, en utilisant un plan d'expérience (Design of 
Experiments DOE), ensuit ce modèle est testé et analysé par la méthode ANOVA pour 
déterminer les effets des différents facteurs et leurs interactions. Les résultats obtenus 
montrent que la fraction volumique des particules présente un effet plus significatif dans 
l’amélioration de la conductivité thermique des nanofluides par rapport aux autres facteurs. 
Contrairement à l'effet de la taille des particules, l'augmentation de la fraction volumique et de 
la température des particules possède un effet positif sur l'amélioration de la conductivité 
thermique des nanofluides. Les valeurs R° et RMSE de 0,99 et 0,0074 respectivement, 
indiquent ainsi que le modèle de TCR fournit une bonne corrélation et que les valeurs prédites 
et expérimentales s’adaptent parfaitement. Enfin, les résultats obtenus par DOE confirment la 
capacité du modèle proposé à prédire la conductivité thermique des nanofluides avec une 


grande précision. 


Le second objectif de ce travail est l’étude de l’application des réseaux de neurones 
artificiels (RNA) pour prédire les propriétés thermophysiques du nanofluide (Al:03/Eau). La 
conductivité thermique TCR (k»f/k) et la viscosité DVR (lf/u) du nanofluide sont 


respectivement modélisées. Les résultats de cette modélisation montrent que la meilleure 


architecture trouvée pour le taux de TCR et DVR est de 15 et 20 neurones dans la couche 
cachée respectivement. Le modèle (k,,/k;) fournit une bonne corrélation entre les valeurs 
prédites et expérimentales avec des valeurs de R° égales à 0,998 et 0,998 pour l'apprentissage 
et la validation. De plus, puisque les valeurs de MSE sont de 0,000018 pour l'apprentissage et 
de 0,0000162 pour la validation pour (lF/uF), donnant ainsi une excellente corrélation entre 


la réponse prédite et observée ; avec des valeurs plus élevées de R2 = 0,999 et de valeurs 
faibles de MSE = 0,000788 pour l'apprentissage et R? = 0,999 et SME = 0,000678 pour la 


validation. 


La comparaison avec les diverses corrélations empiriques prenant en compte les valeurs 
AARD de 0.0019 % pour (khf/k) et de 0.013 % pour (l/u) confirme que le modèle 
proposé à des performances supérieures dans la prédiction du TCR et du DVR des nanofluides 


Alumine/Eau et peut être considéré comme un outil pratique pour la tâche considérée. 


Mots clé: Nanofluides, Nanoparticule, Oxydes métalliques, Conductivité thermique, 


Viscosité dynamique, Plan d'expérience, RNA. 


Abstract 


Nanofluids have recently been considered one of the most popular working fluids in 
heat transfer and fluid mechanics. A precise estimation of the thermophysical properties of 
nanofluids 1s necessary for the study of their heat transfer performance. Thermal conductivity 
coefficient and viscosity are some of the most important thermophysical properties that 


directly influence the application of nanofluids. 


An increase in the thermal conductivity of nanofluids cannot only be attributed to the 
greater thermal conductivity of the suspended nanoparticles. But, the improvement comes 
from a combination of several simultaneous physical factors of varying importance. At this 
nanometric scale, thermal behavior does not necessarily share the models of larger structures, 


which requires doing more research to design appropriate models to explain what we observe. 


In this context, the main objective of this research work is to develop a simple and 
precise numerical model capable of modeling the thermal conductivity of metal oxides / water 
nanofluids, by using an experimental design (Design of Experiments DOE). This model is 
then tested and analyzed by the ANOVA method in order to determine different factors’ 
effects of and their interactions. The obtained results show that the volume fraction of the 
particles has most significant effect in improving nanofluids thermal conductivity compared 
to the other factors. Unlike the effect of particle size, increasing the volume fraction and 
temperature of the particles has a positive effect on improving nanofluids thermal 
conductivity. The R° and RMSE values are of 0.99 and 0.0074 respectively; indicating that 
the TCR model provides a good correlation hence fitting perfectly predicted and experimental 
values. Consequently, the results obtained confirm the capacity of the proposed model of our 


study to predict nanofluids thermal conductivity with great precision. 


The second objective of the present research is to study Artificial Neural Networks 
(ANNS) application to predict thermophysical properties of (ALO;/Water) nanofluid. The 
nanofluid thermal conductivity and its viscosity are thus respectively modeled. The results of 
this modeling show that the best architecture found for the rate of TCR (kr 4 kr) and DVR 
(Un . ur), is of 15 and 20 neurons in the hidden layer respectively. The model of (k» "4 ky) 
provides a good correlation between predicted and experimental values with values of R° 


equal to 0.998 and 0.998 for both learning and validation. In addition, the MSE values are of 


0.000018 for learning and 0.0000162 for validation. An excellent correlation 1s also noticed 
between the predicted and observed response for (u,F/uF), with thus a higher values of 
R2=0.999 and a lower values of MSE=0.000788 for training and R2=0.999 and 
SME=0.000678 for validation. 


Comparison with various empirical correlations, taking into account the AARD% 
values of 0.013 % for (uAF/u) and of 0.0019 % for (khf/k), confirms that the proposed 
model has greater performance in the prediction of TCR and DVR of Alumina/Water 


nanofluids and thus can be considered as a practical tool for the considered task. 


Keywords: Nanofluids, Nanoparticle, Metal Oxides, Thermal Conductivity, Dynamic 


Viscosity, Experimental Design, RNA. 
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Introduction générale 


Dans notre vie quotidienne, le transfert de chaleur et de l’énergie constitue la base de 
nombreux processus industriels. L’épuisement progressif des énergies fossiles conduit à 
améliorer et optimiser les rendements de ces échanges par de nouveaux procédés. Les 
« nanofluides » peuvent s’inscrire dans ce cadre pour être l’outil de transfert thermique et 
énergétique de nos jours et dans l’avenir. Ils apparaissent comme un sujet de recherche 
mondial en raison de leur potentiel, et 1ls peuvent être utilisés dans les échangeurs de chaleur, 
les systèmes énergétiques, les capteurs solaires, les appareils électroniques, en substitution à 


des liquides caloporteurs traditionnels. 


Une estimation précise des propriétés thermophysiques des nanofluides est nécessaire 
pour l'étude de leurs performances de transfert de chaleur. Le coefficient de conductivité 
thermique et la viscosité sont quelques-unes des propriétés thermophysiques les plus 


importantes qui influencent directement l'application des nanofluides. 


L’idée est d’améliorer les propriétés thermiques des fluides caloporteurs en y insérant 
une phase solide de conductivité thermique très élevée, ainsi les nanofluides présentent une 
conductivité thermique significativement plus élevée que les fluides de base qui augmente 
avec l’augmentation de la concentration de nanoparticules. Dans leur développement et leur 
exploitation, 1l convient de trouver un compromis entre de meilleures propriétés thermiques 
tout en maîtrisant leur viscosité (qui ne doit pas être trop importante) et leur stabilité qui doit 
être assurée, sans quoi les propriétés peuvent être modifiées. Différents modèles ont été 
élaborées en vue d’expliquer cette remarquable améliorations de transfert de chaleur, mais la 
confrontation entre ces modèles pour évaluera la conductivité thermique effective des 
nanofluides montre que ces derniers n’approximent pas de façon satisfaisante les expériences 
réalisées sur les nanofluides avec différentes variétés de nanoparticules et de fluides porteurs 


tels que : eau, éthylène glycol, … 


Plusieurs paramètres peuvent jouer un rôle dans l’efficacité des nanofluides comme la 
concentration et la nature des nanoparticules, leur densité, leur type, leur taille, leur aspect, la 


température, le fluide de base, la présence d’un surfactant, etc. 


Parmi les familles de nanoparticules, les nanomatériaux d’oxydes métalliques (Al203, 


CuO, TiO:,...) sont d’un intérêt majeur en raison de leurs excellentes propriétés thermiques 


intrinsèques et qui permettent d’obtenir des nanofluides plus performants que ceux préparés 
avec des nanoparticules métallique par exemple. C’est pourquoi notre travail de recherche se 


focalise sur les nanofluides d’oxydes métalliques. 


L'objectif de ce travail est de préciser les notions qui permettent d’évaluer et 
d’améliorer l’aptitude du modèle postulé à représenter les résultats expérimentaux. En un mot 
ces notions nous aideront à trouver le modèle qui représente le mieux possible les données 


expérimentales et qui devrait nous permettre de faire des prévisions de qualité. 


Il est nécessaire d’utiliser une approche méthodologique qui permet non seulement de 
réduire le coût de l’expérimentation, mais aussi d’établir une organisation optimale des 


expériences et séquences méthodologiques seront décrites dans ce mémoire. 


À cet égard, deux approches mathématique basées sur l’analyse des données 
expérimentales ont été développées afin de prédire les propriétés thermophysiques des 
nanofluides, elles peuvent être classées en deux grandes catégories: les méthodes des plan 
d’expériences, tels que le plan factoriel complet et les méthodes d’intelligences artificielles, 


telle que les Réseaux de Neuronaux Artificiels (RNA) 


Les méthodes des plan d’expériences sont très simples, permettent une bonne 
interprétation physique avec une minimisation du nombre nécessaire d'expériences qui permet 
un gain en temps et en coût financier, mais peu précises et modélisent difficilement le 
couplage. En revanche, la méthode des Réseaux de Neuronaux Artificiels est très précise, très 
souple, peut traiter un élément seul ou plusieurs éléments en réseau ou empilés, les formes 
arbitraires et le couplage. Cependant, elle est plus complexe, permet moins l'interprétation 
physique et nécessite un temps de calcul et un espace mémoire considérable. En outre, RNA a 
été largement utilisée pour prédire les systèmes non linéaires en raison des avantages notables 
de haute précision, de faible coût et de temps, par contre, la méthode de plan factoriel 


complète donne la meilleure estimation linéaire. 


Dans ce contexte, quatre chapitres ont été présentés dans cette thèse. Le premier 
chapitre représente les définitions et les connaissances fondamentales sur les nanofluides y 
compris les différents paramètres thermophysiques qui peuvent influer sur le comportement 


d'un nanofluide. 


Le deuxième chapitre est composé de deux parties. La première se focalise sur l’aspect 


conceptuel et terminologique pour introduire les notions de base sur les plans d’expériences, 
cela comprend les principales techniques mathématiques utilisées pour exprimer les variations 
des réponses en fonction des valeurs des facteurs. Les notions fondamentales des réseaux de 


neurones RNA et leurs propriétés font l'objet de seconde partie de ce chapitre. 


Le troisième chapitre présente la modélisation de la conductivité thermique effective 
des nanofluides d’oxydes métalliques. Ceci pour le développement d’une nouvelle corrélation 
et ceci incluant les paramètres les plus influents sur la conductivité thermique des nanofluides 
tel que : la concentration, la température et le diamètre de nanoparticules. Cette corrélation est 
générée par la méthode de plan factoriel complet. La discussion des résultats obtenue, 
l’évaluation de la corrélation générée qui est effectuée par la comparaison avec d’autres 
corrélations appliquées au même régime, ainsi que la détermination des erreurs absolues et 


réelles, ont été présentés dans ce chapitre. 


Le quatrième chapitre s’intéresse à l’application des réseaux de neurones artificiels 
(RNA) pour prédire les propriétés thermophysiques de nanofluide (A1:0; /Eau). Les résultats 


et les discussions de cette modélisation ont été exposés. 


Cette thèse est terminée par une conclusion générale qui résume le travail de recherche 
entamé ainsi que les principaux résultats obtenus et propose d’éventuelles perspectives pour 


des recherches futures. 


Chapitre I 
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L1 Introduction 

Durant ces dernières années, la consommation d'énergie a considérablement augmenté, ce 
qui a exhorté beaucoup de chercheurs à développer de nouvelles stratégies économiques afin de 
surmonter la menace d’une pénurie d'énergie. Ainsi, plusieurs systèmes énergétiques produisant un 
flux de chaleur très élevé ont été développés dans le domaine industriel, là où les fluides de 
refroidissement conventionnels ne sont pas efficaces (comme l’eau, l’éthylène glycol et l’huile). 
L’une des idées proposées et qui peut améliorer les propriétés thermo-physiques des fluides, est 
l'introduction des particules solides de très bonnes caractéristiques thermiques et de taille 
nanométrique au sein du fluide de base. Cette nouvelle génération de fluides porte le nom de 
Nanofluides (mélange colloïdal de nanoparticules de moins de 100 nm), ce terme a été proposé par 
Choi et Eastman [1] en 1995 et reste couramment utilisé pour caractériser ce type de suspension 
colloïdale. Les nanofluides ont des performances thermiques supérieures qui ne sont pas disponibles 
dans les fluides de transfert de chaleur traditionnels. Ainsi, l'objectif des nanofluides est d'obtenir 
des caractéristiques thermo-physiques les plus élevées aux plus basses concentrations possibles de 
nanoparticules (de préférence <1% en volume) et ceci par dispersion uniforme et une suspension 
stable de nanoparticules (de préférence de diamètre <10 nm) dans les fluides hôtes [2]. 

Ce chapitre se compose de deux parties principales. 

La première, présente une synthèse bibliographique permettant une meilleure appréhension 
des nanofluides, de leur production, leurs domaines d’applications et leurs propriétés physiques et 
thermiques. Par conséquent, elle est consacrée à la présentation d’une revue bibliographique des 
principaux travaux publiés dans le domaine de transfert thermique des nanofluides et de leurs 
propriétés thermo-physiques en général. Les types de nanoparticules présentés sont classés en trois 
groupes : nanoparticules métalliques, nanoparticules d’oxyde métalliques et nanoparticules non 
métalliques. 

La seconde partie est axée sur l’analyse de l’état des connaissances des propriétés thermo- 
physiques et du comportement thermique de ces nanofluides dans les systèmes thermiques. Une 
description de différents modèles utilisés pour l’estimation de ces propriétés thermo-physiques (tel 


que la conductivité thermique, la viscosité dynamique, la masse volumique, etc.) est aussi présentée. 
L2 Les nanofluides 


L.2.1 Définition et propriétés des nanoparticules 
Les nanofluides sont des dispersions de particules de taille nanométrique 


appelées nanoparticules, dans un fluide de base afin d’en améliorer certaines de ces propriétés. Ce 
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type de solution suscite un grand intérêt depuis la découverte de leurs propriétés thermiques 
particulières [3]. 

Le terme « nanoparticule » est utilisé lorsque la particule à au moins une de ses dimensions 
inférieures à 100 nm. Le rayon typique d’un atome étant de l’ordre de 0,1 nm, une nanoparticule 
peut être constituée de quelques dizaines à plusieurs centaines d’atomes. Par exemple, les 
nanoparticules les plus célèbres, les fullerènes, sont constituées de seulement 60 ou 70 atomes de 
carbone (Figure IL.1). Leur très petite taille donne aux nanoparticules des surfaces spécifiques très 


importantes [4]. 


Figure I.1: Illustration de la structure d’un fullerène nC60 [4]. 


La surface spécifique développée par les nanoparticules est 1000 fois plus importante que 
celle d’une microparticule (Figure 1.2). De même, le nombre d’atomes de surface augmente avec la 
diminution de la taille des nanoparticules et la réactivité, étant souvent liée à des propriétés de 


surface, cela explique l’intérêt développé pour les nanoparticules. 


Figure L.2: Illustration de la corrélation entre taille et surface spécifique [4]. 


La définition par l’échelle ne suffit cependant pas à décrire tous les aspects des 
nanoparticules. En effet, leur réactivité est également fortement liée à leur composition qui peut être 


très variée. On distingue par exemple les nanoparticules organiques (constituées majoritairement de 
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polymères), des métalliques (ex : nanoparticules d’or ou d’argent), d’oxydes métalliques (ex : 
Fe:03, Ti0;) ou encore des composites (ex : quantum dots CdSe-ZnS) [4]. 

Les nanoparticules, jonction entre l’échelle atomique et le matériau, ont des particularités 
structurales très spécifiques. En effet, contrairement aux matériaux classiques, les nanoparticules 
voient leurs propriétés changer en fonction de leur taille. Elles sont soumises à la fois aux 
phénomènes de la physique classique et à ceux de la physique quantique lorsqu'elles sont de très 
petites tailles. Les principales propriétés des nanoparticules sont : optiques, électriques, 


catalytiques, thermiques, magnétiques, mécaniques ou encore biologiques [4]. 


L.2.2 Les Nanoparticules et les fluides porteurs 

Une nanoparticule est un assemblage de quelques centaines à quelques milliers d'atomes, 
conduisant à un objet dont au moins l'une des dimensions est de taille nanométrique (1 à 100 nm). 
Les nanoparticules se situent à la frontière entre l'échelle microscopique et l'échelle atomique ou 
moléculaire. Elles ont un intérêt à la fois fondamental qu’appliqué et sont très étudiées de nos 
jours. 

D'une manière générale, les nanoparticules peuvent être classées selon leur forme en trois 
grandes catégories (Figure [.3) [5]. Elles ont généralement la forme sphérique et cylindrique, ou en 


forme de tube comme les nanotubes de carbone (NTC) [6]. 


Figure L.3: Images TEM de quelques exemples de nanoparticules : (a) Nanoparticules sphériques, 
(b) Nanotubes de carbone et Nanofeuilles de (c) graphène [5]. 
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Les nanoparticules les plus utilisées sont classées dans trois grandes familles, à savoir : 
o Les nanoparticules d’oxydes métalliques : 
1. L’oxyde d’aluminium (AlO3). 
2. L’oxyde de cuivre (CuO). 
3. L’oxyde de silicium (SiO). 
4. L’oxyde de Titanium (TiO). 
o Les nanoparticules métalliques : 
1. L’aluminium (Al). 
2. Le cuivre (Cu). 
3. L’or (Au). 
4. L’argent (Ag). 
o Les nanoparticules non métalliques : 
1. Les nanotubes de carbone (CNT). 
2. Le diamant (C). 
Les liquides de base les plus utilisés sont : 
1. L’eau. 
2. L’éthylène glycol (EG). 
3. Les huiles. 
4. Les fluides frigorigènes (R12, R22, etc.) [6,7]. 

Par ailleurs, l'agitation Brownienne, liée à la taille nanométrique des particules, minimise les 
problèmes de sédimentation rencontrés avec des particules de plus grande dimension. 

Les suspensions de nanoparticules dans un fluide (nanofluides), ont des caractéristiques 
thermiques intéressantes comparées aux fluides de base traditionnels [8], ainsi que des qualités 
indéniables d'amélioration du transfert thermique. La conductivité thermique des nanofluides peut 
être significativement plus élevée que celle des liquides purs (Figure 14). Toutefois, la 
quantification de la conductivité thermique des nanofluides et son évaluation nécessitent le recours 
à une instrumentation originale. De plus, Les mesures disponibles de la conductivité thermique de 


différents nanofluides ne sont pas concordantes [7]. 
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Figure I.4: Différence de conductivité thermique entre les fluides de base A (couramment utilisés) 


et nanoparticules B (couramment utilisées) [9]. 


1.2.3 Confection des nanoparticules et des nanofluides 


Les expériences effectuées sur les fluides caloporteurs avec des particules de taille 
micrométrique (10° m) ont montré que certains problèmes pouvaient être engendrés par l’utilisation 
de ces dernières ; tel que le problème de sédimentation des particules, l’augmentation de perte de 
charges etc. Tandis que pour des particules de l’ordre du nanomètre (10° m), ces problèmes sont 
pratiquement insignifiants [1]. 

Tout en supposant que la suspension est stable et durable, l’agglomération de particules est 
négligeable et aucun changement chimique du fluide ; la préparation des nanofluides peut se faire 
par l’une des méthodes suivantes : 

1. La méthode à une seule étape : basée sur la vaporisation d’un matériau solide sous vide 

(Figure LS), puis condensation directe de sa vapeur dans le liquide. 

2. La méthode à deux étapes : elle consiste d’abord à la fabrication d’une nanopoudre sèche, 


puis le mélange et la dispersion de cette dernière dans le liquide. 


Figure LS: Pyrolyse par LASER donnant de la nano poudre à mélanger au liquide ainsi que le 


dépôt direct dans le liquide par pulvérisation cathodique sous vide [6]. 


Les procédés de fabrication des nanoparticules sont de nature physique ou chimique [10, 11]. 
Ils font l’objet de beaucoup de recherches pour améliorer le coût de production qui reste toutefois 
élevé, vu les difficultés de mise en œuvre pour obtenir des particules de taille voulue. Chaque 
méthode a ses inconvénients et certaines sont inaptes à produire de grandes quantités de 
nanoparticules pour une production de masse vu les contraintes inhérentes au procédé. Les 
principaux moyens peuvent être : la mécano-synthèse (broyage à haute énergie), la précipitation 
chimique (réduction de sels métalliques par exemple), le CVD ou PVD, PECVD, la condensation 
des vapeurs par détente brutale dans des tuyères supersoniques, la pulvérisation cathodique à 
radiofréquence, la sublimation par ablation LASER, la pyrolyse LASER, la désintégration 


d'électrodes submergées par passage d'étincelles, etc. (Figure L.6) [3]. 
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Figure L.6: Procédés physique-chimique pour la préparation du nanofluide. 


1.2.4 État de l’art des nanofluides 
Compte tenu de l'importance accordée aux nanofluides, plusieurs études effectuées sur 


leurs propriétés physiques, tels que le coefficient de transfert thermique (la conductivité 
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thermique), la chaleur spécifique et la viscosité dynamique , ont prouvé que les nanofluides 
sont très différents par rapport aux fluides de base en présentant de meilleures caractéristiques 
thermo-physiques. Au fait, l'ajout de micros particules solides au sein des fluides de base 
augmente la conductivité thermique, puisque la conductivité thermique des solides est en 
général plus élevée que celle des fluides [6]. 

Plusieurs travaux ont révélé que la conductivité thermique des fluides contenant des 
nanoparticules peut augmenter de plus de 20% pour le mélange Eau-Al0O3 (1-5 % en volume) [1]. 
Même en très faible concentration (1% en volume), ils permettent d'augmenter la conductivité 
thermique de 150% dans le cas de nanotubes de carbone [12], de 40% pour des nanoparticules 
d'oxyde de cuivre [13] et plus de 20% pour des Oxydes d’Aluminium [12-17]. 

Nous présentons dans ce qui suit les études théoriques et expérimentales portant sur les 
nanofluides depuis leur découverte jusqu’à nos jours. 

De nombreuses recherches ont été menées depuis 2001 sur cette nouvelle classe de 
nanofluides afin de permettre une meilleure compréhension des mécanismes mis en jeu, et mettre 
ainsi au point des fluides caloporteurs plus performants. La forte conductivité thermique des 
nanofluides les désigne, en effet, comme des candidats potentiels pour le remplacement des fluides 
porteurs utilisés dans les échangeurs thermiques en vue d'améliorer leurs performances [14]. 

Certaines limitations susceptibles de réduire les performances des nanofluides utilisées à la 
place des fluides caloporteurs purs, ont été aussi étudiées. L'utilisation des nanofluides comme 
fluide caloporteur dans les échangeurs a montré que l’addition de nanoparticules dans un liquide 
augmente sa viscosité et de ce fait les pertes de charge [15]. De plus, le manque de stabilité dans le 
temps de certains nanofluides peut entraîner l’agglomération des nanoparticules et une modification 
de leur propriété de conduction thermique. Néanmoins, en l’état actuel des recherches, ces deux 
effets restent moins importants lors de l’utilisation de nanofluides que dans le cas de suspension de 
microparticules classiques [16]. 

En dehors des applications industrielles possibles, l’étude des nanofluides présente un intérêt 
scientifique pour la compréhension des phénomènes mis en jeu aux échelles nanométriques. 
Cependant, les propriétés thermiques des nanofluides ne concordent pas avec les prévisions fournies 
par les théories classiques décrivant les suspensions de particules solides dans un liquide. Malgré le 
nombre élevé d’études publiées sur le sujet ces dernières années, aucune théorie ne parvient à 
décrire convenablement les résultats expérimentaux obtenus sur les nanofluides [17]. Il est à noter 


aussi, que les résultats expérimentaux obtenus varient énormément d’une publication à une autre, ce 


qui ne facilite pas la comparaison entre ces derniers et les théories proposées. 


Généralité et revue bibliographique 


Donc, des recherches complémentaires des propriétés thermiques des nanofluides sont 
nécessaires pour valider ou infirmer les différentes théories proposées [7]. 

Maiga et al. [18] ont présenté une formulation numérique pour l’étude de transfert thermique 
par convection forcée des nanofluides (Eau/Al0O3) à l’intérieur d’un tube chauffé, et ce en 
proposant un modèle monophasique. Les résultats d’utilisation de ce modèle ont montré que le 
transfert de chaleur augmente avec l’augmentation des fractions volumiques. 

Namburu et al. [19] ont aussi présenté un modèle monophasique dans le but d’étudier les 
caractéristiques de transfert de chaleur d’un nanofluide dans un tube circulaire où l’écoulement est 
en régime turbulent. Leurs résultats ont montré qu’il ya une augmentation de viscosité et du 
nombre de Nusselt pour de petits diamètres. 

D’autres chercheurs ont proposé des modèles différents et c’est à cet effet que Behzadmehr 
et al. [20] et Mirmasoumi et al. [21] ont utilisé le modèle diphasique d’un nanofluide dans un tube 
en régime de convection forcée. La comparaison des résultats calculés avec les valeurs 
expérimentales montre que le modèle diphasique est plus précis que le monophasé. 

Au cours de la dernière décennie, l’accent des recherches sur les nanofluides est plus évident 
comme le montre la figure I[.7, qui illustre le nombre de publications portant sur les nanofluides 
depuis 2010. Ces études portent notamment sur leur préparation, leur caractérisation, la mesure de 
leurs propriétés physiques et leur utilisation dans diverses applications. Ces données présentées en 
histogrammes ont été obtenues en faisant une recherche du mot « nanofluides » dans la base de 
données Scopus par rapport aux titres, aux résumés et aux mots clés pour la période présentée. Les 
résultats de cette recherche montrent qu’environ 3165 articles ont été publiés en 2019 seulement, et 
cette tendance devrait surement augmenter au cours des prochaines années. 

Ces articles comprennent des applications sur les capteurs solaires [22, 23], les échangeurs de 
chaleur [24, 25], les caloducs [26], le refroidissement des radiateurs [27] et le refroidissement 
électronique [28] ; en plus d’une parade des travaux concernant divers propriétés thermophysiques 


des nanofluides [29, 30]. 
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Figure L.7: Nombre de publications pour les dix dernières années [9]. 


La figure L8 ci-dessous présente le nombre de publications dans le domaine des 


nanofluides des dernières années ainsi que les différents domaines d’utilisation. 
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Figure L.8: Nombre de publications des nanofluides NF, et de leurs conductivités thermiques TC et 


le transfert de chaleur par convection CHT au cours des dernières années (Web of Science) [31]. 


Cette figure illustre le bilan mondial des publications concernant les nanofluides NF, y 
compris leur conductivité thermique (TC), la convection naturelle et le transfert de chaleur par 
convection (CHT), ainsi que les domaines d’application des NF (obtenues à partir de la plateforme 
Web of Science, Clarivate, de 2012 à 2019). Ces publications comprennent tous les types d'articles 


(études expérimentales, théoriques et numériques). Comme le montre cette figure, la recherche sur 
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la NF a continué de croître, à l'exception de la légère déclinaison de l'année 2018 (pour la première 
fois). On constate que les travaux sur la CT (propriété thermo-physique principale de la NE), sont 
les plus fréquemment publiés. En revanche, les publications sur le CHT sont restées très en retard 
par rapport à celles des TC et NF, avec des hauts et des bas pour la période considérée. On constate 
aussi que la convection naturelle, (qui ne nécessite pas de coût supplémentaire d'énergie de 
pompage et qui est également utilisée comme mode de refroidissement clé), n'a pas fait l'objet d'une 
attention considérable de la part des chercheurs. En outre, le nombre des études de convection 
naturelle de NF sont significativement moins que celles des TC et NF, en plus elles sont légèrement 
inférieures à celles de CHT (convection forcée). Ainsi, il existe un besoin évident de travaux plus 
accomplis sur la thermoconvection de NE, et ceci dans diverses conditions thermiques, de 


concentration et dans différentes géométries [31]. 


1.2.5 Revue bibliographique 

Nous présentons dans cette partie des études expérimentales et numériques disponibles dans 
la littérature, sur les propriétés principales des nanofluides et leur impact sur le transfert de chaleur 
tel que la conductivité thermique, la chaleur spécifique et la viscosité dynamique. Diverses 
corrélations ont été proposées et comparées avec les expressions classiques. 

Les premiers travaux traitant le transfert de chaleur en présence des nanoparticules ont 
pratiquement commencé en 1995 avec l'étude de Choi [1], qui ont permis plus tard de déterminer 
les propriétés thermophysiques des nanofluides. Par la suite, plusieurs études expérimentales et 
numériques ont été menées par différents chercheurs pour comprendre et analyser le comportement 
thermo-hydrauliques des nanofluides. Ces études peuvent être classées en trois catégories selon le 
type de matière des nanoparticules : 

" Les nanoparticules d’oxyde métallique. 
" Les nanoparticules métalliques. 
" Les nanoparticules non métalliques. 
1.2.5.1 Les nanoparticules d’oxyde métallique 
L'idée d'insérer, au sein du fluide, des particules de taille nanométriques pour espérer 
augmenter la conductivité thermique effective du mélange est apparue depuis Maxwell en 1881 
[32]. Plusieurs types de matériaux peuvent être utilisés pour l'obtention de nanoparticules 
sphériques. Ces nanoparticules peuvent être fabriquées à l’aide de métaux solides ou d’oxyde 
métallique. Les oxydes métalliques sont faciles à produire, ils sont chimiquement stables et 


faciles à disperser dans l’eau. 
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Ÿ Les nanoparticules d’oxyde d’Aluminium 
L'essentiel des travaux sur ce type de nanofluides a porté sur la mesure de la conductivité 
thermique en fonction des différents paramètres caractérisant les nanoparticules (taille, fraction 


massique et température de l'ensemble). 


Le premier travail publié sur l’amélioration de la conductivité thermique utilisant l’ox yde 
d’aluminium a été mené par Masuda et al. en 1993 [33]. En dispersant des nanoparticules (A1,0;) 
dans l’eau dont la taille est de 13nm de diamètre. Ils ont obtenu une amélioration de la 
conductivité thermique de 30% pour une fraction volumique @ = 4.3%. Ce résultat a été suivi 
par une autre étude publiée par Eastman et al. [8] qui ont trouvé une amélioration de 30% de la 
conductivité thermique de nanoparticules (AI,0;) de taille 33nm en suspension dans l’eau pour une 
fraction volumique de 5%. 

Les travaux de Wang et al. [34] ont porté sur l’étude numérique de la convection 
naturelle dans une cavité bidimensionnelle pour différents nombres de Grashof et pour 
différentes concentrations en volume. Ils ont obtenu une amélioration de la conductivité thermique 
de 40% pour (Ethylène Glycol /AI,0;) de fraction volumique 8% et une augmentation de 14% 
pour (Eau/AI,0;) avec une fraction volumique de 4.5%. 

Lee et al. [35] ont mesuré la conductivité thermique des nanofluides d’oxyde de cuivre 
(CuO) et d’oxyde d’aluminium (Al,0;) dans deux liquides de base : l’eau et l’éthylène glycol 
obtenant ainsi une amélioration de 10% de la conductivité thermique pour le nanofluide 
(Œau/AL0,) de fraction volumique 4.3% et une amélioration de 20% pour (Ethylène Glycol/AL,0O;) 
d’une fraction 6 = 5%. Ces derniers résultats diffèrent considérablement de ceux obtenus par les 
précédents auteurs malgré leur utilisation de nanoparticules de même taille. 

En 2003 Das et al. [36] furent les premiers à investiguer l'influence de la température sur 
l'accroissement de la conductivité thermique des nanofluides. Ils ont relevé une augmentation de 
la conductivité thermique lorsque la température varie de 21°C à 51°C pour les deux nanofluides 
(Œau/AI,0;) et (Eau/CuO) avec des diamètres de particules de 38,4nm pour le premier et de 
28.6nm pour le second. Les résultats expérimentaux dévoilent que la conductivité thermique 
augmente linéairement avec l’augmentation de la température et dépend de la fraction volumique 
des nanoparticules dans l'eau. 

Xie et al. [37] ont examiné l’influence du pH de la solution fluide de base pour 
différentes tailles de nanoparticules sur l’amélioration de la conductivité thermique du nanofluide 


d’(ALO;) dispersées dans l’eau, l’éthylène glycol et l’huile. Ils ont découvert que 


l’amélioration de la conductivité thermique diminue avec l’augmentation de pH. 
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En 2006, Pracher et al. [38] ont mesuré la conductivité thermique à 85°C d'un 
nanofluide constitué d'eau et de particules d'AI,O, de diamètre 10nm avec une fraction volumique 
de 0.5%. Ils ont noté qu’à cette température la conductivité du nanofluide était égale à deux fois 
celle de l'eau. 

Li et Peterson [39] ont utilisé des nanoparticules ALLO; de diamètre 36nm avec une 
fraction volumique de 10% dans l’eau pure à une température de 34.7°C. Ils ont relevé une 
amélioration de 30% de la conductivité thermique du nanofluide. Sur la base de ces résultats 
expérimentaux, Li et Peterson ont proposé une relation permettant le calcul de la conductivité 


thermique du nanofluide (Eau/A1,0O;), à savoir : 


IT = 0,7644819 + 0,018688867T — 0.4621 (Li) 
f 


Où : 
kar, k, @ et T désignent respectivement les conductivités thermiques du nanofluide et du fluide 

de base, la fraction volumique des nanoparticules et la température. 

Ÿ” Les nanoparticules d’oxyde de cuivre 

En utilisant des nanoparticules de CuO, de diamètre 36nm avec une fraction volumique de 
5% dans de l’eau, Eastman et al. [8] ont obtenu une amélioration de la conductivité thermique 
de 60%, deux fois supérieure à la conductivité thermique obtenue dans les mêmes conditions 
avec des nanoparticules d'AI,O;,. 
Lee et al. [35] ont trouvé, en mesurant la conductivité thermique des nanofluides avec des 
nanoparticules de CuO de diamètre 18,6 nm dans de l'eau et dans l’éthylène glycol, que la 
conductivité thermique du nanofluide utilisé augmente linéairement en fonction de la fraction 
volumique des nanoparticules. 
Das et al. [36] ont refait avec des nanoparticules de CuO les mêmes expériences développées 
pour des nanoparticules d'ALO, et ils sont parvenus à des résultats analogues, à savoir que la 
conductivité thermique du nanofluide étudié varie toujours linéairement avec la température. 

Li et Peterson [39] ont mené des expérimentations avec des nanoparticules CuO de 
diamètre 29nm dans l’eau pure avec une fraction volumique 10% et une température de 34.7 €, 
et ils ont relevé une amélioration de 52% de la conductivité thermique du nanofluide. De plus, 
ils ont proposé comme pour le cas des nanoparticules d'ALO, la relation suivante pour le 


nanofluide (Eau/CuO) : 


BIS _ 376108 @ + 0.017924T — 0.30 (2) 
k (0) 


f 


Généralité et revue bibliographique 


Ÿ” Nanoparticules de dioxyde de Titane 

Les nanoparticules d’AI,O; et de CuO, plus faciles à trouver dans le commerce, ont été 
les nanoparticules les plus étudiées. Par contre, d’autres nanoparticules ont été plus tardivement 
étudiée, telle que l’oxyde de Titane (TiO,), qui ont fait l’objet d’investigation par Pack et Choïi 
[40]. Les résultats obtenus ont montré une amélioration de la conductivité thermique de 
10.7 % concernant les nanoparticules de 27 nm dispersées dans l’eau à 4,35 % de fraction 
volumique. Cette valeur est beaucoup plus faible que les 32% obtenues pour le nanofluide 
(Œau/ALO;) avec la même concentration des nanoparticules. Murshed et al. [41] qui ont mesuré la 
conductivité thermique des nanoparticules (Ti0,) sphérique et cylindrique, ont trouvé que le 
nanofluide contenant des nanoparticules sphériques de 15 nm conduit à une faible amélioration de 
la conductivité thermique. Tandis que, pour des nanoparticules cylindriques de 40 nm, ils 
obtiennent une amélioration allant jusqu’à 33% pour une fraction volumique de 5%. 

Le tableau I.1 regroupe des résultats expérimentaux, donnant la conductivité thermique des 
nanofluides obtenus par plusieurs auteurs, et ceci pour des nanoparticules d’oxydes métalliques de 
différentes tailles dispersées dans différents liquides porteurs (l’eau, l’éthylène glycol et l’huile) et 
pour différentes valeurs de fractions volumiques. Ce tableau révèle des discordances entre les 


résultats obtenus par les différents auteurs et avec le même type de nanoparticules et fluide de base. 


1.2.5.2 Les nanoparticules métalliques 
Ÿ”_ Les nanoparticules de cuivre 

Les nanoparticules métalliques ont des conductivités thermiques plus grandes que celles des 
oxydes métalliques correspondantes. La conductivité thermique du cuivre à la température ambiante 
est environ 700 fois plus grande que celle de l’eau et environ 3000 fois plus grande que celle de 
l’huile du moteur. Donc, s’attendrait normalement à trouver des valeurs de conductivité thermique 
de même ordre que celles obtenues avec des oxydes métalliques correspondants et ceci pour des 
fractions volumiques beaucoup plus faibles. Xuan et Li [42] ont présenté une étude démontrant que 
la conductivité thermique augmente avec l’augmentation de la fraction volumique. Le rapport de la 
conductivité thermique du nanofluide (Eau/Cu) et fluide de base augmente de 75% et celui de 
(Huile/Cu) de 45% quand la fraction volumique des nanoparticules varie entre 2.5% à 7.5%. 
Théoriquement, Xuan et al. [42] ont utilisé le modèle de Hamilton-Crosser pour étudier la 
conductivité thermique des mêmes nanofluides et ils ont obtenu des résultats très proches des 


résultats expérimentaux. 
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Tableau I.1: Résultats donnant les valeurs maximales des conductivités thermiques obtenues pour 


différents nanofluides avec des nanoparticules sphériques d’oxydes métalliques et différents fluides 


de base. 
Taille des Fraction L'amélioration 
Nanofluides Nanoparticules (nm) | Volumique ( %) | Maximum de k (%) | Référence 


ALO, + eau 


Masuda et al. 1993 


ALO, + eau 


Eastman et al. 1997 


ALO, + eau 


ALO, + eau 


Pak et cho. 1998 


Wang et al. 1999 


ALO, + eau 


ALO,+ eau 


Lee et al. 1999 


Das et al. 2003 


Xie et al. 2002 


ALO, + eau 


Prasher et al. 2006 


ALO, + eau 


Krishnamurthy et 
al.2006 


CuO +eau 


CuO + eau 


Eastman et al. 1997 


Das et al. 2003 


CuO + eau 


Lee et al. 1999 


CuO + eau 


ALO, +EG 


Wang et al. 1999 


Wang et al. 1999 


ALO, +EG 


Lee et al. 1999 


ALO, +EG 


CuO + EG 


Xie et al. 2002 


Liu et al. 2006 


CuO + EG 


Wang et al. 1999 


ALO, +huile 


Wang et al. 1999 


ALO, + huile 


ALO, + huile 


TiO, + eau 


Xie et al. 2002 


Wang et al. 1999 


10.7 Pak et Choi. 1998 


TiO2+ eau 


33 Murshed et al. 2005 


En 1997 Eastman et al. [8] ont observé une amélioration de 45% de la conductivité 
thermique par l’ajout de moins de 0.055% en volume des nanoparticules de cuivre de 35 nm dans de 
l'huile. Également en 2001, ils ont trouvé une amélioration de 40% de la conductivité thermique 
pour le nanofluide (Cu/Ethylène Glycol) contenant moins de 0.3% en volume de nanoparticules de 


diamètre inférieur à 10 nm. 
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Ÿ _ Les nanoparticules d’Or et d’Argent 

La taille moyenne des nanoparticules d’or (Au) est de 10-20 nm, tandis que, la taille 
moyenne des nanoparticules d’argent (Ag) est de 60-80 nm. 

Patel et al. [43] ont étudié les nanofluides (Au/Thiolate/Toluène) et (Au/Citrate/Eau) à une 
température de 60 °C ; les résultats obtenus montrent qu’il y’a une amélioration de 8.8% de la 
conductivité thermique pour le premier (avec une fraction volumique de nanoparticules ® —0.011%) 
et une amélioration de 8,3% de la conductivité thermique pour le second (avec un 6 = 0.00026%). 
Par contre, avec une solution de (Ag- Citrate/Eau), ils sont parvenus à une amélioration de 4.5% 
pour @=0.001%. Pour le nanofluide (Ag/Eau) la conductivité thermique obtenue est plus faible que 
celle du nanofluide (Au/Eau), bien que l'argent soit meilleur conducteur de la chaleur que l'or. Cette 
nette amélioration de la conductivité thermique lors de l'utilisation des nanoparticules d'or par 
rapport aux particules d'argent pourrait être attribuée à la taille des nanoparticules d'or qui sont de 
l'ordre de 3 à 8 fois plus petites que celles d'argent. Kumar et al. [44] ont rapporté, quant à eux qu’il 
y à une amélioration de la conductivité thermique de 20% pour @ = 0.00013% d’or dans l’eau. 

Le tableau [2 regroupe les résultats expérimentaux, donnant la conductivité thermique des 
nanofluides obtenus par plusieurs auteurs avec des nanoparticules métalliques de différentes tailles 
dispersées dans différents liquides porteurs (l’eau, l’éthylène glycol et l’huile) et ceci pour 


différentes valeurs des fractions volumiques de nanoparticules. 


Tableau IL.2: Valeurs maximales de la conductivité thermique des nanofluides utilisant des 


nanoparticules métalliques obtenues par différents auteurs. 


Taille des Fraction L'amélioration 
Les Nanofluides nanoparticules | volumique maximum de k% Référence 
Cu + huile 35 0.055 45 Eastman et al. 1996 
Cu + eau 100 TS 75 Xuan et Li 2000 
Cu + huile 100 7.5 45 Xuan et Li 2000 


Cu + EG 10 0.3 40 Eastman et al. 2001 
Au + toluène 15 0.011 8.8 Patel et al. 2003 
Au + eau 15 0.00026 8.3 Patel et al. 2003 


Ag +eau 70 0.001 45 Patel et al. 2003 
Fe + EG 10 0.55 18 Hong et al. 2006 


1.2.5.3 Les nanoparticules non métallique : Les nanotubes de carbone 


Ÿ” Structure 


Les travaux d’Harold Kroto, Richard Smaley et Robert Curl (Prix Nobel de Chimie en 1996) 
ont abouti à la découverte d’une nouvelle forme allotropique du carbone, les fullerènes. L’étude 


de ces nouvelles structures formées de 60 atomes de carbone arrangés sous formes de sphères avait 
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très vite suscité un vif intérêt dans le milieu des matériaux nanostructurés. C’est en étudiant ces 


nouvelles structures, synthétisées alors par arc électrique, qu’en 1991 Sumio lijima observa pour la 
première fois les nanotubes de carbone (NTC) multi-feuillets noté MWNT (Multi-Walled 
NanoTube). Ils sont constitués de deux ou plusieurs couches de Graphene. En 1993, il découvrit 
également des nanotubes mono-feuillet. Le NTC mono-feuillet noté SWNT (Single Walled 
NanoTube), se présente sous la forme d’un Graphene, enroulé sur lui-même, décrivant une 
géométrie tubulaire. Si le NTC n’a pas subi de détérioration, une demi-sphère de fullerène forme 
généralement l’une de ses deux extrémités, tandis que la seconde peut être occupée par une 
nanoparticule métallique utilisée comme catalyseur pour la synthèse des NTC. Suivant l’orientation 
des motifs hexagonaux de carbone par rapport à l’axe du tube, il est possible de différencier les 


divers types de SWNT (Figure L.9). 


Figure L.9: Différentes structures des nanotubes identifiées [7]. 


Ÿ_ Propriété thermique des nanotubes de carbone 

Les travaux théoriques et expérimentaux ont permis la détermination de la conductivité 
thermique de ces nanotubes de carbone (3000 W/m°K pour les nanotubes de carbone multi- 
feuillets ou multiwall (MWNTCSs) et 6000 W/m°K pour les nanotubes de carbone single-wall 
(SWNTCSs) [45]. Ce sont les études portant sur les nanofluides à base de nanotubes de carbone 
qui ont montré les augmentations de conductivité thermique les plus élevées. Ces dernières sont 
de l’ordre de plus de 160% pour des concentrations en volume de 1% (Choï en 2001). Une manière 
simple de représenter l’effet de nanotubes de carbone en solution est de considérer les nanotubes en 
suspension comme étant de nano-ponts thermiques susceptibles de conduire rapidement l’énergie 


d’un point à l’autre de l’échantillon comme le montre la figure [.10. 
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Figure L.10: Explication schématisée de la création d’un réseau de ponts thermiques (trait rouge) 
créé par une dispersion de nano-fibres (trait noir) dans un liquide (Beige) qui permet la propagation 


de la chaleur d'une zone chaude (ovale rouge) vers une zone froid [7]. 


Cette augmentation de conductivité thermique peut être prédite par la théorie du milieu 
effectif qui prévoit de fortes augmentations de conductivité thermique pour des fluides contenant 
des nanoparticules ayant un grand rapport d’aspect (Keblinski 2005). Cette théorie permet de 
calculer l’augmentation de la conductivité thermique du liquide en fonction du rapport d’aspect 
des nanotubes utilisés. Elle prédit que la conductivité thermique du nanofluide devrait augmenter 
fortement avec le rapport d’aspect des nanotubes. 

L'augmentation de conductivité thermique prévue par cette théorie (jusqu'à 6000%) est 
largement supérieure à celle obtenue expérimentalement (150%). Cela pourrait-être expliqué par 
la forte résistance inter-faciale entre le nanotube et le liquide, résultant de la faible affinité entre 
les liquides et les nanotubes de carbone qui provoquent une déplétion des molécules de liquide 
environnant autour du nanotube. Ceci a pour effet de réduire le couplage thermique entre les 
nanotubes et le liquide environnant, ce qui réduit la conductivité thermique du mélange. Afin 
de limiter cet effet, des travaux sont actuellement en cours pour obtenir une meilleure adéquation 
entre le nanotube et le fluide porteur. Ces travaux de recherche permettent d'une part d’améliorer 
la stabilité de la suspension et d'autre part l'obtention de nanofluides meilleurs conducteurs de la 
chaleur. En revanche, les molécules greffées à la surface du nanotube peuvent également induire 
une résistance inter-faciale supplémentaire [46]. 

La principale difficulté rencontrée avec les nanofluides à base de nanotubes de carbone la 
dispersion de manière stable dans le temps les nanotubes de carbone dans le liquide de base 
Wang et al. [47]. Le rapport d’aspect des nanotubes, les méthodes de préparation et 
d'ensemencement sont différents d’une publication à une autre, entraînant une grande différence 
entre les valeurs des conductivités thermiques mesurées pour les nanofluides avec des nanotubes 
de carbone. Ainsi la conductivité thermique augmente de 10% pour une concentration de 1% de 
nanotube dans l’eau (Xie 2003), de 38% pour une concentration de 0.6% de nanotubes dans 


l’eau (Assael 2005) ou encore de 150% pour une concentration de 1% de nanotubes dans 


l’Ethylène Glycol (Choi 2001). 
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L'utilisation de nanotube de carbone en suspension Choi et al. [2], où de multi couches de 
nanotubes de carbone de taille de 25 nm et une longueur de 50um sont dispersées dans de l’huile 
Polyoléfine (PAO) à la température ambiante. Une amélioration de rapport de conductivité 
thermique de plus de 160% a été observée avec un ® =1%, cela indique que l’interaction entre 
les fibres conduit à une plus grande valeur de la conductivité thermique. 

D'autres études ont conduit à des résultats différents pour des nanotubes de carbone 
dispersés dans de l’eau. Xie et al. [48] y ont dispersé des (MNTC). L'amélioration de la 
conductivité thermique du nanofluide obtenue est de 20% pour une concentration de 1% de 
volume. 

Wen et Ding [49] ont mesuré l’amélioration de (MNTCSs) en suspension dans l’eau en 
fonction de la température et de la concentration. Ils ont trouvé, qu’en dessous de 30°C, la 
conductivité thermique du nanofluide augmente linéairement avec l’augmentation de la 
température. L'amélioration maximum obtenue est de 31% pour ® = 0.84% et à la température de 
60°C. La même expérience a été menée par Ding et al. [50] qui ont trouvé une amélioration de 
80% pour 6 = 1% et à une température de 30°C. Assael et al. [51] ont obtenu 38% pour une fraction 
volumique de 0.6% de nanotubes dans l’eau ou encore de 150% pour une concentration de 1% de 


nanotubes dans l’éthylène glycol. 


Tableau [.3: Valeurs mesurées de la conductivité thermique des nanofluides avec des nanotubes de 


carbone pour différents fluides de base et pour différentes valeurs de la fraction volumique. 
Les Nanofluides 


Lu 
MWCNTSs/PAO 160 
MWCNTSs/eau 6 Xie et al. 2003 


Référence 
Choi et al. 2001 


MWCNTSs/Decène 


Xie et al. 2003 


MWCNTS$s/eau 


Wen et Ding 2004b 


MWCNTs/PAO 


1 
0.84 
0.35 


Yang et al. 2006 


DWCNTSs/Eau 
MWCNTSs/Eau 
MWCNTSs/Huile synthétique 
MWCNTSs/EG 


1 

1 

1 

1 8 Assael et al. 2004 et 2005 
2 30 Liu et al. 2005 

12.4 Liu et al. 2005 


L.2.6 Domaines d’applications des nanofluides 


Les nanofluides peuvent être employés pour améliorer le transfert de chaleur et delà 
l'efficacité énergétique dans plusieurs systèmes thermiques. Dotés de caractéristiques thermiques 
exceptionnelles, les nanofluides constituent une solution attractive dans de nombreuses applications. 
Selon leur utilisation, les nanofluides sont classés en nanofluides chimiques, nanofluides 


tribologiques, nanofluides surfactants et de revêtement, nanofluides caloporteurs, nanofluides 
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environnementaux (nettoyage de la pollution), nanofluides de procédé d’extraction, nanofluides 
pharmaceutiques et biologiques et nanofluides médicaux (interaction tissu-cellule). 
Dans cette section, quelques exemples d’applications seront présentés tel que le transport, la 


micro-électronique, la défense, le nucléaire, l'espace et le biomédicale [52]. 


1.2.6.1 Refroidissement des systèmes thermiques 

Le refroidissement des systèmes thermiques a été étudié par plusieurs chercheurs, comme 
Tzeng et al. [53] qui ont dispersé des nanoparticules de CuO et Al:03 dans l’huile de 
refroidissement du moteur de transmission, leurs résultats expérimentaux montrent que CuO a le 
meilleur effet de transfert de chaleur. Également, Akhavan-Behabadi et al. [54] ont révélé une 
amélioration significative du transfert de chaleur par condensation de l’écoulement du nano- 
réfrigérant R60O0/POE/CuO dans un tube horizontal lisse. 

Les nanofluides ont un point d’ébullition élevé et par conséquent peuvent être utilisés pour 
augmenter la température de fonctionnement du liquide de refroidissement, puis rejeter davantage 
de chaleur par le biais du système de refroidissement existant ,comme le montre les études 
effectuées par Yang et al. [55]. Ces derniers ont mené des expériences sur un tube ondulé afin de 
décrire les performances en ébullition du nanoréfrigérant MWCNT/R141b. De même des 
expériences sur les nanoréfrigérant R600a-CuO-POE dans un tube horizontal ont été menées par 


Darzi et al. [56]. 


1.2.6.2 Refroidissement des systèmes électroniques 

Un refroidissement efficace est l’un des principaux défis dans le domaine de l’électronique. 
Il à été démontré que les microprocesseurs atteignaient le double de leur capacité de calcul 
lorsqu'on utilise un refroidissement à base de nanofluides. Les dispositifs électroniques avancés 
font face à des défis de gestion thermique dus au niveau élevé de génération de chaleur et à la 
réduction de la surface disponible pour l’évacuation de la chaleur. Parmi les études effectuées, on a 
Tsai et al. [57] qui ont employé un nanofluide à base d’eau pour refroidir une unité centrale de 
traitement dans un microordinateur, ils ont remarqué qu’au même volume de charge, il y a une 
réduction significative de la résistance thermique du caloduc avec nanofluide par rapport à l'eau. De 
leur côté, Ma et al. [58] ont étudié l’effet du nanofluide sur les possibilités de transport de la chaleur 
d'un caloduc d’oscillation, et ils ont démontré qu’à la puissance d’entrée de 80 W et seulement pour 
une concentration volumique de w=1%, la différence de température entre le vaporisateur et le 
condensateur est de 16.6°C. 

1.2.6.3 Refroidissement des systèmes d'applications militaires 

Les exemples d'applications militaires incluent le refroidissement de l'électronique de 
puissance et les armes d’énergie dirigée. Ces derniers impliquent des flux thermiques très élevés 
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(Q>500 à 1000 W/cm?). Les nanofluides sont susceptibles de fournir le refroidissement requis dans 
de telles applications ainsi que dans d’autres systèmes militaires, notamment les véhicules 


militaires, les sous-marins et les diodes laser de haute puissance [52]. 


1.2.6.4 Refroidissement des systèmes spatiaux 

Des recherches plus poussées sur les nanofluides conduiront certainement au développement 
de la prochaine génération de dispositifs de refroidissement incorporant des nanofluides pour des 
systèmes électroniques à très haute chaleur, offrant ainsi la possibilité de simplifier les exigences de 
refroidissement pour les applications spatiales. You et al. [59] et Vassallo et al. [60] ont effectué 
des études pour montrer que la présence des nanoparticules dans le fluide de refroidissement en 
électronique joue généralement un rôle très important dans les applications de l'espace où la densité 


de puissance est très élevée. 


1.2.6.5 Refroidissement des systèmes nucléaires 

La technologie des nanofluides peut être aussi utilisée dans le secteur de l’énergie nucléaire. 
Les nanofluides peuvent être utilisés comme fluide de refroidissement des réacteurs à eau sous 
pression et également pour les liquides de refroidissement des systèmes de refroidissement 
d'urgence du cœur des réacteurs à eau sous pression et des réacteurs à eau par ébullition. 
L'utilisation d’un nanofluide dans les accumulateurs des systèmes de refroidissement d’urgence du 
cœur et l’injection en sécurité peuvent ainsi augmenter les marges de température maximale. Les 
nanofluides peuvent aussi jouer le rôle de fluide de refroidissement pour la rétention du noyau 


fondu lors d’accidents graves dans les réacteurs à eau légère à haute puissance [52]. 


1.2.6.6 La biomédecine et l’activité antibactérienne 

Les nanofluides et les nanoparticules ont beaucoup d’applications en biomédicale. Par 
exemple, pour éviter quelques effets secondaires des méthodes traditionnelles de traitement du 
cancer, les nanoparticules à base de fer ont été employées comme transporteur des médicaments. 
Pareillement, pour une chirurgie plus sûre et pour produire un refroidissement efficace autour de la 
région chirurgicale, Jordan et al. [61] ont effectué une étude où les nanofluides pourraient être 
employés pour produire une température plus élevée dans le but de tuer les cellules cancéreuses 
sans affecter les cellules saines voisines. 

Des recherches sur l’activité antibactérienne de nanoparticules de CuO ont montré qu’elles 
possédaient une activité antibactérienne contre quatre souches bactériennes. La taille des 
nanoparticules était inférieure à celle des pores de la bactérie, ce qui lui donnait la propriété unique 


de traverser la membrane cellulaire sans aucune entrave. On pourrait émettre l’hypothèse que ces 
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nanoparticules formaient des complexes stables avec des enzymes vitales à l’intérieur des cellules, 


ce qui entraverait le fonctionnement cellulaire, entraînant ainsi leur mort [52]. 


1.2.6.7 D’autres applications 

Les nanoréfrigérants et les nanolubrifiants sont des types spéciaux de nanofluides. Ces 
nano-additifs améliorent les performances du système de réfrigération, et peut aider à augmenter 
l’efficacité du système. Les nanofluides ont potentiellement des applications dans la majorité des 
industries de transformation, telles que les matériaux, la production chimique, la nourriture, la 
boisson, le pétrole et le gaz. Dans l’industrie de l’énergie renouvelable, les nanofluides sont utilisés 
pour augmenter le transfert de chaleur à partir des capteurs solaires aux réservoirs de stockage. Le 
remplacement de l’eau de refroidissement et de chauffage par des nanofluides pourrait permettre 
d’atteindre de grandes économies d’énergie associées à d’importantes réductions des émissions 
polluantes. Le stockage de l’énergie thermique sous forme de chaleur sensible et latente est devenu 
un aspect important de la gestion de l’énergie, l’accent étant mis sur l’utilisation et la conservation 
efficaces de la chaleur résiduelle et de l’énergie solaire dans l’industrie et les bâtiments. Les 
performances thermiques des matériaux à changement de phase à base de nanofluides indiquent un 
potentiel de substitution aux matériaux à changement de phase classiques dans les applications de 


stockage à froid [52, 62]. 
L3 Caractéristiques thermophysiques des nanofluides 


Les propriétés thermophysiques des solutions telles que la masse volumique, la chaleur 
spécifique, l'expansion thermique, la viscosité et la conductivité thermique sont profondément 
modifiées par l’ajout des nanoparticules. De nombreux paramètres caractérisant ces nanoparticules 
peuvent avoir un effet notable sur les valeurs des paramètres thermophysiques du nanofluide 
obtenus (la nature des nanoparticules, leur taille, la fraction volumique 6, la conductivité du fluide 
de base et celle des nanoparticules, la température du milieu, etc.). De nombreuses recherches 
théoriques et expérimentales ont été menées, lors de ces deux dernières décennies, en vue de mieux 
caractériser les propriétés physiques des nanofluides, cependant ceci est encore loin d'être 
satisfaisant. Les résultats expérimentaux de l’amélioration de la conductivité thermique en fonction 


de ces paramètres sont résumés dans les tableaux (IL.1-IL.4) donnés en Annexe A. 


La conductivité thermique et la viscosité dynamique sont les deux propriétés 
thermophysiques des nanofluides les plus étudiées comme l’atteste le nombre de modèles 


développés pour le calcul de ces deux grandeurs. Relativement, la masse volumique et la chaleur 


spécifique des nanofluides suscitent peu d’intérêt, notamment sur le plan expérimental [52]. 
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L.3.1 Conductivité thermique des nanofluides 

La conductivité thermique, À (W/m.K), exprime la capacité d’un matériau de conduire ou de 
transmettre de la chaleur. Elle est définie comme la quantité de chaleur transférée par unité de 
surface et de temps sous un gradient de température de 1K par mètre [63]. 

La conductivité thermique est la propriété la plus importante des nanofluides car la 
principale raison de l’ajout de nanoparticules dans le fluide de base est justement d’augmenter sa 
conductivité thermique. C’est pour cela, qu’une grande partie des travaux traitant les nanofluides 
porte sur la quantification de l'augmentation de leur conductivité thermique. 

D'ailleurs, l'objectif de ces travaux est de mesurer cette conductivité thermique afin de 
choisir les nanoparticules et le fluide porteur le plus approprié et ceci pour une conductivité 
thermique la plus élevée possible. Les travaux ont également porté sur la façon de prédire 
l'augmentation de conductivité thermique par rapport à celle des nanoparticules et du fluide porteur, 
en fonction de la taille et de la forme de ces nanoparticules. Et ceci en proposant différentes 
relations prédisant le rapport de conductivité thermique du nanofluide sur celle du fluide à savoir: 
knf 

kp 

1.3.1.1 Méthodes de mesures de la conductivité thermique 

Différentes techniques de mesure de la conductivité thermique des nanofluides ont été 
utilisées, la majorité des expériences de mesure de la conductivité thermique dans les nanofluides 
sont réalisées à l'aide de la méthode du fil chaud transitoire [64-67]. 

D'autres techniques courantes sont employées comme la méthode d'oscillation de la 
température, la technique de déviation du faisceau optique [36], la méthode 3-oméga [68] et la 
technique de l'onde thermique [69]. Certaines de ces méthodes ne sont peut-être pas forcément 
valides lorsqu'elles sont appliquées aux nanofluides, car aucune d'entre elles n'a été initialement 
conçue et destinée aux nano-suspensions. Des différences significatives de valeurs de la 
conductivité thermique mesurée par ces différentes techniques existent, comme l’a signalé 


Buongiorno et al. [70]. 


1.3.1.2 Les paramètres influençant la conductivité thermique des nanofluides 

Dans la littérature, on trouve un nombre important d’études montrant que la conductivité 
thermique des nanofluides dépend de la nature du liquide de base, de la nature du matériau 
constitutif des nanoparticules, de la taille, de la forme et de la concentration volumique des 


nanoparticules, de la température, de la présence d’additifs ainsi que du pH du mélange comme le 


montre la figure L.11[5, 10, 71]. 
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Figure L.11: Paramètres affectants la conductivité thermique des nanofluides [72, 73]. 


L'influence des différents paramètres sur la conductivité thermique des nanofluides est 


présentée dans la partie I de l’ Annexe A. 


1.3.1.3 Modèles de la conductivité thermique 


S 


Etant donné que les nanoparticules à haute conductivité thermique augmentent 
remarquablement la conductivité thermique de nanofluides, plusieurs chercheurs ont développé 
différents modèles d’estimation de cette dernière. En plus de nombreuses expériences ont aussi été 
menées afin de valider ces modèles analytiques, et ceci en les comparant aux résultats obtenus 
expérimentalement. Ces modèles sont des outils de prédictions et d’appréhension des mécanismes 
de transfert thermique au sein des nanofluides [3]. 

Il existe un panel important de modèles de calcul de la conductivité thermique apparente 
dans un nanofluide; et dans ce qui suit nous présentons succinctement quelques-uns de ces modèles. 


Ces modèles recouvrent en réalité deux grands types d’approches, d’une part par analogie électrique 
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et d’autre part par résolution exacte de l’équation de chaleur. En fait, seule l’expérimentation 
physique permet d'obtenir une évaluation précise de la conductivité des nanofluides [7]. 

a) Modèles de type Maxwell 

Ÿ”_ Modèle de Maxwell 

Maxwell [32] a été l’un des premiers à étudier analytiquement la conduction thermique dans les 
particules en suspension en ignorant les effets d’interactions entre ces particules. L’équation de 
Maxwell est obtenue en résolvant l’équation de Laplace pour le champ de température, et elle est 


donnée par : 
Rs — kf 


ks +2kf —p(ks kr) f (L3) 


knp=ks + 3@ 


Notons que le modèle de Maxwell est satisfaisant pour les suspensions contenant des particules 
sphériques de faibles concentrations volumiques, toutefois 1l ne tient pas compte de l'effet de la 
taille n1 de la forme des particules. 

Ÿ_ Modèle de Hamilton-Crosser 

Le modèle de Hamilton-Crosser [74] est une extension du modèle de Maxwell. II concerne les 
particules de forme quelconque. Pour prendre en compte cette diversité de forme, Hamilton 
introduit un facteur géométrique appelé sphéricité (y) qui fait référence à la géométrie sphérique. 
Ce coefficient est défini comme étant le rapport de la surface de la sphère ayant le même volume 
que les nanoparticules et la surface d’une nanoparticule. 
La conductivité thermique apparente du milieu est donnée par l’expression suivante : 


Kknf _ ks +(n-1)kf-(n-1)(kr = ks )@ a à) 
kp  Ks+(n-1)k+(kp_ ks )p ‘ 


Où n est un facteur de forme empirique donné par : 
n = - , (n =3 pour les particules sphériques et n = 6 pour les particules cylindriques. 


Pour y =1 (particules sphériques), le modèle de Hamilton-Crosser est identique au modèle de 
Maxwell. 

b) Modèles des nano-couches 

Y” Modèle de Yu et Choi 

Une autre expression pour calculer la conductivité thermique a été introduite par Yu et Choi 
[75]. Ils ont proposé de modéliser les nanofluides comme étant un liquide de base et des particules 
solides séparés par une couche nanométrique, cette couche agit comme un pont thermique entre le 
fluide et les nanoparticules. De plus, ils ont supposé que la conductivité thermique de la couche 
nanométrique est plus grande que celle du liquide. Et ils ont proposé la relation suivante : 


Knf _ Ks+2kp+2(ks- kg )A+B)5e (L5) 


kf ks +2kf-(ks - kf )A+B)3@ 
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Où B est le rapport de l’épaisseur de la couche nanométrique sur le rayon des particules [3]. 
Ÿ” Modèle de Bruggeman 
Le modèle proposé par Bruggeman [76] semble mieux estimer certains résultats expérimentaux 


comparé aux autres modèles [77], surtout pour le cas des nanoparticules sphériques de différentes 


concentrations. 
k 
= = [Ge — 1)ks + (2—3p)kr] +LVA (L.6) 
Avec : 
2 
A= | — 1} (£) + (2 — 36)? + 2(2 + 99 — 9?) (x) (1.7) 


c) Modèles basés sur le mouvement Brownien 
Le mouvement Brownien est une description mathématique du mouvement aléatoire d'une 
particule de taille importante immergée dans un fluide et qui n'est soumise à aucune autre 
interaction qu'avec les petites molécules de ce fluide. Lorsqu'on considère un nanofluide, on a un 
transport d'énergie directement par les nanoparticules ainsi que par la micro-convection qui est due 
au mélange du fluide autour de ces nanoparticules. 
Plusieurs modèles ont été développés pour la détermination de la conductivité thermique des 
nanofluides basée sur le mouvement Brownien. Deux de ces modèles sont cités ci-dessous. 
Ÿ” Modèle de Choi et Jang 
Choi et Jang [78] ont modélisé la conductivité thermique des nanofluides tenant en compte 
l'effet du mouvement Brownien des particules. Le modèle proposé ne dépend pas seulement de la 
conductivité des nanofluides, mais aussi de la taille et de la température des nanoparticules. Le 
transport d'énergie a été considéré comme étant composé de quatre types de transport : conduction 
thermique du fluide de base, conduction de chaleur dans les nanoparticules, collisions entre les 
nanoparticules (due au mouvement Brownien) et micro-convection causée par le mouvement 
aléatoire des nanoparticules. Toutefois, les collisions entre les particules sont considérées comme 
négligeables par rapport aux autres formes de transport. 


Selon cette étude l'expression suivante a été proposée : 
: d 
kas = kr — p) + k5p + 3C a kr ReäPr;p (LS) 
Où C; est une constante de proportionnalité, d; le diamètre de la molécule de fluide, d,, le diamètre 
des nanoparticules, Pr; le nombre de Prandtl du fluide de base. 
k;, est définie comme : 


ki =Bk} (19) 
Où B est le nombre de Reynolds défini comme une constante: 
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CRM.dp 


(1.10) 
Vf 


Re; —= 


Avec CR. Comme la vitesse du mouvement aléatoire des nanoparticules et v; la viscosité 
cinématique du mouvement aléatoire du fluide de base. Cr my. , qui peut être déterminée de la façon 


suivante : 


Où 1; est le parcours moyen libre des molécules de fluide de base et D, est le coefficient de 
diffusion de la particule qui peut être calculé par l'expression d'Einstein suivante [78]: 


__ kgBT 
Un 37Ufdp 


(1.12) 


Do 


Où k3 est la constante de Boltzmann avec T en Kelvins, et y; la viscosite dynamique du fluide de 
base [78]. 
Y_ Modèle de Koo et Kleïnstreuer 


Koo et Kleinstreuer [79] ont considéré que la conductivité thermique des nanofluides est 


composée de deux parties : 

ka = Kstatistique + Kkprownien (1.13) 
Où : Kstatistique est l'augmentation de la conductivité thermique due aux nanoparticules 
et Kprownien tient compte de l'effet du mouvement brownien. Pour la partie statique on a la formule 


de Maxwell et en conséquence, l'expression suivante a été proposée : 


_ keT = 
Kkprownien = S: 10*BppsCp,f as PP: T) (L.14) 
Où p,, pf sont la densité de la particule et du nanofluide respectivement, T la température en 


Kelvin et c,,f est la capacité calorifique du fluide de base. Un élément additionnel B a été introduit 
pour tenir compte de l'effet du mouvement des particules. En général, B est lié au mouvement des 
particules, il dépend donc de la fraction volumétrique mais aussi de la température, de la forme et 
des propriétés de matériaux (Koo et Kleinstreuer [79]). 

Koo et Kleinstreuer, supposent que B dépend seulement de la concentration volumétrique et 
que tout autre dépendance est incluse dans le paramètre F'; cela signifie que le terme Ê devient plus 
effectif avec l'augmentation de la fraction volumétrique. Le paramètre F qui a été obtenu en 
utilisant les résultats de Das et al. [36] pour le nanofluide à base de CuO, a été introduit afin 


d'inclure les autres dépendances du modèle. Ces deux paramètres sont déterminés 


expérimentalement. Dans le travail de Koo et Kleinstreuer F a la forme : 
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F = (-134.63 + 1722.39) + (0.4705 — 6.049)T (1.15) 
d) Modèles basées sur les grappes 
Les nanoparticules sont connues pour leurs tendances à former des grappes (Prasher et al 
[38]). Une approche pour modéliser ces grappes est donnée par la théorie des fractals (Wang et al 
[47)), celle-ci suppose que la formation des grappes est due au transport rapide de la chaleur sur une 
distance (relativement grande) car la chaleur est conduite plus rapidement par des particules solides 
en comparaison au liquide. 
Le paragraphe suivant présente un modèle théorétique fondé sur les grappes de 
nanoparticules. Le modèle considère également l'effet du mouvement brownien. 
Ÿ_ Modèle de Xuan et al 
Xuan et al. [80] ont proposé l'expression suivante de la conductivité thermique, basé le 
mouvement brownien ainsi que sur les grappes des nanoparticules : 


pe, Fo 2Kf F2 Kg )9 | PrECpp [KT (16) 
T kp+2kf-(Kp- kf JP 2kp A3rraurs ‘ 


knf 

Ici, r est le rayon apparent des grappes de nanoparticules, qui est déterminé 
expérimentalement, T la température en Kelvin, a; la viscosité dynamique du fluide de base 
(comme calculée par Li et Xuan [42). Le premier terme de l'expression de droite correspond à la 
formule de Maxwell pour la conductivité thermique et le deuxième ajoute l'effet du mouvement 
aléatoire des nanoparticules. 

Afin d’examiner la contribution de ce dernier terme, des valeurs pour des nanoparticules de Cu 
de 50 nm avec un @ = 0:03, cette contribution est de 11% lorsque les grappes se forment et de 17% 
quand les grappes ne se forment pas. Dans le cas où @ = 0:04, les contributions sont de 14% et 
24% respectivement. Il a été aussi indiqué que le mouvement Brownien des nanoparticules devient 
plus efficace avec l'amélioration de la température. D'autre part, au fur et à mesure que les 
nanoparticules et les grappes deviennent plus grandes, le mouvement aléatoire devient plus lent et 
cela diminue l'amélioration de la conductivité thermique [10]. 

Ÿ_ Modèle de Chen et al 

Un autre modèle qui considéré l'effet des grappes est le modèle proposé par Chen et al. [81]. La 
formule proposée pour le calcul la conductivité thermique est une modification de la formule de 
Hamilton et Crosser [74] : 


at-1)kr-(n-1)(kr-ka)P ci 


k 
kn _ #f kat(n-Dks+(ks-ka)pe ui 
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Où k et @,, sont la conductivité thermique et la fraction volumétrique des grappes, 


respectivement. Les valeurs de n sont égales à 3 pour les particules sphériques et 5 pour les 
particules cylindriques. 


: 3—-D 
Pi = P (=) (1.18) 


Th 


Avec TA et7, sont les rayons de grappes et de nanoparticules respectivement. D correspond 
à l'indice fractal, avec une valeur prise de 1,8. Les valeurs 7, et 7, sont égales à 2,75 et 3,34 pour 


le (TiO/Eau) et (TiO-/Ethylène Glycol) de forme sphérique. L'estimation de k,, est plus complexe 
puisqu’il dépend de la forme géométrique des particules. Bruggeman [76] a calculé l'expression 
pour les particules sphériques pendant que Nan et al. [5] l'ont calculé pour les particules 


cylindriques. 


e) Les nanofluides contenant des nanotubes de carbone 
Généralement, les particules des nanotubes de carbone peuvent être considérées comme des 
sphéroïdes allongées qui satisfont a >>b-c. Par conséquent, de nombreux modèles développés pour 
les suspensions sphéroïdes allongées peuvent être utilisées pour estimer la conductivité thermique 
effective des nanofluides à base de nanotubes de carbone. Grâce à cette approche, Nan et al. [82], 


ont proposé une équation approximative obtenue à partir de l’équation de Fricke [83,84] et ceci en 


utilisant la condition que kR est beaucoup plus grande que K, : 


1 ks 
Ress me AU T7 (19) 
Nan et al. [85] dans leurs travaux ont négligés l’effet de la résistance de contact : 


1.3.1.4 Résultats expérimentaux et modèles de la conductivité thermique 


La figure I.12 présente une comparaison des résultats de k,;/kf (rapport de conductivité 
thermique du nanofluide sur celle du fluide) en fonction de @ pour des nanofluides (ALO3/Eau), 
obtenus par différents modèles. Il ressort de cette comparaison que le modèle Hamilton-Crosser 
donne des prédictions plus élevées. Tandis que celles du modèle de Yu-Choi se situent entre celles 


du modèle de Hamilton-Crosser et de Maxwell. 


La figure 1.13 présente le rapport conductivité thermique k,/k- en fonction de @ pour 
des nanofluides (ALO-Eau), obtenus à partir de différents travaux expérimentaux. Ils sont 


comparés aux résultats des modèles de Hamilton-Crosser et Yu-Choi et ceci pour une valeur de 
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(y = 0.7). On remarque que les résultats obtenus par ces modèles sont loin de ceux obtenus 
expérimentalement par Massuda et al. [33] et Eastman et al. [8]. Par contre, ils sont en bon accord 


avec ceux obtenus par Xie et al. [37] et Lee et al. [35] et ceci pour les faibles valeurs de o. 


La figure I.14 présente la conductivité thermique du nanofluide (CuO/Eau) prédite par le 
modèle de (Hamilton-Crosser) comparée aux résultats expérimentaux trouvés dans la littérature. 
On constate que les résultats du modèle de Hamilton-Crosser sont en très bon accord avec ceux 
obtenus par Lee et al. [35] ainsi que ceux de Das et al. [36]. Par contre ils diffèrent sensiblement 


de ceux mesurés par Eastmen et al. [9]. 


Maxwell ALOs+eau 


----YUetChoi 
— — Bruggeman 


Figure I.12: Comparaison des résultats prédits du modèle de H-C avec ceux d’autres modèles pour 


les nanofluide (A1:03/Eau) [7]. 


H-C (W=0.7) ALO;3+eau 
=  Xieetal. 
=== =Yu-Choi 
+ Massudaetal. 
e  Leeetal. 
À Eastmanetal. 


Figure L.13: Comparaison des résultats expérimentaux de nanofluides (A1:03/Eau) et ceux obtenus 


à partir des modèles de H-C et de Yu-Choi [7]. 
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——H-C 
A Leeetal. CuO+eau 


= Dasetal. 
e Eastman et al. 


Figure IL.14: Comparaison entre les résultats prédits du modèle de H-C avec ceux obtenus 


expérimentalement pour les nanofluides (CuO/Eau) [7]. 


La figure L.15 montre l’influence de la fraction volumique des nanoparticules sur la 
conductivité thermique du nanofluide (Cu/Eau). On remarque que les résultats prédits par le 
modèle de H-C sont satisfaisants (si l'on prend une sphéricité y = 0.3), par contre pour w=0.7 


les résultats sont loin de l'être. 


H-C (W=0.7) Cu+eau 
+ Eastmanetal. 
— H-C (W=0.3) 


Figure L.15: Comparaison entre les résultats du modèle de H-C avec des résultats expérimentaux 


pour les nanofluides (Cu/Eau) et pour deux valeurs de w [7]. 


Les suspensions à base de nanotubes de carbone donnent de plus grandes conductivités 
thermiques comparées à celles des suspensions de nanoparticules métalliques. Ceci s’explique 
évidement par les valeurs élevées de la conductivité thermique de ces nanotubes de carbones. 

La conductivité thermique des nanofluides (Nanotubes de Carbone/Eau) augmente avec les 
fractions volumiques et ceci pour @ variant de 0% à 8%. En utilisant comme conductivité thermique 
du nanotube de carbone (k, = 2500 W.m!. K'!), on obtient pour une fraction volumique de 8% une 


amélioration de la conductivité thermique de 74% figure I.16. 
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—H-C 


—— Maxwell NTC+eau 
——-Yu-Choi 
Bruggeman 


Figure I.16: Comparaison des résultats prédits par le modèle de H-C et des résultats obtenus par 


d’autres modèles pour les nanofluides (NTC/Eau) [7]. 


—H-c NTC+eau 


A Xieet al. 
#* Assaelat al. 


Figure L.17: Comparaison des résultats prédits par le modèle de H-C et des résultats obtenus 


expérimentalement pour des nanofluides (NTC/Eau) [7]. 


Sur la figure [.17, sont présentés les valeurs du rapport k,;/k;obtenues à partir du 
modèle de Hamilton-Crosser et celles obtenues des résultats expérimentaux d’Assael et al (avec o 
entre 0% à 6%) [51] et de Xie et al (avec @ entre 0% et 2%) [48]. Les conductivités thermiques 
obtenues pour les nanofluides (NTC/Eau) sont plus grandes que celles obtenues des nanofluides 
avec des nanoparticules sphériques. Pour les nanofluides à base de nanotubes de carbones le 
modèle de Hamilton-Crosser n'estime de manière satisfaisante la conductivité thermique que pour 
les faibles valeurs de o. 

Dans le tableau 1.4, sont présentées les valeurs calculées du rapport k,/k, à partir du 
modèle H-C en fonction de différents p% pour (NTC/Eau) et (NTC/EG). On constate qu’il y a une 
amélioration de 54% de la conductivité thermique pour une fraction volumique de 6 % pour les 


deux nanofluides (NTC/Eau) et (NTC/EG). On remarque aussi que cette augmentation est quasi 
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linéaire avec la fraction volumique. Assael et al. [51] ont remarqué qu’il y a une amélioration de la 
conductivité thermique du nanofluide (NTC/Eau) de 28% pour une fraction volumique de 6%, alors 
que Choi et al. [74] ont trouvé une amélioration de 150% pour une fraction volumique de 1% dans 


les nanofluides (NTC/EG). 


Table [.4: Valeurs calculées de (k,/k) à partir du modèle H-C en fonction de 6 % pour deux 
nanofluides (NTC/Eau) et (NTC/EG) [7]. 


(knç/ks) 
po % H-C (NTC+eau) H-C (NTC+EG) 

0 1 1 

1 1.0864 1.0865 
» 1.1746 1.1748 
3 1.2645 1.2649 
4 1.3564 1.3568 
5 1.4502 1.4507 
6 1.5459 1.5466 
7 1.6438 1.6446 
8 1.7437 1.7446 


Le tableau LS ci-dessous, montre les différentes valeurs calculées par le modèle Hamilton- 
Crosser du rapport (k,;/k;) en fonction de @ variant de 0 à 8%, et ceci pour cinq nanofluides à 


savoir (Al,03/Eau), (CuO/Eau), (Cu/Eau), (Ag/Eau) et (Au/Eau). 


Table LS: Valeurs de (k,f/kf) calculé à partir du modèle de H-C en fonction de , pour différents 


nanofluides et avec l’eau comme fluide porteur. 


Rapport des conductivités thermiques (kng /k f) 


Ag 1 1.0430 1.0869 1.1317 | 1.1774 | 1.2241 1.2718 1.3205 1.3702 


NTC 1 1.0864 1.1746 1.2645 | 1.3564 | 1.4502 1.5459 1.6438 1.7437 


L.3.2 La viscosité dynamique 
La viscosité caractérise l’aptitude d’un fluide à s’écouler. La viscosité est une propriété 
importante de l’écoulement des fluides. Sa connaissance est fondamentale pour toutes les 


applications impliquant le transport de fluides, en plus le transfert de chaleur par convection dépend 
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directement de la viscosité des fluides. Il a été démontré précédemment que l’ajout de 
nanoparticules dans un fluide de base permet d’améliorer la conductivité thermique du nanofluide 
obtenu, mais cela peut aussi conduire à une augmentation défavorable de la viscosité dynamique 
[S21: 

La majorité des études traitant de la viscosité des nanofluides ont eu recours, pour son 
calcul : de la viscosité du fluide de base, de la fraction volumique 6 et aux relations de Stokes- 
Einstein ou de Brinkman. Ces relations ont été établies dans des cas précis, à savoir pour des 
particules métalliques de forme sphérique, en faible concentration, parfaitement dispersées et sans 
interaction inter-particules. Or les nanoparticules utilisées dans les différentes expérimentations ne 
sont ni sphériques ni de même taille, en particulier lorsqu'il s'agit de nanofluides avec des nanotubes 
de carbone. 

La dispersion de nanoparticules dans un fluide de base modifie notablement sa viscosité, 11 
s'ensuit alors une modification de la structure des écoulements convectifs qui dépend fortement de 
la viscosité. Pour les suspensions constituées de nanoparticules sphériques, on dispose de différents 
modèles et de formules semi-empiriques décrivant l’évolution de la viscosité relative en fonction de 


la fraction volumique. 


La viscosité relative est le rapport de la viscosité du nanofluide sur celle du fluide de base … [71]. 
f 


1.3.2.1 Les Paramètres affectants la viscosité dynamique des nanofluides 

L’ajout de particules solides dans un liquide provoque une augmentation, plus ou moins 
importante, de la viscosité du mélange obtenu. Là encore, de nombreux facteurs peuvent être à 
l’origine des disparités constatées. Pour des nanoparticules identiques (nature et taille) dans un 
même fluide de base ; la qualité de la dispersion, la viscosité du fluide de base, la fraction 
volumique des nanoparticules et la température peuvent agir considérablement sur la viscosité, 


comme le montre la figure I[.18. 
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Figure L.18: Paramètres affectant la viscosité des nanofluides [72]. 


En augmentant la fraction volumique des nanoparticules dans un fluide porteur, on 
augmente sa conductivité thermique ce qui est intéressant pour le transfert de chaleur, mais cela 
augmente nécessairement la viscosité du nanofluide obtenu. Lorsque la viscosité dynamique 
augmente, la capacité du fluide à s’écouler diminue. La viscosité des liquides tend généralement à 
diminuer lorsque la température de ce dernier augmente. Des corrélations semi-empiriques basées 
sur des résultats expérimentaux sont utilisées pour le calcul de la viscosité dynamique [86, 87]. 

1.3.2.2 Les différents modèles de viscosités 

De nombreuses études théoriques et expérimentales ont été menées afin d’estimer ou de 
mesurer la viscosité dynamique des nanofluides, et ceci pour quantifier l'importance d’un certain 
nombre de paramètres liés à la nature des nanoparticules, leur taille, la température, etc. Plusieurs 
modèles théoriques ont été développés pour prédire l’évolution de la viscosité dynamique des 
suspensions sous certaines conditions. Dans ce qui suit, nous nous intéresserons aux modèles les 
plus utilisés dans le cas des nanofluides [3]. 

a) Modèle d’Einstein 

L'historique des études sur la viscosité effective des mélanges de particules-liquide est 
presque aussi long que celui de la conductivité thermique effective. Einstein [87] fut le premier à 
calculer la viscosité effective d’une suspension de sphères en se basant sur les équations 


hydrodynamiques. Il a déterminé l’énergie dissipée autour d’une seule particule dans un fluide en 
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associant son énergie au travail réalisé pour la faire avancer relativement au fluide. A partir de ce 
résultat, 1l a trouvé une expression décrivant une dépendance linéaire de la viscosité des suspensions 
dans un fluide avec leur concentration : 

Un = FO + 2.59) (1.21) 
Où y» et 1 dénotent respectivement la viscosité du nanofluide et la viscosité du fluide de base. 

@ est la fraction volumique des nanoparticules. 

La formule d’Einstein a depuis été vérifiée expérimentalement et est considérée satisfaisante 
pour des suspensions très diluées de petites particules sphériques (Browniennes ou non 
Browniennes). Cette relation n’est valable que pour des concentrations volumiques inférieures à 
1%. Il est à noter que ce modèle ne tient pas compte des effets de la taille des particules et des 
interactions inter-particules [3,7]. 

b) Modèle de Brinkman 
La formule de Brinkman [89] complète le modèle d’Einstein jusqu’à une concentration 


volumique inférieure à 4%, et elle est favorablement acceptée par de nombreux chercheurs. 


—, hf 
Pa — (1-p)25 (1.22) 


La formule de Brinkman décrit une évolution non linéaire de la viscosité dynamique avec la 


concentration volumique mais ne tient pas compte de la collision entre les particules [3, 52, 7]. 


c) Modèle de Batchelor 

En prenant en compte l’effet du mouvement Brownien et en considérant une suspension 
isotrope de nanoparticules solides et sphériques, Batchelor [90] a développé l’équation de viscosité 
d’Einstein. Il a montré aussi que la viscosité n’est pas une fonction linéaire de @, mais une fonction 
non linéaire de @ avec la présence d'un terme en @?. 
Son modèle est donné par : 
Uns = Up + np + ku p°) (1.23) 
Où # est la viscosité intrinsèque et k;, est le coefficient de Huggins. La valeur de 7 et k,, est de 2.5 
et 6.5 respectivement pour des particules sphériques. Cette équation permet le calcul de la viscosité 


de suspension de nanoparticules avec @ >10% [3, 52, 7]. 


d) Modèle Maiga et al 
D’autres relations ont été proposées dans la littérature qui se limite à des applications bien 
précises. Maiga et al. [18] ont proposé une corrélation à partir de résultats de mesures, à savoir : 
Uns = Uf(123p° + 7.3 + 1) (1.24) 
En comparant cette corrélation aux modelés d’Einstein et de Brinkman, on trouve que ces deux 


dernières approches sous-estiment les valeurs expérimentales de la viscosité des nanofluides. 
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e) Modèle Pack et Cho 

Pack et Cho [40] ont développé un modèle de calcul de la viscosité des nanofluides basé sur 
la fraction volumique des nanoparticules AO; disperses dans l’eau en prenant la température 
ambiante comme référence. Ils ont constaté que la viscosité du nanofluide diminue lorsque la 
température augmente. Ce modèle est donné par : 

Uns = (533.9? + 39.11 + 1) (1.15) (1.25) 
1.3.2.3 Viscosité dynamique des nanofluides 

La compréhension des problèmes de convection est strictement liée à la viscosité des 
nanofluides [2]. Et les données expérimentales de la viscosité effective sont limitées à certains 
nanofluides, tels que : 

o  AbO3/Eau ; (Pak et Cho [40], Das et al. [36], Putra et al. [91], Li et al. [92], Heris et al. [931). 
o  CuO/Eau ; (Heris et al [93], Kulkarni et al. [94]). 

Oo TiO,/Eau ; (Pak et Cho [40]). 

o MWCNT/Eau ; (Ding et al. [501]). 

o AbO;/ Octane ; (Liu et al. [66]). 

Afin de minimiser la divergence entre les mesures expérimentales et les calculs, les 
chercheurs ont proposé des équations pour des applications spécifiques telle que : A1:03/Eau et 
ALO;/Ethylène glycol (Maïga et al. [181), Ti0,/Eau (Tseng et Lin [95]) et CuO/Eau (Kulkarni et al. 
[38]). Le problème avec ces équations, c’est qu’elles ne se réduisent pas à l’équation d’Einstein 
pour de très faibles concentrations volumique, et par conséquent elles manquent de base physique. 
Nguyen et al. [96] ont montré que la viscosité augmente avec la concentration de particules et 
décroit avec la température. Ils ont remarqué aussi que pour une fraction volumique de particules 
inférieure à 4%, la taille n’a quasiment pas d’influence sur la viscosité. Alors que pour des fractions 
supérieures, la taille des particules joue un rôle important. Prasher et al. [38] ont aussi étudié 
l’influence de la taille des particules sur la viscosité dynamique pour différents nanofluides et à des 
concentrations identiques ; et selon les courbes obtenues rien ne permet de tirer de conclusions 
probantes [71]. 

Pak et Cho [40] se sont également penchés sur la variation de la viscosité de deux 
nanofluides et ont montré que leur comportement ne peut se prédire par des modèles empiriques 
standards pour les viscosités de suspension, comme celui développé par Das et al. [97]. Lee et al. 
[35] ont rapporté une augmentation substantielle de la viscosité pour l'eau et les nanofluides à base 


d'éthylène glycol pour les nanoparticules d’Al:0:. Ils ont également évoqué que cette augmentation 


pourrait compenser le transfert de chaleur observé dans les nanofluides [2]. 
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Das et al. [36] ont mené une étude similaire pour un nanofluide A1:0;/Eau, ils ont confirmé 

en particulier, le comportement Newtonien du nanofluide pour des concentrations de 1% à 4%. Ils 
ont également montré une tendance linéaire de la viscosité par rapport au taux de cisaillement, 
même si les valeurs de viscosité sont encore plus élevées que celles de l'eau pure [2]. 
Il existe des divergences dans la littérature sur la nature de la rhéologie des nanofluides ; certains 
auteurs croient que les nanofluides sont Newtoniens [34, 38] alors que d’autres ont vu que les 
nanofluides sont non Newtoniens [50]. Dans ce contexte, Chen et al. [81] ont réalisé des études à la 
fois expérimentales et théoriques de l'éthylène-glycol avec des nanoparticules de TiO, pour des 
concentrations de 0,5% à 8,0% en poids et à une gamme de température variant de 20 à 60 °C. 
Leurs résultats ont montré que dans ces conditions le comportement du nanofluide est Newtonien 
(21: 

Sergis et Hardalupas [98] ont également conduit des analyses statistiques de convection de 
chaleur des nanofluides ; et ils ont trouvé que le mode de convection semble être augmenté avec 
l'augmentation de la température, de la concentration volumétrique des nanoparticules et de la 
dimension décroissante des particules. Malheureusement, les expériences sur la convection couvrent 
seulement des nanoparticules d'oxyde avec des concentrations élevées de nanoparticule (> 1% du 
volume). En plus, les expériences de convection sont aussi nécessaires pour pallier à ces lacunes, 
particulièrement de particules métalliques (qui ont une conductivité thermique beaucoup plus 


élevée) et pour de faibles concentrations (< 1% du volume) [99]. 


L.3.3 La masse volumique 

La masse volumique d’une substance, notée p, correspond à la masse de cette substance 
dans une unité de volume. La masse volumique des nanofluides est proportionnelle à la fraction 
volumique des particules solides. Le calcul de la masse volumique effective p, d’un nanofluide est 
assez simple et assez conforme aux valeurs expérimentales. Elle est souvent estimée sur la base du 
principe de la règle de mélange en supposant le nanofluide comme un milieu homogène (bonne 
dispersion des nanoparticules dans le fluide) [40, 1001. 

La masse volumique du nanofluide est donnée par la formule suivante : 
Pnf = À — P)pf + PPp (1.26) 
Où Ph, Ph et pF dénotent respectivement la masse volumique du nanofluide, la masse volumique 
des nanoparticules et la masse volumique du fluide de base. 

Les valeurs des rapports de la masse volumique des nanofluides sur la masse volumique du 
fluide de base pour cinq nanofluides (eau+A1,03), (eau+CuO), (eau+Cu), (eau+Ag) et (eau+Au) en 


fonction de la fraction volumique (@) qui varie entre 0% et 8% ; sont présentées dans le tableau L.6. 
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Table L.6: Rapport de masses volumiques des différents nanofluides en fonction de la fraction 


volumique [7] 


Pn/Pf 

p%| Al203/Eau CuO/Eau Cu/Eau Ag/Eau Au/Eau 
0% 1 1 1 1 1 
1% 1.0296 1.0546 1.0797 1.0951 1.1833 
2% 1.0593 1.1092 1.1594 1.1903 1.3666 
3% 1.0890 1.1638 1.2391 1.2855 1.5500 
4% 1.1186 1.2184 1.3188 1.3807 1.7333 
5% 1.1483 1.2730 1.3985 1.4759 1.909167 
6% 1.1780 1.3276 14782 1.5711 2.1000 
7% 1.2077 1.3823 1.5579 1.6663 2.2834 
8% 1.2373 1.4369 1.6376 1.7615 2.4667 


L.3.4 La chaleur spécifique 

La chaleur spécifique C, est la capacité calorifique par unité de masse d’une substance ou 
d’un système homogène. La chaleur spécifique correspond à l’apport de chaleur nécessaire pour 
élever la température de l’unité de substance de 1K lors de la transformation considérée. Elle reflète 
la capacité d’un matériau à accumuler de l’énergie sous forme thermique, pour une masse donnée, 
quand sa température croit. Une capacité thermique élevée signifie qu’une grande quantité 
d’énergie peut être stockée moyennant une augmentation relativement faible de la température. La 
dispersion de nanoparticules solides aux fluides de base provoque généralement une diminution de 
la chaleur spécifique du nanofluide obtenu, car la chaleur spécifique des nanoparticules solides est 
plus faible que celle des fluides de base [3, 52]. 

Il existe deux modèles pour le calcul de la chaleur spécifique des nanofluides. La chaleur 
spécifique peut être estimée à l’aide de la règle des mélanges ou sur la base de l’hypothèse de 
l'équilibre thermique neutre des particules solides et du fluide de base. 

En se basant sur la règle de mélange d’une suspension homogène, Pak et Cho [40] ont 
proposé l’expression suivante pour le calcul de la chaleur spécifique de nanofluides : 

Cheff = PCpp+ A —-p)Cf (1.27) 
OùC, eff; Cppet Cp fpdénotent respectivement la chaleur spécifique du nanofluide, la chaleur 
spécifique des nanoparticules et la chaleur spécifique du fluide de base. Et @ est la fraction 


volumique des nanoparticules. 
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De son côté, Xuan et Roetzel [101] ont développé un modèle de la chaleur spécifique des 
nanofluides basé sur l’hypothèse de l’équilibre thermique entre les particules solides et le fluide de 
base de la manière suivante : 


PPpCpp+(i-p)PpfbCp fb 
_—_— 1.28 
peff PPp+(1-p)pfp CESR) 


Où p}, pfn dénotent respectivement la masse volumique des nanoparticules et la masse volumique 
du fluide de base. 
Cette expression a été jugée mieux adaptée pour le calcul de la chaleur spécifique des 


nanofluides par certains chercheurs [52]. 
L.3.5 Coefficient d'expansion thermique 


La variation de masse volumique sous l'action de la température est caractérisée par 


le coefficient d'expansion thermique Br = — n (5) » appelé également coefficient de dilatabilité. 


Pour calculer la valeur de ce coefficient pour les nanofluides de très nombreux auteurs ont 
proposé par analogie avec la relation [7, 52]: 


Pnf = À — P)pf + PPp (1.29) 


L’expression suivante : 


Bas = (À — p)B + pr (1.30) 


Où BF et B, désignent respectivement les coefficients de dilatabilité du fluide porteur et celle des 


nanoparticules. 


L.4 Conclusion 


Dans ce chapitre, nous avons d'abord donné une description générale des nanofluides 
(synthèse, caractérisation, production et applications). Une revue bibliographique détaillée de leurs 
propriétés thermophysiques et de leurs performances thermiques et hydrodynamiques a ensuite été 
abordée. Ainsi, nous avons montré à travers cette analyse que les performances des nanofluides 
dépendent fortement de nombreux paramètres tels que la structure et la qualité des nanoparticules 


ainsi que le type du fluide de base. 


Par conséquent, 1l est nécessaire d'évaluer rigoureusement les propriétés thermophysiques 
des nanofluides afin de mieux comprendre leur comportement d'écoulement et leurs propriétés 


thermiques en vue d'une éventuelle application dans les systèmes thermiques. 
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D’après cette étude bibliographique, il est difficile de comparer les différents résultats 
expérimentaux des nanofluides car chaque expérience a été menée pour des nanoparticules de 
différentes tailles et avec de préparations d’échantillons distinctes. Or la taille des particules et la 
préparation des échantillons sont des facteurs déterminants qui influent fortement sur les résultats 
obtenus. En plus L’agglomération des particules ne nous permet pas d’étudier de manière 
rigoureuse l’influence de la taille de ces dernières sur l’augmentation de conductivité thermique du 
fluide. 

D’autres approches théoriques ont tenté de comprendre l'origine et l'importance de cette 
augmentation de la conductivité thermique des nanofluides. Néanmoins aucune d’entre elles n’est 
parvenue à décrire convenablement tous les comportements obtenus expérimentalement. 

Un meilleur contrôle des propriétés physico chimiques de ces fluides serait donc nécessaire pour 


pouvoir les étudier de façon précise et ainsi comprendre leur comportement. 
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Chapitre II 


Approche mathématique | ? he 


Introduction 

L’expérimentateur, quel que soit son domaine d’étude, est toujours confronté au problème 
difficile de l’organisation optimale de ses essais. C'est-à-dire comment obtenir les bonnes 
informations dans les meilleurs délais et pour le moindre coût. 

Les scientifiques ont proposé différentes méthodes qui permettent l’optimisation du choix 
des essais et de celui de leur enchaînement en cours de l’expérimentation. Celles-ci sont plus 
efficaces que les démarches habituelles de conduite des essais, permettant ainsi d’arriver aux 
mêmes résultats qu'avec trois à quatre fois moins d’essais. 

Ce chapitre examine l’approche mathématique et la modélisation des transferts de chaleur de 
nanofluides dans le but de développer un nouveau modèle mathématique de la conductivité 
thermique des nanofluides en incluant les paramètres les plus influents. Ce modèle est généré par la 
méthode des plans d’expérience présentée dans la première partie de ce chapitre. 

En ce qui concerne la deuxième partie, elle présente la technique des Réseaux de Neurones 
Artificiels (RNA). Celle-ci est une technique de traitement de données, ayant des fonctionnalités 
d’identification, de contrôle ou de filtrage aboutissant à des solutions plus efficaces. 

L'objectif de ce chapitre est de préciser les notions qui permettent d’évaluer et d’améliorer 
l’aptitude du modèle postulé à représenter les résultats expérimentaux. En un mot, ces notions nous 
aideront à trouver le modèle adéquat, qui représente le mieux possible les données expérimentales 


et qui devrait nous permettre de faire des prévisions de qualité. 
IL.1 Théorie des plans d’expériences 


IL.1.1 Introduction 

La théorie des plans d’expériences optimaux est une branche importante des statistiques, à 
l'interface avec l’optimisation, qui a de nombreux champs d’applications. Son but est de trouver les 
valeurs qu’un expérimentateur doit donner aux variables de contrôle des expériences à sa 
disposition, avant de les réaliser. Ces variables de contrôle peuvent prendre différentes formes 
(nombre de fois qu’on va réaliser une expérience, taux d’échantillonnage d’un appareil de mesure, 
temps pendant lequel on enregistre des résultats, etc.), et affectent les données mesurées. 
L’estimation que l’expérimentateur fait des quantités qu’il souhaite mesurer va donc dépendre de 


ces variables [1]. 


IL.1.2 Domaine d’application des plans d’expériences 
La méthodologie des plans d’expériences est destinée à optimiser la conduite et 
l'interprétation de l’expérimentation en faisant varier simultanément les différents facteurs et en 


fournissant des outils d’exploitation des résultats. C’est donc une approche rationnelle pour 
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construire la série d’essais la plus légère possible en prenant en compte les contraintes telles que 


celle qui relèvent du coût ou du temps par exemple [2]. 


Les plans d'expériences sont utilisés dans les études industrielles en recherche et développement, vu 
qu’ils s’inscrivent dans une démarche générale d’amélioration de la qualité. Ils interviennent dans 
de nombreux domaines industriels tels que : 
Ÿ”_ Industries chimiques, pétrochimiques et pharmaceutiques ; 
Ÿ”_ Industries mécaniques et automobiles ; 
Ÿ”_ Industries métallurgiques. 
Leur utilisation vise à la : 
Ÿ” Détermination des facteurs clés dans la conception d’un nouveau produit ou d’un nouveau 
procédé ; 
Ÿ”_ Optimisation des réglages d’un procédé de fabrication ou d’un d’appareil de mesure ; 
Ÿ”_ Prédiction par modélisation du comportement d’un procédé [3]. 
IL.1.3 Modélisation et logiciels 
Les plans d'expériences permettent d'organiser au mieux les essais qui accompagnent une 
recherche scientifique ou des études industrielles. Ils sont applicables à de nombreuses disciplines et 
à toutes les industries à partir du moment où l'on recherche le lien qui existe entre une grandeur 
d’intérêt y, tel que le rapport d'amélioration de la conductivité thermique ou autre, et des variables 
x; telle que la température, et ceci dans un but d'optimisation [2, 4]. 


Il faut penser aux plans d’expériences si l’on s’intéresse à une fonction du type : 


y = f(x) (1) 


Raison pour laquelle il existe des logiciels pour les traiter comme Mini-Tab ou 
principalement JMP (John's Macintosh Product), pour ne citer que les plus connus [5]. 

Les corrélations sont basées sur une exploitation statistique de résultats expérimentaux 
souvent associées à des considérations physiques tandis que les modèles sont, en premier lieu, basés 
sur des considérations phénoménologiques enrichies par des données expérimentales, ce qui leur 


donne, a priori, un caractère plus universel. 
IL.1.4 Méthodologie des plans d’expériences 
11.1.4.1 Notion d'espace expérimental 


Un plan d’expériences peut être défini comme une suite d’essais organisés à l’avance, de 
manière à déterminer, au moyen d’un nombre réduit d’essais et avec un maximum de précision, 


l'influence de multiples paramètres (facteurs) sur une (ou plusieurs) propriété (s) étudiée(s) 


(réponses). 
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De manière générale, la méthode des plans d’expériences cherche à déterminer et à établir les liens 
existant entre deux types de variables : 

e La réponse : grandeur physique étudiée ; 

e Les facteurs : grandeurs physiques modifiables par l’expérimentateur, sensées influer sur 

les variations de la réponse. 

Plus précisément, la méthode des plans d’expériences vise aussi bien à comprendre les relations 
liant la réponse avec les facteurs, que les facteurs entre eux [6]. 
La grandeur d’intérêt, qui est généralement notée y, porte le nom de réponse. Les variables qui 
peuvent modifier la réponse sont appelées facteurs. On parle donc des facteurs qui influent une 
réponse. Les termes facteur et réponse sont universellement employés dans le domaine des plans 
d'expériences [5]. 

La construction des plans et l’interprétation des résultats dépendent en grande partie des 
types de facteurs utilisés dans l’étude. On distingue plusieurs types de facteurs, nous retiendrons les 
types de facteurs suivants : les facteurs continus, les facteurs discrets, les facteurs ordonnables et les 
facteurs booléens. 

Un expérimentateur qui lance une étude s'intéresse à une grandeur qu'il mesure à chaque 
essai. Cette grandeur s'appelle la réponse, c'est la grandeur d'intérêt. La valeur de cette grandeur 
dépend de plusieurs variables. Au lieu du terme variable on utilisera le mot facteur. Le premier 
facteur peut être représenté par un axe gradué et orienté, figure (II-1). La valeur donnée à un facteur 
pour réaliser un essai est appelée niveau. Lorsqu'on étudie l'influence d'un facteur, en général, on 
limite ses variations entre deux bornes. La borne inférieure qui est le niveau bas et la borne 
supérieure qui est le niveau haut [5, 6]. 

Domaine du facteur 


Facteur 


AN 


Niveau haut 


Figure IL.1: Domaine du facteur [5]. 


S'il y a un second facteur, 1l est représenté, lui aussi, par un axe gradué et orienté. On le 
définit, comme pour le premier facteur par son niveau haut, son niveau bas et son domaine de 
variation. Ce second axe est disposé orthogonalement au premier. On obtient ainsi un repère 


cartésien qui définit un espace euclidien à deux dimensions. Cet espace est appelé l'espace 


expérimental (figure IL.2) [5, 6]. 
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Espace expérimental 


Facteur 2 


Figure IL.2: Espace expérimental [5]. 


Le niveau x, du facteur 1 et le niveau x2 du facteur 2 peuvent être considérés comme les 
coordonnées d'un point de l'espace expérimental. Une expérience donnée est alors représentée par 
un point dans ce système d’axes. Par conséquent, un plan d'expériences est représenté par un 
ensemble de points expérimentaux, comme montré dans la figure IL.3 ci-dessous. Cet ensemble de 
points définit ce qu’on appelle le domaine d'étude, qui représente l'espace expérimental choisie par 


l'expérimentateur pour faire ses essais [6]. 


Figure IL.3: Domain d’étude [5]. 


Les définitions qui ont été données s'appliquent bien aux variables continues. Mais il existe 
d'autres types de variables comme les variables discrètes. On peut encore parler d'espace 
expérimental mais il n'aura pas les mêmes propriétés que l'espace des variables continues. Il y a 
également les grandeurs ordonnables comme, par exemple, des distances qui peuvent être courtes, 
moyennes et longues. Là aussi, la notion d'espace expérimental existe toujours mais cet espace 
possède des propriétés différentes des deux premiers [6]. 


11.1.4.2 Notion de surface de réponse 


Les niveaux x; représentent les coordonnées d'un point expérimental et y est la valeur de la 
réponse en ce point. On définit un axe orthogonal à l'espace expérimental et on l'attribue à la 
réponse. La représentation géométrique du plan d'expériences et de la réponse nécessite un espace 
ayant une dimension de plus que l'espace expérimental. Un plan à deux facteurs utilise un espace à 


trois dimensions pour être représenté : une dimension pour la réponse, deux dimensions pour les 
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facteurs. À chaque point du domaine d'étude correspond à une réponse. L’ensemble de tous les 
points du domaine d'étude correspond à un ensemble de réponses qui se localisent sur une surface 
appelée la surface de réponse (figure 11.4). Le nombre et l'emplacement des points est le problème 
fondamental des plans d'expériences. On cherche à obtenir la meilleure précision possible sur la 


surface de réponse tout en limitant le nombre d’expériences [6]. 


Facteur! 


Figure IL.4: Illustration de la surface de réponse [4]. 


11.1.4.3 Modélisation mathématique 

On choisit a priori une fonction mathématique qui relie la réponse aux facteurs en prenant un 
développement limité de la série de Taylor-Mac Laurin. Les dérivées sont supposées être constantes 
et le développement prend la forme d'un polynôme de degré plus ou moins élevé : 
Y= GE a FLE test TE aiX;? + if. N EXT 0 Xy (IL2) 
Où : 

e  yest la réponse ou la grandeur d'intérêt. Elle est mesurée au cours de l'expérimentation et 
elle est obtenue avec une précision donnée. 

e x; représente le niveau attribué au facteur à par l'expérimentateur pour réaliser un essai. 
Cette valeur est parfaitement connue. On suppose même que ce niveau est déterminé sans 
erreur (hypothèse classique de la régression). 

° do, Gi, &j, 4j; Sont les coefficients du modèle mathématique adopté à priori. Ils ne sont pas 


connus et doivent être calculés à partir des résultats des expériences. 


L'intérêt de modéliser la réponse par un polynôme est de pouvoir calculer ensuite toutes les 
réponses du domaine d'étude sans être obligé de faire les expériences. Ce modèle est appelé 


"modèle postulé" ou "modèle a priori" [7]. 


Deux appoints doivent être apportés au modèle précédemment décrit. Le premier 
complément est le "manque d'ajustement". Cette expression traduit le fait que le modèle à priori est 


fort probablement différent du modèle réel qui régit le phénomène étudié. Il y a un écart entre ces 


deux modèles qui est le manque d'ajustement (Lack of fit en Anglais). 
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Le second complément est la prise en compte de la nature aléatoire de la réponse. En effet, si 
l'on mesure plusieurs fois une réponse en un même point expérimental, on n'obtient pas exactement 
le même résultat. Les résultats sont dispersés, les dispersions ainsi constatées sont appelées erreurs 
expérimentales. Ces deux écarts, c.à.d. manque d'ajustement et erreur expérimentale, sont souvent 


réunis dans un seul écart (notée e) [6]. 
Le modèle utilisé par l'expérimentateur s'écrit alors : 


Y= dt EGr FL Mr tt TE aiX;? + if. XX Xe + (IL3) 


I1.1.4.4 Types de plan d’expériences 

Un plan d’expériences définit un ensemble d’expériences réparties dans le domaine 
expérimental. Il existe de nombreux types de plans d’expériences et le choix de l’un d’eux peut 
s’avérer difficile. Le principe de ce choix est de sélectionner le plan le plus avantageux, présentant 
le moins d’expériences possible pour obtenir le plus d’informations possible. Ce choix dépend de 
l'objectif fixé, des méthodes d’analyse utilisées, des contraintes expérimentales, des ressources, du 
comportement étudié (lié aux facteurs et aux variables d’intérêt) et de l’expertise à disposition. 
Nous nous contenterons ici d’expliquer brièvement les types de plan d’expériences les plus utilisées 
[8]. 

A. Plan factoriel complet à deux niveaux 

Ce sont des plans très adaptés vu qu’ils sont plus rapides ainsi simples à exploiter. Ils sont 
notés 2°, alors qu’ils comportent k facteurs étudiés et le 2 signifie le nombre de niveaux par 


facteurs. Toutes les combinaisons de niveaux sont effectuées au cours de l’expérimentation. 


Ces plans peuvent être utilisés indistinctement pour les variables continues et pour les 
variables discrètes, pour un plan factoriel complet à 3 facteurs noté 2°. Les points expérimentaux 


sont représentés aux sommets d’un cube, comme représenté dans la Figure ILS ; 


Figure IL.5: Disposition des points expérimentaux d’un plan 23[6]. 
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La matrice des expériences représentée sur le Tableau II-1 traduit l'emplacement de ces 


points en indiquant leurs trois coordonnées en variables codées [3]: 


Tableau IL.1: Matrice d’un plan d’expériences 23 et résultats expérimentaux. 


N° de Zi 2 Z3 Réponse 
l'essai x, x x, y 

1 1 1 1 Yi 

2 2] =) 1 _ 

3 -1 1 -1 Ya 

4 1 1 1 Ya 

5 1 -1 1 Ys 

6 1 2 1 Ye 

7 1 1 -1 Y7 

8 1 1 1 Ye 


B. Plans composites centrés 
Un plan composite est constitué de trois parties : 
— Un plan factoriel dont les facteurs prennent deux niveaux. 
— Au moins un point expérimental situé au centre du domaine d’étude. 
— Des points axiaux en étoile. Ces points expérimentaux sont situés sur les axes de chacun des 
facteurs. 
Ces plans sont souvent utilisés lorsque le plan demande une expérimentation séquentielle, car ils 
peuvent intégrer des informations provenant d'une expérience factorielle correctement planifiée [2]. 
C. Plans factoriel à niveaux mixtes 
Il y a autant de facteurs que l'on veut et chaque facteur prend le nombre de niveaux 
nécessaires à la bonne exécution de l'étude. Par exemple, un facteur peut prendre 3 niveaux, un 


autre 4 niveau et un troisième 6. Ces plans sont très utilisés pour les variables discrètes [7]. 
11.1.4.5 Analyses statistique de l’équation de régression 


a. Effets significatifs (test de Student) : 
Cas où chaque essai est répété m fois ; on détermine la moyenne arithmétique des résultats 


des essais parallèles [6] : 


ÿ; = Ad j = 1,2,..,N (1.4) 


m 
Puis la variance de sondage : 


2 — Zu=1iu-Ÿi) A nn (LS) 


t m-1 


Où m est le nombre de répétition pour chaque essai et N est le nombre d’expériences. Si les 


variances de sondage sont homogènes, on calcule alors la variance de reproductibilité : 
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2 LS 
Sep = “+ (IL.6) 


b. Validation de l’équation de régression 
Avant d’accepter le modèle postulé, il est impératif de vérifier que ce dernier décrit 
correctement les variations de la réponse en fonction des facteurs et cela en effectuant les tests 


suivants : 


b.1. Recherche de biais du modèle : Afin de chercher la nature des erreurs qui peuvent être 
engendrées par les variables non contrôlées, on a recours à la comparaison de la variance 
résiduelle (caractérisant la différence entre les valeurs mesurées et celle calculées) avec la variance 
de reproductibilité (qui caractérise la dispersion maximale) [9, 10]. L’égalité ou non des deux 
variances résiduelles et reproductible se basent sur le test de Fisher qui s’énonce comme suit : 


ç2, = LmOi-P)? 


SiF = % € F(0.95,N— I Di u IL7 
LE = ep < ( 79, N — LI, No — ) D (Hi-Yo)2 ( ° ) 
SŸep — foi 


b.2. Test de signification du modèle : Lorsque le modèle est sans biais, il est possible 
d’effectuer le test de signification de la régression. En effet le reste de Student permet d’alléger le 
modèle mathématique en éliminant les facteurs et les interactions qui n’ont pas de conséquence 
significative sur la réponse [11]. Nous pouvons juger si la part de la dispersion de la réponse 
expliquée par le modèle est grande devant la dispersion résiduelle, en comparant les variances 
correspondantes par le test de Fisher : 


_ 2LO-ÿ)2/1-1 
EL Gi-m)2/N-1 (IL.8) 


y représente la moyenne des essais donnés par : 


in — re Vi 
ÿ = un (IL.9) 


b.3. Coefficient de détermination : Le coefficient de détermination appelé aussi coefficient 


de corrélation R° est défini par : 


_ DC 7E2E 


2 
— A Gi-p)2 (IL.10) 


Plus le coefficient de détermination est proche de 1, plus le modèle est représentatif c’est-à-dire, 


meilleur. 
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b.4. Analyse des résidus : Lorsqu’on utilise le modèle pour estimer la réponse, on mesure 
généralement un écart appelé résidu définit comme la différence entre la valeur mesurée et la valeur 


prédite par le modèle, soit : 


e = ÿ, — y; (L.11) 


Si la moyenne des erreurs entre la surconsommation donnée par le modèle et celle mesurée n’atteint 


pas 5% alors l’équation du modèle simule parfaitement l’expérience. 


11.1.4.6 Fonction désirabilité 

La fonction désirabilité est très utile quand il faut trouver le meilleur compromis entre 
plusieurs réponses. Cette fonction a été proposée par Derringer et Suich et elle figure dans de 
nombreux logiciels de plans d’expériences [5]. Elle est très facile à utiliser et rend de grands 
services. Cette fonction varie de 0 à 1. La valeur 1 correspond à une satisfaction maximale et la 


valeur 0 à un refus complet. 


Décollement 


Peu de 
décollement 


(La réponse est 
satisfaisante) 


Décollement 

de plus en plus 
acceptable 

(La réponse est 
plus ou moins 
satisfaisante) 


Trop de 
décollement 


(La réponse ne 
convient pas) 


d=1 Désirabilité 


Figure IL.6: La fonction désirabilité varie de 0 à 1 quand la réponse s’améliore progressivement 


[5]. 
IL.2 Les Réseaux de Neurones Artificiels 


IL.2.1 Introduction 

Les techniques de l’intelligence artificielle sont connues actuellement pour leur grande 
potentialité de pouvoir résoudre les problèmes liés aux processus industriels, notamment le contrôle, 
l'estimation et l’identification des paramètres des systèmes variant. Parmi ces techniques, on trouve 
les Réseaux de Neurones Artificiels (RNA), qui est une technique de traitement de données ayant 
des fonctionnalités d’identification, de contrôle ou de filtrage plus efficaces [12]. 

Depuis une vingtaine d'années, l'utilisation des réseaux de neurones artificiels (RNA) s'est 


étendue dans des domaines très divers de l'industrie et des services. Un Réseau de Neurones 
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Artificiel est un modèle de calcul dont la conception est très schématiquement inspirée du 
fonctionnement des vrais neurones L’origine des réseaux de neurones vient de l’essai de 
modélisation du neurone biologique par Warren Mc Culloch et Walter Pitts [13]. Ils supposent que 
l'impulsion nerveuse est le résultat d’un calcul simple effectué par chaque neurone et que la pensée 
née grâce à l’effet collectif du réseau de neurones interconnectés. Et comme définition on peut dire 
que : « Les réseaux de neurones artificiels sont des réseaux fortement connectés de processeurs 
élémentaires fonctionnant en parallèle. Chaque processeur élémentaire calcule une sortie unique 
sur la base des informations qu'il reçoit et toute structure hiérarchique de réseaux est évidemment 
un réseau » [14]. 

Dans cette partie, nous présentons dans un premier temps des généralités et les concepts de 
base des RNA. Nous aborderons ensuite le processus d’apprentissage d’un RNA qui se réalise 
généralement par le choix du type d’apprentissage et de l’algorithme de mise à jour. Comme type, 


seul l’apprentissage à rétro-propagation du gradient d’erreur sera détaillé. 


IL.2.2 Domaines d'application des réseaux de neurones artificiels 
Aujourd’hui, les réseaux de neurones artificiels ont de nombreuses applications dans des secteurs 
très variés [13] : 
1. Traitement d’images : reconnaissance de caractères et de signatures, compression d’images, 
reconnaissance de forme, cryptage, classification, etc. 
2. Traitement du signal : filtrage, classification, identification de source, traitement de la 
parole.….…etc. 
3. Contrôle : commande de processus, diagnostic, contrôle qualité, asservissement des robots, 
systèmes de guidage automatique des automobiles et des avions..….etc. 
4. Défense : guidage des missiles, suivi de cible, reconnaissance du visage, radar, sonar, lidar, 
compression de données, suppression du bruit...etc. 
Optimisation : planification, allocation de ressource, gestion et finances, etc. 
Simulation : simulation du vol, simulation de boîte noire, prévision météorologique, recopie 


de modèle. .….etc. 
II.2.3 Fonctionnement des réseaux de neurones 


L'origine des réseaux de neurones vient de l’essai de modélisation mathématique du cerveau 
humain, en supposant que l’impulsion nerveuse est le résultat d’un calcul simple effectué par 
chaque neurone et que la pensée née grâce à l’effet collectif d’un réseau de neurone interconnecté 


(Figure IL.7) et (Tableau IL.2) [15]. 


Approche mathématique 


(a) 


Neurones Biologiques 


nl 
ù 
à 
À 
ke À 
S = 
S'e 
Ÿ à 
& à 
Ÿ 
ÿ € 
à 
& ÿ 
à 
>= 
Ÿ 
L2 


inconscient 


Figure IL.7: Structure d’un réseau de neurone biologique et le cerveau humain [12]. 


Tableau IL.2: Analogie entre le neurone biologique et le neurone formel [15]. 


Neurone artificiel Neurone biologique 
Poids de connexion Synapses 
Signal de sortie Axones 
Signal d'entrée Dendrite 
Fonction d'activation Soma 


Un réseau de neurones est un assemblage de constituants élémentaires interconnectés 
(appelés "neurones" en hommage à leur modèle biologique), qui réalisent chacun un traitement 
simple mais dont l’ensemble en interaction fait émerger des propriétés globales complexes. Chaque 
neurone fonctionne indépendamment des autres, de telle sorte que l'ensemble forme un système 
massivement parallèle. L'information est stockée de manière distribuée dans le réseau sous forme de 
coefficients synaptiques ou de fonctions d'activation, 1l n'y a donc pas de zone de mémoire et de 
zone de calcul, l'une et l'autre sont intimement liés. 

Un réseau de neurone ne se programme pas, il est entraîné grâce à un mécanisme 
d'apprentissage. Les tâches particulièrement adaptées au traitement par réseau de neurones sont : 
l’association, la classification, la discrimination, la prévision ou l'estimation et la commande de 
processus complexes. Les réseaux de neurones artificiels consistent en des modèles plus ou moins 
inspirés du fonctionnement cérébral de l’être humain, en se basant principalement sur le concept de 


neurone [12]. 
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IL.2.4 Le modèle neurophysiologique [14] 

Le cerveau humain contient environ 100 milliards de neurones. Ces neurones permettent, 
entre autre, de lire par exemple un texte tout en maintenant une respiration régulière permettant 
d’oxygéner le sang et actionnant le cœur qui assure une circulation efficace de ce sang nourrissant 


ainsi les cellules, etc. Chacun de ces neurones est par ailleurs fort complexe. 


Figure IL.8: Modèle d'un neurone biologique. 


Essentiellement, 1l s’agit de tissu vivant et de chimie. Les spécialistes des neurones 
biologiques (neurophysiologistes) commencent à peine à comprendre quelques-uns de leurs 
mécanismes internes. On croit en général que leurs différentes fonctions neuronales, y compris 
celle de la mémoire, sont stockées au niveau des connexions (synapses) entre les neurones. C’est ce 
genre de théorie qui a inspiré la plupart des architectures de réseaux de neurones artificiels. 
L’apprentissage consiste alors soit à établir de nouvelles connexions, soit à en modifier des 
existantes. 

Un neurone est une cellule particulière comme la montre la figure IL.8. Elle possède des extensions 
par lesquelles elle peut distribuer des signaux (axones) ou en recevoir (dendrites). 

Dans le cerveau, les neurones sont reliés entre eux par l’intermédiaire des axones et des 
dendrites. En première approche. On peut considérer que ces sortes de filaments sont conductrices 
d'électricité et peuvent ainsi véhiculer des messages depuis un neurone vers un autre. Les dendrites 
représentent les entrées du neurone et son axone sa sortie. 

Un neurone émet un signal en fonction des signaux qui lui proviennent des autres neurones. On 
observe en fait au niveau d'un neurone, une intégration des signaux reçus au cours du temps, c'est à 
dire une sorte de sommations des signaux. En général, quand la somme dépasse un certain seuil, le 


neurone émet à son tour un signal électrique. 


IL.2.5 Modélisation d’un neurone formel [14, 15] 
Les réseaux de neurones artificiels sont des modèles qui peuvent être décrits par leurs 


composants, leurs variables descriptives et les interactions des composants. 
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Ÿ” Structure 
Un "neurone formel" est une fonction algébrique non linéaire et bornée, dont la valeur dépend 
des paramètres appelés coefficients ou poids. Les variables de cette fonction sont habituellement 
appelées "entrées" du neurone et la valeur de la fonction est appelée sa "sortie". 
Un neurone est donc avant tout un opérateur mathématique, dont on peut calculer la valeur 
numérique par quelques lignes de logiciel. La représentation graphique d’un neurone est comme 


indiqué dans la figure IT.9. 


Noyau du neurone 


Figure IL.9: Modèle d'un neurone artificiel [15]. 


À partir des observations de neurone biologique, découle le modèle du neurone formel 
(proposé par W. M. Culloch et W. Pitts en 1943) où : 
+ x; : Représentent les vecteurs d’entrées, elles proviennent soit des sorties d’autres neurones, soit 
de stimuli sensoriels (capteur visuel, sonore.….). 
+ _W;;: Sont les poids synaptiques du neurone j, ils correspondent à l’efficacité synaptique dans les 
neurones biologiques (w;; > 0 pour synapse excitatrice; w;; < 0 pour synapse inhibitrice). Ces 
poids pondèrent les entrées et peuvent être modifiés par apprentissage. 
+ Biais : L’entrée prend souvent les valeurs -1 ou +1 qui permet d’ajouter de la flexibilité au réseau 
en permettant de varier le seuil de déclenchement du neurone par l’ajustement des poids et du biais 
lors de l’apprentissage. 
+ _ Noyau : Il intègre toutes les entrées ainsi le biais ; et il calcul la sortie du neurone selon une 
fonction d’activation (souvent non linéaire) pour donner une plus grande flexibilité d’apprentissage. 
Le modèle mathématique d’un neurone artificiel est illustré dans la figure IL.9. Un neurone est 
essentiellement constitué d’un intégrateur qui effectue la somme pondérée de ses entrées. Le 
résultat de cette somme est ensuite transformé par une fonction de transfert f qui produit la sortie y 
du neurone. 

En suivant les notations présentées dans la section précédente, les n entrées du neurone 


correspondent au vecteur x = [x1,xX2,x3..x, |! alors que w = [w11,W21,W31..wn1]T représente le 


vecteur des poids du neurone. 
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La sortie s de l’intégrateur est donnée par l’équation suivante : 

S = Da Wiaxi + D = Winx + WoiX2 + WaaXs + + + WniXn (IL.12) 
Que l’on peut aussi écrire sous forme matricielle : 

s=wTx+b (IL. 13) 


Cette sortie correspond à une somme pondérée des poids et des entrées plus ce qu’on nomme le 
biais b du neurone. Le résultat s de la somme pondérée s’appelle le niveau d’activation du 
neurone. Le biais b s’appelle aussi le seuil d’activation du neurone. Lorsque le niveau d’activation 
atteint ou dépasse le seuil b, alors l’argument de f devient positif (ou nul) ; sinon il est négatif. 

Ÿ”_ Fonctions d'activations 

Différentes fonctions de transfert peuvent être utilisées comme fonction d’activation dans les 

RNA (tableau IL3). Les fonctions d'activations les plus utilisées sont les fonctions "seuil" (en 
anglais "hard limit"), qui sont linéaire et sigmoïde. Comme son nom l’indique, la fonction seuil 
applique un seuil sur son entrée. Plus précisément, une entrée négative ne dépasse pas le seuil et la 
fonction retourne alors la valeur 0 (ce 0 peut être interprété comme significativement faux) ; alors 
qu’une entrée positive ou nulle ne dépasse pas le seuil et la fonction retourne la valeur 1 (vrai). 

Utilisée dans le contexte d’un neurone, cette fonction est illustrée à la figure IL.10.a. On 
remarque alors que le biais b dans l'expression de y = hard lim(wTx — b) détermine 
l'emplacement du seuil sur l’axe wTx, où la fonction passe de 0 à 1 permettant ainsi de prendre des 


décisions binaires. 


Figure IL.10: Fonctions d’activations : (a) du neurone « seuil » ; (b) du neurone « linéaire », et (c) 


du neurone « sigmoïde ». 


La fonction linéaire est très simple, elle affecte directement son entrée à sa sortie : 
y = S. Appliquée dans le contexte d’un neurone, cette fonction est illustrée à la figure IL.10.b. Dans 
ce cas, la sortie du neurone correspond à son niveau d’activation dont le passage à zéro se produit 
lorsque wTx = b. 


La fonction de transfert sigmoïde est quant à elle illustrée à la figure IL.10.c. Son équation 


est donnée par : 
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ÿ = (IL.14) 


E 1+exp”$ 
Elle ressemble soit à la fonction seuil soit à la fonction linéaire selon que l’on est loin ou 
près de b, respectivement. La fonction seuil est amplement non linéaire car il y a une discontinuité 
lorsquewTx = b. De son côté, la fonction linéaire est tout à fait linéaire et elle ne comporte aucun 
changement de pente. La sigmoïde est un compromis intéressant entre les deux. Notons finalement, 


que la fonction "tangente hyperbolique (tanh)" est une version symétrique de la sigmoïde. 


Tableau II.3: Différentes fonctions d'activations utilisées dans les RNA. 


Nom de La fonction Relation entrée/sortie | 1cône | Nom MATALB 
_ y =0 si s<0 
EEE 
y=1 si s>0 


Linéaire 


Linéaire saturée 


Linéaire saturée symétrique 3 


s s 
s s< 

s si £<s= 

y=1 si s> 

y =0 si s<0 

Linéaire positive poslin 
py=s si s2z20 
1 
Sigmoïde y 
” l+exp 


Tangente hyperbolique 


= s s maximum 
Compétitive 
autrement 


I1.2.6 Architecture des réseaux de neurones 


L'architecture est un concept très important qui joue un rôle déterminant dans la 
classification des RNA. Dans la littérature on utilise souvent le mot structure comme architecture. 
Chaque architecture a sa propre organisation (il s’agit de la façon dont ils sont ordonnés et 
connectés) qui est adapté à des applications bien spécifiques. 

La majorité des réseaux de neurones utilise le même type de neurones. Quelques 
architectures plus rares se basent sur des neurones dédiés. L’architecture d’un réseau de neurones 
dépend de la tâche à apprendre. Un réseau de neurone est en général composé de plusieurs couches 
de neurones, des entrées jusqu’aux sorties. On distingue deux grands types d'architectures de 
réseaux de neurones : les réseaux de neurones non bouclés et les réseaux de neurones bouclés [12, 
15]. 

11.2.6.1 Les réseaux de neurones non bouclés [16] 

Un réseau de neurones non bouclé réalise une (ou plusieurs) fonctions algébriques de ses 


entrées, par composition des fonctions réalisées pour chacun de ses neurones. Un réseau de 
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neurones non bouclé est représenté graphiquement par un ensemble de neurones "connectés" entre 
eux, avec l’information circulant des entrées vers les sorties sans "retour en arrière". Si l’on 
représente le réseau comme un graphe dont les nœuds sont des neurones et les arêtes des 
"connexions" entre ceux-ci, le graphe d’un réseau non bouclé est acyclique. Le terme de 
"connexions" est une métaphore, car dans la très grande majorité des applications, les réseaux de 
neurones sont des formules algébriques dont les valeurs numériques sont calculées par des 
programmes d’ordinateurs et non par des objets physiques (circuits électroniques spécialisés). 
Néanmoins, le terme de connexion, issu des origines biologiques des réseaux de neurones, est passé 
dans l’usage, car il est commode quoique trompeur. Il a même donné naissance au terme de 
connexionnisme. 
>» Réseaux de neurones monocouches 

La structure d’un réseau monocouche est telle que des neurones organisés en entrée sont 

entièrement connectés à d’autres neurones organisés en sortie par une couche modifiable de poids 


(figure IL.11). 


Couche Couche 
d’entrée de sortie 


Figure IL.11: Schéma d'un réseau de neurones monocouche. 


» Réseaux de neurones multicouches 

Les neurones sont arrangés par couche. Il n'y a pas de connexion entre neurones d'une même 
couche et les connexions ne se font qu'avec les neurones de couches avales. Habituellement, chaque 
neurone d'une couche est connecté à tous les neurones de la couche suivante et ainsi de suite. Ceci 
nous permet d'introduire la notion de sens de parcours de l'information (de l'activation) au sein d'un 
réseau et donc définir les concepts de neurone d'entrée et de neurone de sortie. Par extension, on 
appelle couche d'entrée l'ensemble des neurones d'entrée et couche de sortie l'ensemble des 
neurones de sortie. Les couches intermédiaires n'ayant aucun contact avec l'extérieur sont appelées 
couches cachées. 

La figure IT. 12 représente un réseau de neurones non bouclé qui a une structure particulière, 


très fréquemment utilisée. Il comprend des entrées, deux couches de neurones cachés et des 
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neurones de sortie. Les neurones de la couche cachée ne sont pas connectés entre eux. Cette 


structure est appelée Perceptron multicouches [15]. 


Couche Couche ‘oucl Couche de 
d'entrée cachée 1 ée sortie 


Figure IL.12: Schéma d'un réseau de neurones non bouclé (Perceptron multicouches). 


On note aussi que les réseaux multicouches sont beaucoup plus puissants que les réseaux 
simples à une seule couche. En utilisant deux couches (une couche cachée et une couche de sortie), 
à condition d’employer une fonction d’activation sigmoïde sur la couche cachée, on peut entraîner 
un réseau à produire une approximation de la plupart des fonctions avec une précision arbitraire 
(cela peut cependant requérir un grand nombre de neurones sur la couche cachée). Sauf dans des 
rares cas, les réseaux de neurones artificiels exploitent deux ou trois couches. 

Les réseaux de neurones non bouclés sont des objets statiques si les entrées sont 
indépendantes du temps et si les sorties le sont également. Ils sont utilisés principalement pour 
effectuer des tâches d'approximation de fonction non linéaire, de classification ou de modélisation 
de processus statiques non linéaire [15, 16]. 

11.2.6.2 Les réseaux de neurones bouclés 

Contrairement aux réseaux de neurones non bouclés dont le graphe de connexions est 
acyclique, les réseaux de neurones bouclés peuvent avoir une topologie de connexions quelconque, 
comprenant notamment des boucles qui ramènent aux entrées la valeur d'une ou plusieurs sorties. 
Pour qu'un tel système soit causal, 1l faut évidemment qu'à toute boucle soit associé un retard. Un 
réseau de neurones bouclé est donc un système dynamique, régi par des équations différentielles ; et 
comme l'immense majorité des applications sont réalisées par des programmes d'ordinateurs, on se 
place dans le cadre des systèmes à temps discret où les équations différentielles sont remplacées par 
des équations aux différences finies. 

Il s’agit donc de réseaux de neurones avec retour en arrière (feedback network ou recurrent 


network), (Figure IL.13) [15, 16]. 
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Figure IL.13: Schéma de réseau de neurones bouclé. 


IL.2.7 L'apprentissage 

Parmi les propriétés des réseaux de neurones les plus fondamentales c’est sûrement la 
capacité d’apprendre de son environnement et d’améliorer sa performance à travers un processus 
d'apprentissage. 

L'apprentissage est une étape très importante du développement d'un réseau de neurones 
durant laquelle le comportement du réseau est modifié itérativement jusqu'à l'obtention du 
comportement désiré, et ce par l'ajustement des poids des neurones à une source d'information bien 
définit. L'apprentissage consiste également à extraire des conformités parmi les données utilisées 
pour l'entrainement du réseau, mais l'objectif essentiel de l'apprentissage est la résolution du 
problème par la prévision, la classification, etc [15]. 

Et comme définition on peut dire que : « L'apprentissage est une phase du développement d'un 
réseau de neurones durant laquelle le comportement du réseau est modifié jusqu'à l'obtention du 
comportement désiré. L'apprentissage neuronal fait appel à des exemples de comportement », [14]. 

Pour un RNA, l'apprentissage peut être considéré également comme étant le processus de la 
mise à jour des poids (connexion) au sein du réseau et ceci dans le but d'ajuster la réponse du réseau 
à l'expérience et aux exemples. Il existe plusieurs types de règles d'apprentissage qui peuvent être 
rassemblées en deux catégories qui sont : l'apprentissage supervisé et l'apprentissage non supervisé. 

Des exemples d'utilisations de différentes règles d'apprentissage, d’architectures et 
d’algorithme qui sont relié à ces types supervisé et non supervisé sont illustrés dans les tableaux 
IL.1 et IIL2 dans l’annexe B. 

11.2.7.1 Stratégies d'apprentissage 

L’apprentissage au sein des différentes architectures dépend de l’architecture du réseau et de 
l’environnement du problème. Les deux règles d’apprentissage pour corriger les poids d’un neurone 
(règle de Hebb et de Widrow) ne concernent qu’un seul neurone. Ces règles peuvent servir pour 
mettre à jour les poids d’un neurone, de certains réseaux de neurones, mais ne peuvent être 


généralisées et s’appliquer à n’importe quelle architecture. Chaque architecture possède ses 


spécificités et nécessite une règle d’adaptation des poids qui lui est propre. 
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L'apprentissage n’est pas modélisable dans le cadre de la logique déductive, celle-ci en effet 
procède à partir de connaissances déjà établies dont on tire des connaissances dérivées. Or il s’agit 
ici de la démarche inverse, par observations limitées tiré des généralisations plausibles. 

La notion d’apprentissage recouvre deux réalités (figure IL.14) [14] : 

* La mémorisation : le fait d’assimiler sous une forme dense des exemples éventuellement 
nombreux. 

+ La généralisation : le fait d’être capable, grâce aux exemples appris, de traiter des exemples 


distincts, encore non rencontrés, mais similaires. 


Généralisation 


Mémorisation 


Figure IL.14: Principe de l'apprentissage. 


11.2.7.2 Algorithmes d'apprentissages 

Définir l'architecture d'un RNA pour un problème quelconque revient à faire un choix 
judicieux de la taille du réseau, du nombre total de couches et de neurones, de distribution des 
données et des fonctions de transfert. Le choix de ces paramètres dépend de l'utilisateur. 

Il n'existe pas dans la littérature pour le moment des données suffisantes qui peuvent 
déterminer clairement les paramètres à adopter pour résoudre un problème donné. Les différentes 
techniques présentées par différents auteurs sont destinées à résoudre des cas particuliers. 

Un ensemble de règles bien définies permettant de réaliser un tel processus d’adaptation des 
poids constitue ce qu’on appelle l’algorithme d’apprentissage du réseau. 

Il y a plusieurs types d'algorithmes d’apprentissage tels que : l’apprentissage par correction 
d’erreur, l’apprentissage à rétro-propagation du gradient d’erreur et l’apprentissage par 
renforcement [15]. 


Dans la suite de ce chapitre, nous nous sommes intéressés à l’algorithme de la rétro- 


propagation du gradient d'erreur (back-propagation). 


IL.2.8 Le Perceptron Multicouches (MLP) 
11.2.8.1 Structure du réseau MLP 
Une seule couche de neurones ne pouvant réaliser que des séparations linéaires, l'idée vient 


alors de rajouter des couches dites cachées pour réaliser un réseau de neurone multicouche. Dans 


une couche, les neurones ne sont pas connectés entre eux. 
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Un réseau à couches est une extension du célèbre perceptron avec une ou plusieurs couches 
intermédiaires appelées "cachées". Le perceptron multicouches (Multi Layered Perceptron - MLP) 
sont les réseaux de neurones les plus connus. Un perceptron est un réseau de neurones artificiel du 
type “feedforward”", c'est-à-dire à propagation directe [12]. Le schéma de la figure IL.15 montre un 
réseau à trois couches possédant trois entrées et une sortie. La première est celle des entrées (elle 
n’est cependant pas considérée comme couche neuronale par certains auteurs car elle est linéaire et 
ne fait que distribuer les variables d’entrées). La deuxième est dite couche cachée et constitue le 
cœur du réseau de neurones. La troisième, constituée ici par un seul neurone est la couche de sortie, 
on numérote les couches de l’entrée vers la sortie donc sur le schéma de gauche à droite, la couche 
cachée a trois neurones, un tel réseau est ainsi noté par 3-3-1. 

On remarque sur la figure IL.14, des termes x$° en entrée des neurones (le terme en exposant 
représente, non pas la fonction puissance, mais plutôt l’indice (”) de la couche du réseau de 


neurones). 


Couche de 


Couche ÿ 
: : sortie 
d'entrée 


Figure IL.15: Réseau de neurones de type perceptron à une couche cachée. 


Le perceptron multicouche (MLP) est très utilisé en identification et en contrôle. Avec une 
couche cachée, il constitue un "Approximateur universel", qui peut approximer n’importe quelle 
fonction entrées/sorties (sous réserve de mettre suffisamment de neurones dans la couche cachée et 
d'utiliser des sigmoïdes pour les fonctions) d’activation. Bien entendu, les théorèmes 
mathématiques ne démontrent pas qu’un réseau à une seule couche cachée est optimal. 
Malheureusement, il n’existe pas de règle générale qui donne le nombre de neurones à retenir pour 
la couche cachée. Le résultat le plus important est certainement le théorème de Hecht-Nielsen, où 
toute fonction continue peut-être implémentée exactement comme un réseau de neurones à trois 
couches ayant n cellules en entrée, 2n+1 cellule en couche cachée et m cellules de sortie. Il faut 
bien noter que ce théorème ne donne aucune indication quant au nombre de connexion (le réseau 
n’est pas toujours totalement connecté), et ne garantit pas que ce nombre de neurones est optimal 


mais suffisant. Enfin, le réseau de neurones fait partie de Réseaux Adaptatifs Non-linéaires, et cela 
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signifie que ses agents (neurones) s’organisent et modifient leurs liens mutuels lors d’une procédure 


fondamentale qu’est l’apprentissage. Pour une tache précise, l’apprentissage du réseau de neurones 
consiste donc à adapter les différents poids [12, 14, 15]. 

11.2.8.2 L’algorithme de la rétro-propagation du gradient d'erreur (back-propagation) 

L’un des algorithmes les plus répandus est celui de la "rétro-propagation" ou "back- 
propagation". Cet algorithme change les poids d’un réseau dont l’architecture est fixée par 
l’opérateur, à chaque fois qu’un exemple y; = f (x;) est présenté. Ce changement est fait de sorte à 
minimiser l’erreur entre la sortie désirée et la réponse du réseau à une entrée x;. Ceci est réalisé 
grâce à la méthode de descente de gradient. A chaque itération le signal d’entrée se propage dans le 
réseau dans le sens entrée-sortie, une sortie est ainsi obtenue, l’erreur entre cette sortie et la sortie 
désirée est calculée puis par rétro-propagation "error back-propagation", des erreurs intermédiaires, 
correspondant à la couche cachée sont ainsi calculées et permettent l’ajustement des poids w;; (t) 
de la couche cachée [151]. 
L’algorithme de rétro-propagation du gradient comporte donc deux phases : 

1. Propagation : à chaque étape, on présente au réseau un exemple en entrée. Cette entrée 
est propagée jusqu’à la couche de sortie. 

2. Correction : À coup sûr, le réseau ne fournira pas exactement ce que l’on attendait. On 
calcule donc une erreur (en général la somme quadratique moyenne des erreurs pour tous 
les neurones de sortie) que l’on rétro-propage dans le réseau. Ce processus est 
interrompu dès que l’erreur globale est estimée suffisante. 


La figure IL.16 montre le processus de rétro-propagation. 
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Figure IL.16: Apprentissage des réseaux de neurone par l’algorithme de rétro-propagation. 
IL.2.9 Conclusion 


Les plans d'expériences permettent d'organiser au mieux les essais qui accompagnent une 
recherche scientifique ou des études industrielles. Ils sont applicables à de nombreuses disciplines et 
à toutes les industries à partir du moment où l’on cherche le lien qui existe entre une grandeur 


d'intérêt «y» et des variables « x; ». L'utilisation des plans d’expériences permet de limiter le 
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nombre d'essais et admet de quantifier l'impact des paramètres expérimentaux sur le phénomène 


étudié. Pour cela, il faut suivre des règles mathématiques et adopter une démarche rigoureuse. 


Les réseaux de neurones formels, tels que nous les avons définis, possèdent une propriété 
remarquable qui est à l'origine de leur intérêt pratique dans des domaines très divers. Le choix 
d'utiliser tel ou tel type de réseau de neurones dépend de l'application mais aussi des capacités de 
traitement du processeur sur lequel ils s'exécutent. L'apprentissage est vraisemblablement la 
propriété la plus intéressante des réseaux neuronaux, cependant il existe plusieurs algorithmes 
utilisés pour faire l’apprentissage des réseaux multicouches, et en général, l’algorithme 
d’apprentissage de rétro-propagation (RP) est très utilisé. 

De même, pour la conclusion de ce chapitre, on rappellera ce qui suit : 

L'application de la méthode des plans d’expériences ainsi que la méthode des réseaux des 
neurones dans le domaine de transfert de chaleur des nanofluides, permet de générer des modèles 


mathématiques décrivant l’influence des différents facteurs sur la conductivité thermique. 
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Chapitre I / 


Modélisation de la conductivité thermique effective des nanofluides par l’utilisation des 
plans d'expériences à Plan factoriel complet 


IIL.1 Introduction 


Afin d'améliorer le taux de transfert de chaleur des fluides, de nombreux chercheurs 

ont tenté d’augmenter leur conductivité thermique. L’amélioration des propriétés thermiques 
des fluides caloporteurs se fait en y insérant une phase solide de conductivité thermique très 
élevée ; ainsi les nanofluides obtenus présentent une conductivité thermique significativement 
plus élevée que les fluides de base [1-S]. Cette amélioration augmente avec l’augmentation de 
la concentration de nanoparticules. Toutefois, dans leur développement et leur exploitation, il 
convient de trouver un compromis entre des propriétés thermiques améliorées tout en 
maîtrisant leur viscosité (qui ne doit pas être trop importante dans les échangeurs par 
exemple) et leur stabilité (qui doit être assurée), sans quoi les propriétés du fluide peuvent être 
modifiées. Plusieurs paramètres peuvent jouer un rôle dans l’efficacité des nanofluides 
comme la concentration et la nature des nanoparticules, leur densité, leur type, leur taille, leur 
aspect, la température, le fluide de base, la présence d’un surfactant, etc. 
Les propriétés thermophysiques, à savoir la conductivité thermique, la viscosité dynamique, la 
capacité thermique et la masse volumique, des fluides sont fortement modifiées par l’ajout de 
nanoparticules. Cette modification dépend du type de matériau des fluides de base et des 
particules, leur taille, leur concentration, la température du milieu, etc. 

Les propriétés thermophysiques des nanofluides peuvent être contrôlées par les 
spécifications chimiques et physiques du fluide de base et des nanoparticules. 

La conductivité thermique des nanofluides est l'un des facteurs qui jouent un rôle 
important sur les performances de transfert de chaleur. Ainsi, l'utilisation d'un type approprié 
de nanofluides est très importante pour améliorer la vitesse de transfert de chaleur [6, 7]. 
Cependant, les résultats des différentes recherches illustrent de larges écarts et incohérences 
dans les valeurs de la conductivité thermique des nanofluides en raison de différents 
paramètres, tels que la différence dans les techniques de production et de stabilisation des 
nanofluides, la caractérisation des nanofluides ainsi que les méthodes de mesure de la 
conductivité thermique [8]. 

Il existe également une forte corrélation entre la conductivité thermique des 
nanofluides et divers paramètres, à savoir : la fraction volumique des nanoparticules, le 
matériau des particules, leur taille et leur forme, le type et la température du fluide de base, la 
capacité thermique spécifique, le type de surfactant et la densité [9-17]. 

L'objectif général de l'utilisation des nanofluides est d'avoir des propriétés 


thermophysiques optimales à des fractions volumiques minimales de particules. De plus, il a 
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été rapporté que la conductivité thermique des nanofluides augmente avec l'augmentation de 
la concentration de nanoparticules [7, 18-20]. 

Il a également été révélé que le taux d'amélioration de la conductivité thermique 
diminue avec celui du fluide de base, alors qu'il est directement proportionnel à la 
conductivité thermique des nanoparticules [21, 22]. Néanmoins, pour une certaine 
concentration de nanoparticules, 1l est montré que pour une conductivité thermique plus faible 
le coefficient de transfert thermique (sera) est amélioré. Ceci est en accord avec les modèles 
conventionnels comme celui de Maxwell. D'autre part, il a été montré par Alrashed et al [5] 
que l'amélioration de la conductivité thermique des nanofluides est indépendante du fluide de 
base. 

Un autre facteur important affectant la conductivité thermique des nanofluides est la 
taille des nanoparticules, en effet 1l a été confirmé que la conductivité thermique la plus 
élevée est obtenue par une distribution uniforme de petites particules [21]. On suppose 
également qu’une augmentation de la surface spécifique des nanoparticules pourrait conduire 
à une augmentation de la conductivité thermique effective des nanofluides, car le transfert de 
chaleur entre les particules et le fluide de base se produit à l’interface particule/fluide. De 
plus, le mouvement Brownien des nanoparticules pourrait être intensifié par la réduction de la 
taille des particules. On estime donc que ce dernier est l’un des facteurs les plus importants 
pour améliorer le transfert de chaleur des nanofluides [21, 23]. 

La forme des nanoparticules est un autre paramètre affectant le degré d'amélioration de 
la conductivité thermique des nanofluides [24]. Le taux de rehaussement des particules 
allongées est supérieur à celui des particules sphériques [25]. En général, une diminution de la 
sphéricité des nanoparticules conduit à une augmentation de la conductivité thermique des 
nanofluides. 

Généralement, la conductivité thermique des nanofluides est plus sensible à la 
température et l'augmentation de la température entraîne une amélioration de la conductivité 
thermique [4, 20, 26]. De nombreuses études ont montré l'existence d'une température 
optimale à laquelle la conductivité thermique des nanofluides est maximale [7]. 

De plus, le pH a un effet sur la conductivité thermique des nanofluides car son taux 
d'amélioration augmente avec l'acidité [17]. 
Le taux d'amélioration de la conductivité thermique peut aussi être affecté par les 


tensioactifs et les additifs utilisés afin d’empêcher l'agglomération des particules [14, 27]. 
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Un résumé des études expérimentales sur la conductivité thermique des nanofluides peut être 
trouvé dans la référence [28]. Et une revue complète et détaillée des différents modèles peut 
être également trouvée dans les références [29, 30]. 

En raison des écarts importants et incohérences des résultats obtenus par les 
différentes méthodes, des recherches supplémentaires doivent être effectuées afin d’élucider 
les effets de la conductivité thermique des nanofluides. 

En conséquence, le présent travail de recherche vise à trouver une nouvelle corrélation qui 
prend en compte tous les paramètres qui influent sur la conductivité thermique en utilisant des 
données expérimentales élargie et étendues. 

Donc, dans ce chapitre, un modèle mathématique est développé pour prédire la 
conductivité thermique des nanofluides d’Oxydes Métalliques/Eau. Quatre facteurs 
d'influence majeurs sur la conductivité thermique des nanofluides ont été sélectionnés ; et un 
plan d'expérience (Design of Experiments DOE) avec la méthode ANOVA sont utilisés pour 
déterminer les effets de ces facteurs et leurs interactions. Une méthode factorielle complète à 
deux niveaux est considérée pour chaque facteur [2], puis les résultats sont analysés. Cette 
modélisation est basée sur des données expérimentales recueillies dans la littérature, puis les 
résultats prédits ont été comparés à ceux d'autres modèles ainsi qu'aux résultats 
expérimentaux. 

La nouveauté dans le présent travail est l’utilisation seulement des facteurs 
significatifs, nécessitant ainsi des données expérimentales et du temps de calcul moindres. En 
outre, les conditions de fonctionnement sont à basse température et à fraction volumique 
moyennes, qui diffèrent de la majorité des modèles empiriques et statistiques utilisés (qui 
sont en générale appliqués pour des plages de températures plus élevées). Par conséquent, le 
but de cette étude est de trouver un nouveau modèle de prédiction de la conductivité 
thermique de trois types de nanofluides pour des basses températures et des fractions 
volumique moyenne. 


IIL.2 Position du problème 


IIL.2.1 Caractéristiques des matériaux 

Les conductivités thermiques des fluides ordinaires sont plutôt faibles, donc afin 
d'améliorer le transfert de chaleur des fluides, différents types de nanoparticule de 
conductivité thermique élevée sont utilisées. En raison des caractéristiques élevées de 


conductivité thermique et de performance de transfert de chaleur, les nanofluides sont 
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devenus un choix très attrayant pour des applications pratiques dans les systèmes d'ingénierie 
thermique. 

Les données expérimentales utilisées dans cette étude sont extraites de la littérature 
[31-35] ; et les détails sont présentés dans le tableau IIL.1. 

Les nanofluides utilisées dans le présent travail sont les nanofluides d'oxyde 
métallique, ce sont les types les plus utilisées dans le transfert de chaleur en raison de leur 
faible coût. Les nanoparticules d'oxyde métallique sont produites par des méthodes 
conventionnelles, elles sont plus stables que les nanoparticules métalliques [36]. Trois types 
de nanoparticules d'oxyde métallique sont étudiés : les Oxydes d'Aluminium (Alumine) 
AbLO;3, Oxydes de Titane (Titania) TiO2 et Oxydes de Cuivre CuO [37-39]. Les propriétés 
thermophysiques correspondantes du fluide de base (Eau) et des phases solides sont données 


dans le tableau IIL2 [40]. Ce dernier présente les propriétés thermophysiques utilisées dans 


cette étude. 


Tableau IIL.1: Caractéristiques des nanofluides à base d'eau adoptées dans ce travail. 


Nanoparticules | Fraction volumique | Diamètre de NP Température Chercheurs 
(%) (nm) (K) 
CuO 0-18 29, 28.6 température ambiante Aminian [31] 
ALO; 14 20, 25, 28 
CuO 0-5 28.5 Différente Kazemi-Beydokhti et al. [32] 
ALLO; 0-5 20, 25 température 
TiO; 0.2-2 21 288-308 Duangthongsuk et al. [33] 
TiO; 0-8 21 288-338 Yiamsawasd et al. [34] 
ALO; 14 30+10 274-313 Longo et al. [35] 
TiO; 1-6 30-50+10 


Tableau IIL.2: Propriétés thermophysiques de l'eau et des nanoparticules sélectionnées. 


Propriétés thermophysiques u P C K 
[kg/m.s] [Kg/m°] [J/kg. K] [W/m.K] 
Fluide de base Eau (H,0) 0.000912 997.047 4184 0.6063 
Alumine (ALO;) / 3970 765 40 
D Oxydes de Titane (TiO:) / 4250 686.2 8.9538 
Oxydes de Cuivre (CuO) / 6500 535.6 20 


IIL.2.2 Méthodologie 
La méthode de plan factoriel complet est une technique puissante, rapide et 


avantageuse, utilisée pour obtenir la meilleure optimisation. Cette conception nous permet de 
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déterminer les effets sur la réponse de chacun des niveaux de facteurs. De plus, elle prend en 
considération toutes les combinaisons entre niveaux de facteurs ainsi que les effets des 
interactions entre facteurs. Cette technique est utilisée pour réduire le nombre d'expériences, 
le temps et le coût global du processus avec une meilleure réponse. L'importance de cette 
méthode est montrée par différentes études utilisant des conceptions expérimentales [41, 42]. 
IIL.2.3 Choix de la réponse, des facteurs et du type de plan correspondant 

Les différents facteurs pris en compte dans ce travail sont les plus influents sur la 
conductivité thermique des nanofluides ; à savoir la fraction volumique de nanoparticules, la 
température, le diamètre des nanoparticules et les types de nanoparticules. La grandeur 
d’intérêt (réponse) est alors le rapport de conductivité thermique du nanofluide sur celle du 
fluide et on la note comme : (kng " kr). On constate que les trois premiers facteurs sont des 
facteurs continus tandis que le dernier facteur est un facteur discret qui dépend du type de 
nanoparticules. 

En ce qui concerne la méthode des plans d’expériences, on a utilisé les plans factoriels 
à niveaux mixtes avec le modèle correspondant de premier ordre. 

IIL.2.4 Modèle mathématique 

Dans cette étude, une méthode des plans d’expériences est adoptée pour optimiser la 
conductivité thermique des nanofluides. La technique (DOE) est réalisée pour identifier les 
variables clés qui affectent le taux de conductivité thermique TCR (ie. k,F/kf) et pour 
déterminer à quels niveaux ces facteurs doivent être conservés dans le but de réduire les 
erreurs de mesure. La méthode de plan factoriel complet est utilisée afin de permettre l’étude 
des effets causés par les facteurs indépendants ainsi que les interactions entre ces facteurs [43, 
44]. 

L'utilisation des variables centrées réduites présente l'intérêt de pouvoir généraliser la 
théorie des plans d'expériences quels que soient les facteurs ou les domaines d'études retenus. 
Remplacer les variables naturelles par les variables codées va permettre d'avoir pour chaque 
facteur le même domaine de variation (entre le niveau bas codé -1 et le niveau haut codé +1). 

Le tableau IIL3 montre les facteurs sélectionnés et leurs domaines de variation. Le 
domaine de variation des facteurs doit permettre de couvrir la zone optimale requise de la 
conductivité thermique des nanofluides, qui présente le domaine réel d'utilisation des facteurs 
dans cette partie. Le choix de l’intervalle de chaque facteur n’est pas aléatoire, il est selon les 


données expérimentales existant dans cette zone de transfert de chaleur. 
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Le logiciel JMP (John's Macintosh Product), conçu pour faciliter au maximum la 
réalisation de ces opérations en assurant le suivi complet des calculs des plans d’expériences, 
est utilisé. Les données expérimentales sont données dans le tableau IIL.2 

La matrice d'expériences récapitulant l'ensemble des essais comprend plusieurs 
colonnes : la première identifie les essais, les suivantes indiquent les états des facteurs 
exprimés en "unités codés" ; chaque colonne étant représentative d'un facteur. La dernière 
colonne indique les réponses obtenues pour chaque essai. Les lignes inférieures rappellent les 
valeurs des niveaux en unités naturelles. Le tableau IIL.4 présente une matrice d'expériences 
pour 3 facteurs à 2 niveaux. 

Un plan factoriel complet à deux niveaux pour trois facteurs indépendants est 
considéré : la fraction volumique (X:), la température (X2) et le diamètre des nanoparticules 
(X3); avec deux niveaux factoriels codés comme —1 et +1 (niveaux bas et haut 
respectivement). Les types de nanoparticules (X4) sont codés comme L pour (A1O3), L2 pour 
(TiO)) et La pour (CuO). Un plan factoriel complet 2 (c'est-à-dire que trois facteurs varient 
sur deux niveaux) pour trois types de nanoparticules est utilisé. Par conséquent, un total de 
vingt-quatre essais expérimentaux est suggéré pour analyser l'interaction de chaque niveau 
(tableau IIT.4). Le rapport (kr kr) est considéré comme le facteur de réponse (facteur 
dépendant). Les résultats sont analysés et tracés avec des intervalles de confiance à 95% à 
l'aide du logiciel JMP. 

IL.3 Résultats et discussion 


La corrélation obtenue prend la forme d'un polynôme et le modèle mathématique 
généré par la méthode des plans d’expérience est de la forme suivante : 
Y = 


L2 — —0.2098875 + X, * (X: * 


1.2871875 + 0.04789 * X, —-0.214930833* X, -0.0628004 * X, + 
| L3 — —0.0198875 


L1 — 0.228885 


0 a a 4 
X4 
L1 — 0.036285 


-0.013795) + X, * (X;) * 0.064382916) + X, * (X;) Le — —0.0118625 + X, * (X,) 


L3 — —0.0244225 
a12 A3 


414 X4 
L1 — —0.237329166 L1 — 0.0559529166 
E — 0.2153258366 + X, * A) L2 — 0.058566666 +X;, * 
L3 — 0.022003333 L3 — —0.114519583 


Q24 X4 434 X4 
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L1 — —0.055947916 L1-— —-0.0313325 

X3 * (X4)4L2 — —0.071206666 + X,* {X * (X4) L2 — 0.0226025 IL. 1) 
L3 — 0.1271545833 L3 — 0.0087300 
| 


da Xa 4124 X4 


Le modèle obtenu a pour forme générale : 


Y = 1.083716666 + (0.032505833 + C,)X; + (0.0061325 + C>)X2 + (—0.008745416 + C3)X3 + 
(—0.005212916)X:X3 + (Ca) X1X2 + (C:)X1X3 + (0.00418375 + Cé)X1X2X3 + Cy (IIL2) 


@—0.025 T—-296.5 dp—0.000000025 


Avec : X, = ——— = — 
1 0.015 ” “2 Le. 274 0.000000004 


, €t Ci, C2, Ca, Ca, Cs, Ce et C, sont des 


constantes dépendant du type de nanoparticule et sont déterminées par leurs propriétés 
physiques. Les valeurs de ces constantes pour chaque nanoparticule sont données dans le 


tableau IILS. 


Tableau IIL.3: Les facteurs sélectionnés à différents niveaux du plan factoriel. 


Facteurs niveaux bas (-1) niveaux haut (+1) 
(X1) fraction volumique (%) 1 4 
(X2) Température (K) 295 298 
(X3) diamètre des nanoparticules (nm) 21 29 


Les modèles partiels de chacun des trois types de nanofluides sont : 


ALO;: Ÿ = 1.0640665 + 0.036458749 X, + 0.00429625 X, + 0.00165375 X; (IIL.4) 
— 0.005212916 X,X; + 0.00021375 X,X, + 0.000518333 X,X; 


+ 0.00179125 X,X2X3 


TiO: : Ÿ = 1.080645 + 0.036027499 X, + 0.000395 X, — 0.0039375 X; (IIL.S) 
— 0.005212916 X,X; + 0.0076375 X,X, + 0.004548333 X,X, 


— 0.00080875 X,X2X; 


CuoO : Ÿ = 1106438749 + 0.02503125 X, + 0.01370625 X, — 0.023656249 X, (IIL.6) 
— 0.005212916 X,X3 — 0.00785125 X,X, — 0.005066666 X,X; 


+ 0.0115675 X,XX3 
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Tableau IIL.4: Matrice d'expérience. 


1.100320 


1 -l 1 L1 1.093 


1 1 1 L1 1.105400 


1 -1 -1 L2 1.096330 


1 1 -1 L2 1.130000 


1 -1 1 L2 1.111920 


1 1 1 L2 1.115060 


1 -1 -1 L3 1.154 


1 1 -1 L3 1.16 


1 -1 1 L3 1.088200 
1 1 1 L3 1.110300 


Modélisation de la conductivité thermique effective des nanofluides par l’utilisation des 
plans d'expériences à Plan factoriel complet 


Tableau IIL.S: Valeurs des constantes. 


Constantes 
Nanoparticules 
Cx10 | C:x10 | C3x10* | Cux10 | Csx10 4 | Cix10* | C,x10* 
ALO; 39.52 -18.36 103.99 | 2.13 5.18 -23.92 |-196.5 
TiO> 35.21 -57.37 | 45.16 76.37 45.48 -49.92 | -30.71 
CuO -74.74 75.73 -149.1 78.51 -50.66 | 73.85 227.22 


IIL.3.1 Test du modèle proposé 

Le modèle obtenu est testé pour validation par la comparaison du rapport de 
conductivité thermique (k,/k;) des valeurs prédites avec celles réelles, tout en tenant 
compte de l'erreur standard (Figure II.1). Il est clairement montré que toutes les valeurs de 
TCR sont dans la ligne centrale (valeurs abscisses égale aux valeurs des ordonnées) et que la 
valeur de RMSE (Root Mean Square Equals) est égale à seulement (0,0074). De plus, la 
proportion de la variation autour de la moyenne exprimée par le modèle (linéaire ou 
polynomiale) représenté par (R?), a pour valeur égale à (0,99). Ceci indique que le modèle de 
TCR fournit une bonne corrélation et que les valeurs prédites et expérimentales s’adaptent 


parfaitement. 


TCR Actual 


TCR Predicted 


Figure III. 1: les valeurs actuelles de TCR en fonction des valeurs prédites. 


(P value=0.0013, R square=0.99, RMSE=0.0074) 
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Le tableau (IIL.6) donne les paramètres estimés avec leurs erreurs types ainsi que le t- 
test correspondant à chaque valeur de variable indépendante. Les résultats obtenus montrent 
que toutes les erreurs types sont de très faibles valeurs (<0,003) et que les probabilités 
correspondantes indiquent que les résultats sont statistiquement significatifs représentés par 
(). 

L'analyse de variance (ANOVA) pour un plan factoriel complet du modèle 
polynomial a été utilisée afin d’évaluer l’effet de chaque paramètre sur le modèle de TCR. La 
méthode statistique d’'ANOVA est utile, non seulement pour déterminer l'influence de 
n’importe quels paramètres apportant une contribution aux résultats expérimentaux mais aussi 
pour interpréter les données expérimentales. Le tableau (11.7) montre les résultats d'ANOVA 
obtenu pour TCR. On remarque que le modèle proposé est significatif donnant une valeur de 
P=0,0013. En outre, le niveau de confiance pour les termes linéaires X, (P<0,0001), X3 
(P<0,0045) et X4 (P<0,0008) sont les plus significatifs sur TCR, tandis que X2 (P<0,0157) est 
toutefois significatif. La plupart des termes carrés et leurs interactions sont également très 
significatives. 

De même, on peut voir à partir du tableau (IL.7) que la fraction volumique des 
nanoparticules est le facteur ayant la plus grande influence sur le modèle de TCR, car elle a 
une valeur F très élevée (F-value = 457.0790). Par conséquent, la fraction volumique est le 
facteur le plus significatif sur ce modèle. Également, le diamètre des nanoparticules a un effet 
significatif sur TCR avec (F-value = 33.0849), ce qui montre que le diamètre est autant 
important. Par contre, la température qui a un effet inférieur à celui des autres paramètres 
principaux avec une valeur F-value de 16,26, a également un effet significatif sur la réponse. 
En ce qui concerne les interactions, on constate que les termes (X2XxX3, X3xX4 et X1XX2XX4) 
ont des effets significatifs sur TCR ; cependant le reste des interactions ont un effet non 
significatif sur la réponse, car la valeur de leurs probabilités (P>0.05). De plus, les valeurs R° 
et R° ajustées de 0,99 et 0,97 respectivement, sont très proches de 1 ; suggérant une forte 
relation entre les variables et qu’il y a une corrélation élevée entre les valeurs expérimentales 


et ajustées. Cela confirme que le modèle obtenu présente une excellente capacité prédictive. 
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Tableau IIIL.6: Prédiction de modèle et les erreurs standards. 


Term Estimation Std Error t Ratio Prob>|t| 
Intercept value a, | 1-08037167 0.00152043 710.57 <0.0001* 
X:(0,01,0,04) 0.03250583 0.00152043 21.38 <0.0001* 
X2(295,298) 0.0061325 0.00152043 4.03 0.0157* 
X3(2,1e-8,2,9e-8) | -0.00874542 | 0.00152043 -5.75 0.0045* 
X4L1] -0.01965042 | 0.00215021 -9.14 0.0008* 
X4IL2] -0.00307167 | 0.00215021 -1.43 0.2263 
XA[L3] 0.02272208 0.00215021 10.57 0.0005* 
X,X*X4[L1] 0.00395292 0.00215021 1.84 0.1399 
XiXx4[L2] 0.00352167 0.00215021 1.64 0.1768 
X1XX4[L3] -0.00747458 | 0.00215021 -3.48 0.0254* 
X2XX3 -0.00521292 | 0.00152043 -3.43 0.0266* 
X,xX4[L1] -0.00183625 | 0.00215021 -0.85 0.4412 
X2xX4[L2] -0.0057375 0.00215021 -2.67 0.0559 
X2xX4[L3] 0.00757375 0.00215021 3.52 0.0244* 
X3xX4[L1] 0.01039917 0.00215021 4.84 0.0084* 
X3xX4[L2] 0.00451167 0.00215021 2.1 0.1039 
X3xX4[L3] -0.01491083 | 0.00215021 -6.93 0.0023* 
X1XX2XxX4[L1] 0.00021375 0.00215021 0.1 0.9256 
XXXXX4[L2] 0.0076375 0.00215021 3:59 0.0238* 
XXXXX4[L3] -0.00785125 | 0.00215021 -3.65 0.0217* 
X1XX3XxX4[L1] 0.00051833 0.00215021 0.24 0.8214 
X1XX3XX4[L2] 0.00454833 0.00215021 2.12 0.1019 
X1XX3XX4[L3] -0.00506667 | 0.00215021 -2.36 0.078 
XX XX X3 0.00418375 0.00152043 2.15 0.0513 
X1XX2XX3xXA[L1] | -0-0023925 0.00215021 -1.11 0.3282 
X1XX2XX3xXA[L2] | -0.0049925 0.00215021 -2.32 0.081 
X1XX2XX 3x X[L3] 0.007385 0.00215021 3.43 0.0264* 


(* Signifiant selon la valeur de probabilité < 0.05) 
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Tableau IIL.7: Analyse de variance ANOVA pour TCR. 


source SS DE MS F-value Prob > F 
Model 0.04270618 19 0.0022477 40.5130 0.0013* 
X:1(0,01,0,04) 0.02535910 î 0.0253591 457.0790 <0.0001* 
X2(295,298) 0.00090258 1 0.0009026 16.2684 0.0157* 
X:(2,1e-8,2,9e-8) 0.00183558 î 0.0018356 33.0849 0.0045* 
x 0.00729494 2 0.0036475 65.7429 0.0009* 
XIXX 0.00067118 2 0.0003356 6.0487 0.0617 
are 0.00065219 î 0.0006522 11.7552 0.0266* 
X2XX4 0.00074922 2 0.0003746 6.7521 0.0522 
X3XX4 0.00280665 2 0.0014033 25.2939 0.0054* 
X1XX2XX4 0.00096015 2 0.0004801 8.6530 0.0352* 
XX XX X 4 0.00037302 2 0.0001865 3.3617 0.1391 
X1XX2XX 3 0.00042009 1 0.0004201 7.5718 0.0513 
XX X2XX XX 4 0,00068150 2 0.0003407 6.1418 0.0603 
Residual 0,00022192 4 0.0000555 _ _ 
Total 0.04292811 23 _ _ = 
R? 0.99483035 = _- = _- 
Adj R° 0.97027453 = = = _ 


(* Signifiant selon la valeur de probabilité < 0.05, SS somme des carrés, DF degré de liberté, 
MS carré moyen). 


0.02272 296.381 
T&) 


Figure IIL.2: Prédiction de TCR par rapport à X1(9@), X2(T), X3(dp) et X4 (Nb). 


[Lignes rouges : valeur prédite de TCR ; Lignes bleu : 95% de confiance sur la moyenne de TCR] 
La figure (IL2) montre les profils de la réponse prédite TCR par rapport aux quatre 
facteurs: (4, T, dp et Np), codés X1, X2, X3 et X, respectivement. 
Ces résultats montrent qu'il existe une corrélation entre la fraction volumique (4) et la 


température (T) et TCR, car ces deux facteurs ont un effet positif sur la réponse. La valeur de 
TCR augmente linéairement et atteint des valeurs maximales quand les facteurs X, et X; sont 


au maximum. En revanche, le diamètre des nanoparticules (dp) a un effet négatif sur la 
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réponse. Cela signifie que la réponse diminue lorsque ce facteur passe du niveau bas au 
niveau haut. En conséquence, il peut être confirmé que CuO est la nanoparticule prééminente 


utilisée pour augmenter la réponse TCR. 


g 1.156126 
OÙ [1.13748, 
F 117477 
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Figure IIL.3: Profiles de prédiction aux paramètres optimaux. 


La maximisation de la désirabilité est effectuée pour obtenir des paramètres de facteur 
optimaux pour un TCR maximum. Le profil de la figure (IL3) affiche les paramètres 
optimaux suivants : 0,042 pour la fraction volumique ; 294,73 K pour la température ; 
21,27 nm pour le diamètre des particules et CuO pour le type de nanoparticules, donnant ainsi 
une réponse estimée égale à 1,162. 

IIL.3.2 Validation du modèle proposé 

Afin de valider le modèle proposé du rapport de conductivité des nanofluides (TCR), 
une comparaison avec des résultats expérimentaux et théoriques est effectuée. Les figures 
IIL.4-TIL.9 montrent les effets de différents paramètres sur le rapport de conductivité thermique 
de l'ensemble de données expérimentales des nanofluides (tirées de la littérature) par rapport 


au modèle obtenu dans notre étude. 


| Working fluid:CuO-Water, T=295 K | (a) | Working fluid:CuO-Water, T=298K | @) 


M Das et al, dp=28.6nm 

©  Khedkar et al, dp=25nm 
— Proposed model, dp=28.6nm 
— - Proposed model, dp=25nm 


© Mintsa et al, dp=29nm 
— Proposed model, dp=29nm 
À Lee et al, dp=23.6nm 
— - Proposed model, dp=23.6nm 


Thermal conductivity ratio (k,,g/kf) 
Thermal conductivity ratio (k,,gkf) 


Particle volume fraction (%) Particle volume fraction (%) 


Figure IIL.4: Effets de la fraction volumique des nanoparticules sur TCR du nanofluide 


CuO/Eau par rapport au modèle proposé pour différents diamètres : (a) 295 K et (b) 298 K. 
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Les effets des fractions volumiques sur le TCR prédit avec différents diamètres de 
nanoparticules de nanofluides CuO-Eau et à des températures de 295 K et 298 K, sont donnés 
dans les figures IIL.4 (a) et (b) respectivement. Comme on peut le voir sur ces figures, le 


modèle proposé est en très bon accord avec les données expérimentales [31, 45, 46]. 


[ Working fluid : Ti0,-Water, T=298 K b) 


| Working fluid : A1,03-Water, T=298 K | (a) 


M  Yiamsawasd et al, dp=21nm 

À  Duangthongsuk et al, dp=21nm 
— Proposed model, dp=21nm 

© Masuda et al, dp=27nm 
— = Proposed model, dp=27nm 


À Zhang et al, dp=20nm 
—— Proposed model, dp=21nm 
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Figure IIL.S: Effets de la fraction volumique des nanoparticules sur TCR du nanofluide de (a) 


AbLO;/Eau et (b) Ti0,/Eau par rapport au modèle proposé pour différents diamètres. 


La figure IILS (a) montre la capacité de prédiction de TCR du modèle obtenu pour un 
nanofluide Al203/Eau à 298 K, comparé aux données expérimentales [32]. Ainsi, le modèle 
proposé dans cette étude a une prédiction précise du rapport de conductivité thermique et ceci 
pour différentes fractions volumiques de particules. 

De même, la figure 5 (b) illustre un bon accord de TCR du nanofluide TiO:/Eau du modèle 
proposé par rapport aux données expérimentales de [33, 34, 45]; et ceci aux mêmes 
conditions de température, de fraction volumique et pour différentes diamètres des 


nanoparticules. 


Working fluid: CuO-Water, ® =3.3 % Working fluid: A1303-Water, ®=3.3 % (b) 


© Mintsa et al, dp=29nm 
— Proposed model, dp=29nm 
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Figure IIL.6: Effets de la température sur TCR du nanofluide de (a) CuO-Eau and (b) ALO;- 


Eau, par rapport au modèle proposé. 
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Les effets de variation de température (T) sur le rapport de conductivité thermique des 
nanofluides de CuO/Eau et Al203/Eau pour des nanoparticules de diamètre de 29 nm sont 
présentés respectivement dans les figures IIL6 (a) et IIL6 (b). Les prédictions de TCR du 
modèle proposé sont dans une précision acceptable pour des données expérimentales 


disponibles [31]. 


Working fluid:CuO-Water, T=295 K 


© Mintsa et al, dp=29nm 
— - Hamilton and Crosser (1962) 
- + Koo and Kleinstreuer (2004) 
— - Liand Peterson (2006) 
— - Chonet al (2005) 
--- neural model (2016) 
= Proposed model 
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Figure IIL.7: Comparaison des résultats des modèles expérimentaux et théoriques de TCR de 


nanofluide CuO/eau avec le nouveau modèle proposé, en fonction de la fraction volumique. 


La figure IIL.7 présente les résultats de TCR de nanofluide CuO/Eau à température 
ambiante de notre modèle comparés aux résultats expérimentaux de référence [31] ainsi 
qu'aux résultats des corrélations développées par Chon et al [47], Koo & Kleinstreuer [48], Li 
& Peterson [49], Hamilton—Crosser [50] ainsi que le modèle Neuronal d’Aminian [31]. 
Comme on peut le voir, le modèle proposé dans cette étude offre un meilleur accord avec les 
valeurs expérimentales que les résultats obtenus par les modèles Li et Peterson, Koo & 
Kleinstreuer et Neural. De même, la corrélation développée par Chon et al [47] correspond 
toutefois parfaitement aux données expérimentales. De plus, on remarque que la relation 
Hamilton—-Crosser donne de bons résultats pour des petites valeurs de fraction volumique de 


particules. 
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| Working fluid: CuO-Water, ® =3.3 % | 
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Figure IIL8.: Comparaison des résultats des modèles expérimentaux et théoriques de TCR de 


nanofluide CuO-eau avec le nouveau modèle proposé, en fonction de la température. 


La variation du rapport de conductivité thermique des nanofluides en fonction de la 
température obtenue par notre modèle est comparée aux résultats de diverses corrélations 
ainsi qu'à l'ensemble de données expérimentales de Mintsa et al [31], et ceci pour CuO/Eau de 
29 nm de diamètre et 4% de fraction volumique (figure IIL8). Comme on peut le voir 
clairement, le modèle proposé dans notre étude présente un meilleur accord avec les données 
expérimentales. Néanmoins, les deux relations développées par Koo & Kleinstreuer [48] et 
Chon et al [47] montrent des résultats relativement bons. Par contre, le modèle de Li & 
Peterson [49] ne prédit pas correctement la tendance à l'amélioration de la conductivité 
thermique obtenue par Mintsa et al [31]. Il est également intéressant de noter que le modèle 
de Hamilton—-Crosser [50] reste une prédiction acceptable de la conductivité thermique 
effective des nanofluides à température ambiante. 


T=298 K 


M CuO-Water, dp=25nm 
—— Proposed model, dp-25nm 
© Al03-Water, dp=20nm 
— - Proposed model, dp=21nm 

à TiO,-Water, dp=21nm 
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Figure IIL.9: Effets de la fraction volumique des nanoparticules sur le TCR des nanofluides 


par rapport au modèle proposé pour différentes nanoparticules. 
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La figure II. illustre une comparaison entre des valeurs expérimentales de TCR avec 
celles obtenue à partir de notre modèle, l’ensemble des données expérimentales utilisées dans 
cette figure proviennent des références [32, 34, 46]. Les résultats prédits par notre modèle 
sont de bonne précision et présentent une bonne corrélation par rapport aux résultats 
expérimentaux. 

IL. 4 Conclusion 

Dans cette étude, une approche par un plan factorielle complet a été réalisée afin 
d’estimer l'amélioration de la conductivité thermique des nanofluides d’Oxyde 
Métallique/Eau avec des nanoparticules d'ALO3, TiO; et CuO. Un nouveau modèle évaluant 
le rapport de conductivité thermique (TCR) des nanofluides Oxyde/Eau a été développé en 
utilisant des données expérimentales abondantes et étendues. Ce modèle est applicable dans 
une plage de fractions volumiques de particules de 1% à 4%, à de basses températures entre 
295 K et 298 K et pour des diamètres de particules de 21 nm à 29 nm. Les facteurs affectant 
le rapport d'amélioration de la conductivité thermique sont identifiés et analysés par la 
méthode ANOVA. Ces facteurs sont : la fraction volumique, la température, le diamètre et le 
type des nanoparticules. Les résultats obtenus montrent que la fraction volumique des 
particules a l'effet le plus significatif sur l'amélioration de la conductivité thermique des 
nanofluides par rapport aux autres facteurs. Une augmentation de la fraction volumique et de 
la température des particules possède un effet favorable sur l'amélioration de la conductivité 
thermique des nanofluides, tandis qu’une taille de particule plus grande tend à la réduire. On 
constate que la fraction volumique des nanoparticules a aussi un grand impact sur les 
propriétés thermo-physiques des nanofluides. Aïnsi, la concentration optimale de 
nanoparticules est essentielle pour déterminer dans la formation de nanofluides qui peut 
améliorer les performances des systèmes thermiques. Enfin, l'évaluation des résultats montre 
une excellente corrélation entre les réponses de TCR prédits et observés, indiquant une bonne 
concordance entre le modèle obtenu et les données expérimentales. Par ailleurs, les résultats 
obtenus confirment la capacité du modèle proposé dans notre étude à prédire la conductivité 


thermique des nanofluides avec une grande précision. 
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Chapitre IV, 


Application des Réseaux de Neurones Artificiels (RNA) pour prédire les propriétés 
thermophysiques de nanofluide (A1203 /Eau) 


IV.1 Introduction 

Comme nous l'avons vu auparavant, plusieurs paramètres peuvent jouer un rôle dans 
l'efficacité des nanofluides comme la concentration, le type ainsi que la forme des nanoparticules, 
les propriétés du fluide de base, la température, etc. La majorité des études sur les propriétés des 
nanofluides ont accordé une attention particulière à ces propriétés, en particulier en ce qui concerne 
l'influence de la viscosité dynamique et de la conductivité thermique sur la capacité de transfert de 
chaleur [1]. 

De nombreuses analyses expérimentales [2-6] et théoriques [7-11] des propriétés 
thermophysiques des nanofluides ont été menées par de nombreux scientifiques. En général, la 
température et la fraction volumique sont considérées comme des facteurs importants dans le 
développement des modèles de calcul des propriétés thermophysiques, en particulier pour le calcul 
de la conductivité thermique et la viscosité [12-15]. L'augmentation de la température conduit à 
une augmentation de conductivité thermique et une réduction de la viscosité dynamique; tandis 
qu'une concentration plus élevée en phase solide entraîne une augmentation à la fois de la 
conductivité thermique et de la viscosité dynamique. 

Outre la température et la fraction volumique, d'autres paramètres tels que la forme et le type 
de particules ainsi que le rapport de mélange influencent également les propriétés thermophysiques 
des nanofluides. La valeur précise de ces dernières sous divers paramètres ne peut pas être prédite 
facilement. De plus, en termes de modèles de suspensions particules-liquides classiques, le 
mécanisme de transport de chaleur pour les nanofluides est difficile à comprendre. 

Actuellement, la plupart des modèles disponibles sont applicables aux nanofluides avec une 
concentration plus faible de nanoparticules (de 0,01 à 0,5% en volume) et une dispersion uniforme. 
Cependant, ces modèles ne présentent pas de corrélation entre les propriétés thermophysiques et la 
forme des nanoparticules. En plus, la mesure des propriétés thermophysiques des nanofluides pour 
une température, une concentration, une taille et une forme spécifiques nécessite beaucoup de 
temps, d'énormes données et d'argent. Afin de surmonter ces problèmes, de nombreux chercheurs 
ont suggéré différentes formes de calcul ou de modèles empiriques pour la prédiction des propriétés 
thermophysiques de ces nanofluides. Et le moyen le plus simple et le plus efficace pour déterminer 
ces propriétés pour divers nanofluides est l’utilisation des approches informatiques douces telles que 
les Réseaux de neurones artificiels (RNA) [16]. 

De ce fait, on constate que la majorité des travaux de recherches montrent que les 
paramètres et mécanismes affectant la conductivité thermique des nanofluides ne peuvent pas être 
pris en compte simultanément par des modèles traditionnels de suspensions de particules [17]. En 


conséquence, plusieurs nouveaux modèles utilisant des Réseaux de Neurones Artificiels Améliorés 
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RNAA (EANN) ainsi que des approches Adaptatives du Système d'Inférence Neuro-Floue (ANFIS) 
ont été proposés ; comme les travaux d'Alrashed et al. [18, 19] qui donnent des résultats d’une 
précision appropriée. Tandis-que Hosseini et al. [20] ont présenté un modèle de conductivité 
thermique des nanofluides en fonction de la conductivité thermique des nanoparticules prenant en 
compte le fluide de base, la température, la fraction volumique, la taille des nanoparticules ainsi que 
les propriétés de l’enveloppe inter-faciale et leurs résultats montrant une bonne compatibilité avec 
les données expérimentales ainsi qu’avec d'autres modèles. Shahsavar et al. [21] qu’ont à eux, ont 
exploré plusieurs expériences avec la méthode des Réseaux de Neurones Artificiels (RNA) afin de 
modéliser la conductivité thermique ainsi que la viscosité des nanofluides. Bagherzadeh et al. [22] 
ont développé l'approche améliorée EANN pour le calcul de la conductivité thermique des 
nanofluides, et les résultats obtenus étaient conforme et ceci en utilisant moins de temps de calcul. 
Sadeghi et al. [23] ont utilisé la technique RNA pour vérifier les résultats expérimentaux des 
caractéristiques thermiques (énergie et exergie) du nanofluide Oxyde de Cuivre/Eau dans un capteur 
solaire à tube sous vide et ceci via les méthodes de perceptron multicouche (MLP) et en fonction de 
la base radiale (RBF). Les auteurs ont observé que le taux d’erreur maximal des modèles MLP était 
inférieur à celui du modèle RBF. Les études numériques (par approche Eulérienne à deux phases), 
avec la technique RNA (algorithme d'apprentissage de rétro-propagation de Levenberg—-Marquardt), 
et expérimentales de Naphon et al. [24], montrent que l'écart moyen des caractéristiques thermiques 
des nanofluides de Titane/Eau dé-ionisée obtenu entre les données calculées et prévues était de 


1,25% dans tous les cas. 


Par conséquent, dans ce chapitre, une modélisation par la méthode des Réseaux Neuronaux 
Artificiels est utilisée pour prédire et évaluer la conductivité thermique effective et la viscosité 


dynamique des nanofluides Al:03/Eau. 


IV.2 Réseaux de Neurones Artificiels 

La procédure de fonctionnement des réseaux de neurones artificiels est brièvement 
expliquée ci-dessous. 

Un réseau neuronal artificiel est une approche prédictive précise et puissante qui s'applique à 
la modélisation des systèmes d'ingénierie et à la reconnaissance des modèles. Un avantage majeur 
de l'utilisation des réseaux de neurones est que ce modèle peut simplement apprendre la relation 
entre les variables d'entrée et de sortie. Les modèles de réseaux neuronaux se révèlent extrêmement 


précis dans des applications telles que le contrôle de processus, la modélisation, la simulation et 


l'identification de système [25]. 
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En outre, la technique RNA a été largement utilisée pour prédire les systèmes non linéaires 


en raison des avantages notables qui sont la haute précision, les faibles coûts et temps [26]. 


Des études récentes montrent que les modèles de prédiction de la conductivité thermique 
effective et la viscosité dynamique des nanofluides par réseau de neurone et plus particulièrement 
le Perceptron Multicouches (MLP), qui est un type de réseaux de neurones à fort potentiel prédictif, 


donnent de très bonnes performances [27-32]. 


IV.2.1 Type de modèle RNA 

Les réseaux de neurones MLP (Perceptron multicouches ou en anglais the multilayer 
perceptron) sont des réseaux artificiels de type "feed-forward", c'est-à-dire à propagation directe. 
Une architecture de réseau de neurones largement utilisée dans le domaine des approches de 
données pilotées. Les réseaux de neurones MLP sont capable de créer un mappage entre les données 
d'entrée et de sortie via un algorithme d'apprentissage pour réduire l'erreur entre la sortie réelle et la 
sortie réseau [29, 31]. 

Dans cette étude, un réseau MLP est adopté pour prédire la conductivité thermique et la 
viscosité du nanofluide d’Alumine/ Eau. Les neurones sont groupés en couches, la première et la 
dernière sont appelées respectivement couche d'entrée et couche de sortie. Entre ces deux couches 
il peut exister une ou plusieurs couches appelées couches cachées comme le montre la figure (1V.1) 


ÉESIR 


y; w) 


Couche L 
(Couche de sortie) 


Couche L-1 
(Couche cachée) 


Couche 2 
(Couche cachée) 


Couche 1 
(Couche d’entrée) 


Figure IV.1: Architecture d'une structure neuronale MLP 
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Avec : 
L est le nombre total de couches, la première est la couche d'entrée, la LieME est la couche de sortie, 


les couches de 2 à L — 1 représentent les couches cachées, N, est le nombre de neurone dans la 


[eme couche avec = 1,2, ..,L. wi; représente le coefficient synaptique de la connexion reliant le 


neurone j de la couche L — 1 au neurone i de la L"Æcouche, le coefficient w{, représente le biais 


ou le seuil du neurone i de la L'”® couche. x; est la i°7€ entrée du réseau, Zi est la sortie du 


neurone i de la couche L et enfin y; est la i®7€ 


0,1.…., Ni. 


sortie du réseau. Avec j = 0,1..,N,_;,i = 


La sortie y(x, w) est une fonction non linéaire du vecteur des variables x, de composantes 
X1, X2, …, Xn €t du vecteur des paramètres w, dont les composantes sont les (n) N, + N; paramètres 


du réseau. 
y = yGw) = g {wo + Di [wi Qt +2 wlx;)]} AV.D 


Avec Wy,,o et w!, qui sont les biais attribués respectivement à la couche de sortie et à la couche 
masquée. Wn, iet wi; désigne la matrice de pondération liant le calque masqué au calque de sortie et 


l'entrée au calque masqué, et f(.) représente la fonction d'excitation (activation) [34]. 
Nous adoptons la forme du MLP de la figure (IV. 1), contenant une seule couche cachée et 
ayant la fonction sigmoïde comme fonction d'activation donné par l'équation (1V.2) et la fonction 


linéaire comme fonction d'activation de la couche de sortie donné par l'équation (1V.3). 


fG) = —— (V2) 


1+e* 


g(x) = x (IV.3) 


IV.2.2 Procédure d’apprentissage de réseaux de neurones artificiels 

Le processus d'apprentissage est une procédure bien établie qui tente d'ajuster les paramètres 
de modèle RNA, c'est-à-dire les biais et les poids de connexion. Au cours du processus 
d'apprentissage, le modèle RNA est doté d'un nombre suffisant de données expérimentales 
contenant à la fois des variables indépendantes et dépendantes des phénomènes considérés. 

Ensuite, un algorithme d'apprentissage approprié est utilisé pour ajuster les paramètres du 
modèle RNA de manière à pouvoir prédire les cibles réelles avec une précision acceptable. La phase 
de formation commence avec des valeurs aléatoires pour les poids et les biais du modèle RNA. Par 
la suite, les signaux numériques de variable (s) indépendante (s) sont envoyés au réseau artificiel et 
sont amenés à le traverser jusqu'à ce qu'ils atteignent la couche de sortie. Enfin, la couche de sortie 


est responsable de la production du ou des signaux de sortie. Un algorithme d'apprentissage est 
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ensuite appliqué pour minimiser la différence entre les valeurs réelles et calculées pour la ou les 
variables dépendantes en régulant les paramètres du modèle RNA. Cet ajustement se poursuit 
jusqu’à ce que l'écart entre les valeurs cibles calculées et réelles atteigne la tolérance prédéfinie. 
Dès que la phase de formation est terminée, les poids et les biais sont ajustés et restent inchangés. A 
ce stade, il est possible d'utiliser l'approche RNA entraînée pour estimer la ou les variables 


dépendantes à partir de nouveaux ensembles de données indépendants [1]. 


‘ 
0 


Données 
d'apprentissage : 
Algorithme 
d'apprentissage 


Données de =: Réseau de neurones 
! validation & \ (MLP) 
‘ detest 


Figure IV.2: Schéma bloc de l'apprentissage d'un réseau de neurones 


IV.3 Méthodologie 


IV.3.1 Données expérimentales 

Dans cette étude, une structure de réseau avec une couche cachée est utilisée pour prédire la 
conductivité thermique et la viscosité du nanofluide (figure IV.3). La fraction volumique, la 
température, le diamètre et la forme des nanoparticules sont saisis en tant que paramètres 
indépendants (Entrées), le rapport de conductivité thermique (k,F/k) et le rapport de viscosité 
(Unr/u) sont considérés comme des paramètres dépendants (Sorties). La Gamme de ces 


paramètres est représentée dans le tableau IV. 1. 


Fraction 
Volumique 


Température 


Diamètre de 
nanoparticule 


Forme de Couche de 
nanoparticle sortie 
CS 


Couche 
d'entrée 


Figure IV. 3: Structure de réseau de neurone utilisé dans cette étude 
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Pour développer le modèle RNA, 210 valeurs de données expérimentales de conductivités 
thermiques effectives et de viscosités du nanofluide Al:03/Eau obtenues à partir de la littérature 
[35] sont utilisées pour concevoir et valider les performances du modèle prédictif. Les détails des 
données ainsi que les propriétés thermophysiques de nanoparticule et du nanofluide utilisés dans 


cette étude sont présentés dans les tableaux 1, 2, 3, 4 et 5 respectivement dans l’annexe B. 


Tableau IV.1: Gammes de paramètres pour le nanofluide AO; / Eau 


Paramètres Gamme 
X1: Fraction Volumique (%) 1-5 
X2: Température (K) 300-320 
X3: Diamètre de nanoparticule (nm) 20-100 
X4: Forme de nanoparticle Lame / plaquette / cylindrique / brique / sphérique 


IV.3.2 Développement du modèle RNA 

La méthodologie du présent travail a commencé avec les données collectées à partir des 
résultats expérimentaux menés dans l'étude précédente [35]. Les données ainsi collectées ont été 
analysées puis les entrées et sorties ont été sélectionnées en conséquence. Les paramètres tels que la 
température, la concentration, la taille et la forme des nanoparticules ont été sélectionnés comme 
entrées, tandis que les propriétés thermophysiques telles que la conductivité thermique et la 
viscosité ont été considérées comme des sorties pour le développement du réseau neuronal. 

L’ensemble de données récupérées ont été subdivisées en deux ensembles, 70% pour 
l'apprentissage et 30% pour la validation. L'ensemble de données de validation à 30%, représentant 
30 points de données sélectionnées au hasard n'a pas été utilisé dans le processus de formation, maïs 
a été utilisé comme un validateur externe. 

Les modèles RNA proposés ont été ajustés par l’algorithme de de Levenberg- 
Marquardt (LM) dans la phase d'apprentissage, en utilisant les ensembles de données 
d'apprentissage divisés, ensuite les modèles ont été testés et validés. 

Dans le présent travail, le logiciel MATLAB R2015b a été utilisé pour le développement des 
modèles. L'algorithme LM est sélectionné pour obtenir la valeur maximale du coefficient de 
détermination (R°). Les modèles proposés sont formés par des algorithmes sélectionnés pour les 
moindres erreurs possibles [16]. L'algorithme LM a donné une meilleure précision et est décrit 
comme un processus itératif qui localise la plus petite fonction multivariée énoncée. 

Pour évaluer les performances de notre modèle, nous avons adopté les indicateurs 


statistiques les plus utilisés dans la littérature, résumés dans le tableau IV. 2 


Application des Réseaux de Neurones Artificiels (RNA) pour prédire les propriétés 
thermophysiques de nanofluide (A1203 /Eau) 


Tableau IV.2: Indicateur d'évaluation des modèles 


Désignation Equation 


L'erreur absolue moyenne (Mean Absolute Error) x 
1 
MAE = % D lil 
1=1 


Erreur quadratique moyenne (Mean Square Error) 14 
MSE = n 2,0" -$> 


Somme de l'erreur quadratique (Sum of Squared Error) , 
SSE = D 'Gi-P 
1=1 


Racine carrée de l'erreur quadratique Moyenne (Root Mean Square Error) Dé 


N 
1 
RMSE = De - | 


_ZaOn = #)? 
DCE 


Le coefficient de détermination (coefficient of determination) 
R?=1 


L'écart relatif moyen absolu (Absolute Average Relative Deviation) 


j" _… 
Yi 


N 

1 

DA 
Da D) 


Où N est le nombre de données expérimentales, y; est la valeur prédite, Ÿ, est la valeur actuelle et ÿ est la moyenne des 


valeurs actuelles. 


Note : 

La MAE est utilisée pour mesurer l’écart entre les prédictions et les mesures. 

La RMSE représente la moyenne des carrés des résidus. 

Rappelons que la réduction de la MAE améliore en moyenne la prédiction, alors que la réduction de 


la RMSE signifie plutôt de limiter les grosses erreurs de prédiction. 
IV.3.3 Choix de l'algorithme d'apprentissage 


L'apprentissage est une étape importante de la prédiction, elle consiste à trouver des valeurs 
adéquates des poids et des biais suivant les conditions initiales choisies. 

Nous adoptons dans cette étude l'algorithme d'apprentissage de Levenberg Marquardt LM), 
qui permet d'obtenir une solution numérique au problème de minimisation d'une fonction souvent 
non linéaire et dépendante de plusieurs variables ; Cet algorithme utilise une méthode dite de 
second ordre consistant à évaluer le gradient de la fonction de coût et de son second dérivé. 

L'algorithme de LM est un standard pour l'optimisation de l'erreur quadratique, connu par 
ses propriétés de convergence rapide pour les problèmes d'approximation de la « fonction de coût ». 
Cet algorithme favori, est surtout évident si un apprentissage très exact est exigé [36-40]. 

Pour approximer une fonction non linéaire par un perceptron multicouche, le choix de 
l'algorithme d'apprentissage est crucial. Dans cette application nous avons écrit un programme en 
langage Matlab basé sur l'algorithme itératif de retro-propagation. Notre choix s'est établi après 


plusieurs tests sur l'algorithme du deuxième ordre de LM, vue le temps d'apprentissage très long du 
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modèle itératif que nous avons développé et la précision des résultats obtenus en utilisant 
l'algorithme de Levenberg Marquardt [33]. 
L'algorithme optimal pour le développement de l'architecture du Réseau Neuronal Artificiel 


est décrit dans l'organigramme de la Figure IV.4. 


Nous avons développé un algorithme qui fait varier le nombre de neurones de la couche 
cachée pendant un nombre d'itérations paramétrable. Le modèle retenu est le premier atteignant une 
erreur quadratique acceptable ou minimale. 

Dans ce but, l'apprentissage a été lancé en normalisant l'ensemble complet de données. 
L'ensemble de données est ensuite divisé en paramètres d'entrée et de sortie (target). La fonction 
d'activation, l'algorithme d'apprentissage et le nombre de couches cachées ont été ainsi sélectionnés. 
Le nombre de deux neurones dans la couche cachée est sélectionné pour initialiser la formation du 
réseau. Le nombre de neurones dans la couche cachée est régulièrement augmenté de cinq à chaque 
répétition d'apprentissage. 

Le nombre optimal de neurones dans la couche cachée est atteints, lorsque le coefficient de 
corrélation est au maximum et l'erreur quadratique moyenne est au minimum [29]. 

L'ensemble de données ainsi obtenu est ensuite comparé à l'ensemble de données 
expérimentales pour obtenir le meilleur ajustement global. L'architecture ainsi obtenue est acceptée 


comme étant le meilleur fournisseur de performances. 


Préparation des données et 


Normalisation 


Définir les entrées et les cibles 


Définir la fonction d'activation, 
l'algorithme d'apprentissage, et 
le nombre de couches 


Définir le nombre de neurones 


Données d'apprentissage et de 
validation 
Calculer les erreurs 


Nombre de Nombre de neurones = 
neurones <30 Nombre de neurones +5 


Sélectionnez la meilleur RNA 
architecture 


Evaluer les données de testet 


validation 


Figure IV.4: Alscorithme de réseau neuronal utilisé dans ce travail 
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IV.4 Résultats et discussion 


IV.4.1 Estimation de la conductivité thermique basée sur le RNA 

Un réseau de neurone artificiel RNA a été développé sur la base des données expérimentales 
de la littérature disponible. Ces données expérimentales ont été divisées en deux groupes comme 
indiqué précédemment. Pour l'apprentissage 70% des données sont prises et 30% ont été utilisés 
pour les tests et validations. 
L'évaluation des modèles à été réalisée en utilisant à la fois les illustrations statistiques et 
graphiques. La méthode d'essai et d'erreur était effectuée pour sélectionner les meilleures valeurs 
pour le nombre de neurones et de couches cachées. Différents nombres de couches et de neurones 
ont été vérifiés et leurs métriques statistiques ont été déterminées. Il a été constaté que le réseau 
avec une (01) couche cachée a donné les meilleurs résultats souhaités. Les valeurs de l’erreur 
absolue moyenne (MAEËE), l'erreur quadratique moyenne (MSE), la somme de l'erreur quadratique 
(SSE), la racine carrée de l'erreur quadratique moyenne (RMSE) et le coefficient de détermination 
(R°?) pour l'apprentissage et la validation sont présentés dans le tableau IV.3 et schématisés dans la 


figure IV. 5. 


Les valeurs MSE, MAE, SSE, RMSE et R° observées par le modèle MLP à différents 
nombres de neurones cachés pour l'ensemble global de données (entrainement + validation) sont 
affichées dans la figure (IV. 5). On remarque que chacun des MSE, MAE, SSE, RMSE prend sa 
valeur minimale à 15 neurones, contrairement à R° qui prend sa valeur maximale au même nombre 


de neurones. Cela signifie que le nombre optimal de neurones dans la couche cachée est de 15. 


Tableau IV.3: Évaluation des performances du réseau de neurones MLP pour les échantillons 


d'apprentissage et de test pour TCR. 


Indicateur Nombre de Neurones 
Apprentissage Test 
5 10 15 20 25 5 10 15 20 25 

MAE 9,22E-03 | 2,32E-02 | 2,46E-03 | 9.62E-03 | 1,50E-02 | 1,61E-02 | 3.17E-02 | 3,12E-03 | 1.26E-02 | 1.94E-02 
MSE 5,99E-04 | 9,45E-04 | 1,18E-05 | 6.78E-04 | 5,03E-04 | 1,98E-03 | 2.11E-03 | 1,62E-05 | 1,78E-03 | 2.22E-03 
SSE 4,50B-02 | 7,09E-02 | 8,82E-04 | 5,09E-02 | 3,77E-02 | 4,94E-02 | 5,28E-02 | 4,05E-04 | 4,45E-02 | 5,55E-02 
RMSE 2.45E-02 | 3,07E-02 | 3,44E-03 | 2,60E-02 | 2,24E-02 | 4,45E-02 | 4,59E-02 | 4,02E-03 | 4.22E-02 | 4,71E-02 
R° 9,30E-01 | 8.82E-01 | 9,99E-01 | 9.19E-01 | 9,48E-01 | 9,30E-01 | 8.82E-01 | 9,99E-01 | 9.19E-01 | 9,48E-01 


104 


Application des Réseaux de Neurones Artificiels (RNA) pour prédire les propriétés 
thermophysiques de nanofluide (A1203 /Eau) 


Train, 
Validation 


: Train, 
Train. Validation 
Validation 


e) R? 


Figure IV.5: Comparaison statistique 


Il peut être compris à partir du tableau IV. 3 que les performances prédictives du MLP pour 
le sous-ensemble d'entraînement sont améliorées et ceci en augmentant le nombre de neurones 
cachés. Toutefois, les performances de l'approche MLP pour l'estimation de l'ensemble de données 
de test ne s'améliorent qu'à 15 neurones cachés et par la suite aucun progrès évident n'est observé. 
On peut conclure que l'augmentation des neurones cachés de plus de 15 ne fait qu'agrandir les 
réseaux MLP et n'a aucun effet positif sur l'amélioration de la précision [1]. Par conséquent, en 
considérant les résultats rapportés dans le tableau IV. 3, une approche MLP avec une seule couche 
cachée ayant 15 neurones est choisie comme la meilleure structure pour l'estimation du rapport de 


conductivité thermique du nanofluide Alumine/Eau. 


Pour simuler le comportement d'un système donné par le réseau de neurones MLP, une 
quantité relativement étendue de données expérimentales est nécessaire. La figure IV. 6 représente 
la variation de Somme de l'erreur quadratique observée (L'erreur entre les prédictions MLP et les 
données réelles du TCR des nanofluides d'Alumine/Eau) en fonction du nombre d'époque. 

Selon le processus d'optimisation, la meilleure architecture trouvée pour le taux de 


conductivité thermique est de 15 neurones dans la couche cachée. 
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Best Validation Performance is 1.62e-05 at epoch 21 


— Train 
— Validation 


Sum Squared Error (sse) 


10 15 
27 Epochs 


Figure IV.6: Apprentissage de RNA utilisant 15 neurones dans une couche cachée. 


Comme mentionné précédemment, la procédure d'ajustement des poids et biais du modèle 
MLP est effectuée par une technique d'optimisation à savoir un algorithme d'apprentissage. Il est 
évident que la technique d'optimisation tente de minimiser la MSE observée ainsi que la SSE par 
une procédure d'itération. Dans la figure IV. 6, le terme d'époque (Epoch) montre le nombre 
d'interactions que l'algorithme d'apprentissage a essayé d'ajuster avec le paramètre MLP ; en 
utilisant la procédure donnée de la section IV.2.2. 

On peut voir sur la figure IV. 6, qu'en augmentant le nombre d'itérations (c'est-à-dire 
l'époque), les poids et biais du modèle MLP convergent vers leur valeurs optimisées, et par 
conséquent, la SSE observée diminue continuellement. La meilleure performance de validation est 


de 1,62x10* et ceci après 21 époques. 


a) b) 


—€— Experimental Data [| —6— Experimental Data 
—+— ANN 7 [[—2— ANN 


20 40 60 d 10 15 
Data Number Data Number 


Figure IV. 7: Valeurs estimées du taux d'amélioration de conductivité thermique 


a) Données d'apprentissage, b) Données de validation. 


La Figure IV.7 montre le contraste entre les valeurs déterminées expérimentalement et les 
valeurs prédites par RNA (ANN) pour l'apprentissage et la validation. On peut voir que la différence 
entre les valeurs déterminées expérimentalement et les valeurs prédites par RNA est assez petite, et 


le modèle est capable de prédire les valeurs de TCR des données invisibles, c'est-à-dire des données 
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de test avec une bonne précision. Il montre également que le modèle s'est bien adapté aux données 


et qu'il n'y a pas de surajustement dans le modèle malgré le nombre élevé d'itérations. 


0.01 


0.005 F 


20 30 40 50 60 70 5 10 15 
Data Number Data Number 


Figure IV. 8: Erreurs de RNA pour TCR 


a) données d'apprentissage, b) données de validation. 


La figure IV. 8 présente les erreurs du modèle MLP (ye» - Ypréd) pour la conductivité 
thermique de nanofluide (A1O:/H20) en fonction du nombre des données expérimentales pour 
l'apprentissage et la validation. On observe que l'erreur maximale pour les données d'apprentissage 
est de 0.9% et celle pour tester les données est de 0.5%. On remarque aussi que les l'erreur pour les 


données d'apprentissage diminuent avec l'augmentation du nombre des données expérimentales. 


a) Training R2=0.99874 b) Validation R2?=0.9861 


ANN Data 


1.1 1.2 | 1.1 1.15 12 
Experimental Data Experimental Data 


Output = 1*Target + 0.0031 Output —— 0.98*Target + 0.024 


Figure IV.9: Graphe des données expérimentales du TCR en fonction des valeurs prédites par le 


meilleur modèle MLP. a) données d'apprentissage, b) données de validation 


La figure IV. 9 montre la comparaison des résultats expérimentaux pour la conductivité 
thermique de nanofluide (Al1:0:/H,0) avec les résultats prédits obtenus en utilisant le modèle MLP 
optimal. Les résultats de l'apprentissage et de la validation sont ainsi présentés. La valeur R° de 
0,99874 dans le processus d'apprentissage montre un coefficient de corrélation élevé. Aussi, comme 
on peut l'observer sur la figure IV.9.a, tous les points de données sont alignés le long de la ligne 


diagonale du graphique, qui montre un modèle de prédiction de haute performance. 
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La Figure IV.9.b exhibe la valeur R° de l'ensemble de données de validation qui est égale à 
0,9861 confirmant autant une précision et une prédiction élevée. 
Cela signifie qu'il existe un excellent accord entre les prédications du modèle MLP et les 


données expérimentales réelles. 


Tableau IV.4: Comparaisons de la précision prédictive de différentes corrélations intelligentes et 
empiriques de TCR expérimental mesuré par Chandrasekar et al. [41] 


(Avec T = 35,5 °C et dp = 36 nm). 


Fraction L'approche intelligente considérée et les corrélations empiriques 

volumique 

d'ALO; Chon et al. Murshed et al. Mintsa et al. Yu and Xie et al. Nan et al. Modèle TCR 

[42] [43] [44] Choi [45] [46] [47] MLP Actuel [41] 
0,01 1,0331 1,0818 1,0175 1,0341 1,0343 1,0292 1.037 1,032 
0,02 1,0545 1,1646 1,0351 1,0672 1,0673 1,0584 1.077 1,074 
0,03 1,076 1,2523 1,0516 1,1023 1,1025 1,0896 1.0963 1,096 
AARD% 0,587 3,406 1,524 0,388 0,386 0,635 0.0019 _ 


Le tableau IV. 4 présente les résultats de l'estimation des données expérimentales TCR 
rapportées par Chandrasekar et al. [41] en utilisant notre modèle intelligent ainsi que diverses 
corrélations empiriques [42—47]. On peut clairement voir que notre modèle MLP proposé prédit les 
données expérimentales avec le plus petit écart à savoir (AARD =0.0019 %). Les corrélations 
proposées par Yu et Choi [45] et Xie et al. [46] qui fournissent relativement de bons résultats avec 


(AARD = 0,39%) et sont ainsi parmi les meilleures corrélations empiriques. 


La figure IV. 10 illustre la comparaison des données expérimentales de Chandrasekar et al. 
et du TCR prédit par notre modèle du nanofluide Al1:0;/Eau en fonction de la fraction volumique 


pour un diamètre de particule de 36 nm et une température de 35,5 ° C [41]. 
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Figure IV.10: Évaluation comparative pour la prédiction de la TCR du nanofluide du modèle MLP 


développé avec des modèles classiques et expérimental mesuré par Chandrasekar et al. [41]. 


Les résultats indiquent que l'augmentation de la fraction volumique opérationnelle de 1 à 3% 
augmente la conductivité relative des nanofluides, les réseaux de neurones MLP sont en parfait 


accord avec les données expérimentales par rapport aux autres modèles. 


IV.4.2 Estimation de la viscosité dynamique basée sur RNA 

Similairement à la modélisation précédente de TCR, la meilleure structure du modèle MLP 
est sélectionnée par des analyses d’essais et d’erreurs. L'évaluation des modèles a été réalisée en 
utilisant à la fois les illustrations statistiques et graphiques. Les résultats de cette procédure ont été 
rapportés dans le tableau IV. 5 et la figure IV. 11. Différents modèles MLP ayant 1 à 25 neurones 


cachés ont été développés, entraînés, testés et leurs performances ont été évaluées. 


109 


Application des Réseaux de Neurones Artificiels (RNA) pour prédire les propriétés 
thermophysiques de nanofluide (A1203 /Eau) 


Train 


Validation 


25 


Train 
Validation 


e) R? 


Train, 
Validation 


Train. 
Validation 


Figure IV.11: Comparaison statistique 


La figure IV. 11 montre la variation de MSE, MAE, SSE, RMSE et R° des données 


d'apprentissage et de validation avec le nombre de neurones dans la couche cachée. Les résultats 


montrent que les valeurs MSE, MAE, SSE et RMSE sont les plus faibles lorsqu'il y a 20 neurones à 


: : : : : 2 : 
la fois en apprentissage et en validation, Contrairement au R°, qui prend dans tous les cas une valeur 


élevée, signifiant que le nombre optimal de neurones dans les couches cachées est de 20. 


Tableau IV.5: Évaluation des performances du réseau de neurones MLP pour les échantillons 


d'entraînement et de test pour DVR. 


Indicateur Nombre de Neurones 
Apprentissage Test 
5 10 15 20 25 5 10 15 20 25 

MAE |5,22E-02 | 9,80E-03 | 1,48E-02 | 9,38E-03 | 3,87E-02 | 2,54E-02 | L0IE-02 | 1,66E-02 | 1,16E-02 | 1,64E-02 
MSE 8,10E-03 | 3,00E-04 | 8,63E-04 | 7,88E-04 | 9,73E-03 | 1,40E-03 | 2,00E-04 | 6,/8E-04 | 4,09E-04 | 5,47E-04 
SSE GAAE-01 |2,23E-02 | 5,31E-02 | 6,31E-02 | 7,/0E-02 | 3,40E-02 | 4,50E-03 | 1,/0E-02 | 1,02E-02 | 1,37E-02 
RMSE | 900E-02 | 1,73E-02 | 2,94E-02 | 2,81E-02 | 986E-02 | 3,74E-02 | 1,41E-02 | 2,60E-02 | 2,02E-02 | 2,34E-02 
R? D91E-01 | 1,00E+00 | 9,99E01 | 9,99E-01 | 9,91E-01 | 99IE-01 | 1,00E+00 | 9,99E-01 | 9,99E.01 | 9,91E-0I 


Les performances du modèle RNA pendant le 


viscosité peuvent être discernées dans la figure IV. 


processus d'apprentissage pour le taux de 


12. Selon le processus d'optimisation, la 
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meilleure architecture trouvée pour le taux de viscosité est de 20 neurones dans la couche cachée. Il 
a été observé aussi que la meilleure performance de taux de viscosité est pour des valeurs de MSE 


égale à 0,0004089 observée à 12 époques. 


Best Validation Performance is 0.0004089 at epoch 12 


Mean Squared Error (mse) 


8 10 
18 Epochs 


Figure IV.12: Apprentissage de RNA utilisant 20 neurones dans une couche cachée. 


Une comparaison de la prédiction du réseau neuronal pour le taux de viscosité avec les 
données expérimentales correspondantes de l'ensemble de données d'apprentissage et de validation, 
est représentée sur la figure IV. 13. Dans la majorité des cas, le taux de viscosité prédit est très 
proche de la valeur observée expérimentalement. La manière dont la précision de la prédiction varie 
en fonction du nombre de points de données expérimentales dans l'ensemble de données 


d'apprentissage est particulièrement intéressante pour ce travail. 


a) b) 


Se 1——— 
—— Experimental Data —€— Experimental Data 
—#— ANN —+— ANN 


10 20 30 40 50 60 70 5 10 15 
Data Number Data Number 


Figure IV.13: Valeurs estimées du rapport de viscosité 


a) données d'apprentissage, b) données de validation 


La figure IV. 14 expose les erreurs de modèle MLP (yep - Ypréa) pour la viscosité dynamique 
relative (DVR) de nanofluide (Al:0;/H:0) en fonction de nombre des données expérimentales pour 


l'apprentissage et la validation. On observe que l'erreur maximale pour les données d'apprentissage 
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est de 6,5% et l'erreur maximale pour tester les données est de 7%. On remarque aussi que les 
erreurs pour les données d'apprentissage diminuent avec l'augmentation du nombre des données 


expérimentales. 


10 20 30 40 50 60 70 5 10 15 
Data Number Data Number 


Figure IV.14: Erreurs de RNA pour DVR 


a) données d'apprentissage, b) données de validation 


a) Training R?—0.9991 : b) Validation R2=0.99912 


2 3 1.2 1.4 1.6 1.8 
Experimental Data Experimental Data 
Output = 1*Target + 6.1e-05 Output —— 1*Target + 0.029 


Figure IV. 15: Graphe des données expérimentales du DVR en fonction des valeurs prédites par le 
meilleur modèle MLP. 


a) données d'apprentissage, b) données de validation. 


La figure IV. 15 évoque la comparaison des résultats expérimentaux pour la viscosité 
dynamique relative (DVR) du nanofluide (Al1203/H:0) avec les résultats prédits en utilisant le 
modèle MLP optimal. Les résultats de l'apprentissage et de la validation sont ainsi présentés. 

La figure IV. 15 montre d'excellents ajustements le long de la ligne d'égalité pour les sous- 
ensembles d'apprentissage et de validation. Les coefficients de régression R° sont également 
présentés avec chaque ensemble de données. La valeur des coefficients de régression proche de 1 


indique une relation étroite entre les données expérimentales et les données estimées par RNA. On 
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peut voir aussi que la valeur des coefficients de régression est supérieure à 0,99 pour tous les jeux 
de données utilisés dans le présent modèle, confirmant une précision de prédiction élevée. 

Plusieurs modèles capables de prédire la viscosité des nanofluides ont été sélectionnés pour 
être comparés à nos résultats. La figure IV.16, compare les résultats de prédiction de la viscosité 
relative des nanofluides (uniquement pour les données de teste), du modèle (MLP) proposé par nos 
Soins avec ceux existants. 

Afin d'évaluer les performances de prédiction des réseaux de neurones MLP pour le 
nanofluide Al:03/Eau avec un diamètre de 36 nm et une faible concentration volumique de 
nanoparticules, la figure IV.16 présente les valeurs des AARD de notre modèle MLP comparé à 


celles d'autres modèles. 


lative 


CEE — Maiga Batch | Chen 
pt ra ( elor | etal 
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[MRMSE) 0,0202 | 11583 | 1,1762 | 1,5319 | 1,5479 | 1,705 


Figure IV. 16: Comparaison du modèle proposé (MLP) avec les modèles existants pour la 


prédiction de la viscosité relative des nanofluides. 
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En raison de la plus faible valeur de l'écart relatif absolu moyen (AARD = 0,013%) de notre 
réseau de neurones MLP par rapport aux autres modèles tels que Heïdari et al. [48], Maiga et al. 
[49], Einstein [50], Batchelor [51], Chen et al [52] et le modèle de Brinkman [53]. Nous 
démontrons ainsi que les réseaux de neurones MLP peuvent être utilisés avec une précision 
satisfaisante pour la prédiction de la viscosité des nanofluides A1,0; / Eau. 

Les valeurs de RMSE de notre model ainsi que des autres modèles sont également comparés 
dans la figure IV.16.B. 

Le modèle MLP proposé par nos soins, montre une amélioration plus claire avec (RMSE = 
0,0202) par rapport aux modèles RNA 1[31, 48], ainsi que par rapport à d'autres modèles de 
prédiction de viscosité [49, 51, 53]. Par exemple, le présent modèle a réduit l'erreur de prédiction 
(RMSE) de 1,705 à environ 0,0202 par rapport à Brinkman, indiquant une amélioration de 84-fois 


dans la prédiction de la viscosité. 


IV.S Conclusion 

Dans la présente partie du travail, une procédure systématique basée sur le Réseau Neuronal 
Atrtificiel RNA est développée pour prédire la conductivité thermique et la viscosité de nanofluides 
(ALO;/ Eau). La fraction volumique, la température, le diamètre et la forme des nanoparticules sont 
saisis en tant que paramètres indépendants (Entrées), le taux de conductivité thermique (k,/k;) et 


le taux de viscosité (1/17) Sont considérés comme des paramètres dépendants. 


Une approche MLP avec une seule couche cachée est choisie comme la meilleure structure 
pour l'estimation du rapport de conductivité thermique et de viscosité dynamique du nanofluide 
Alumine/Eau. Le nombre optimal de neurones dans les couches cachées sont 15 et 20 neurones pour 
le taux de conductivité thermique et de viscosité dynamique du nanofluide considéré 


respectivement. 


Il ressort clairement de ces résultats que toutes les valeurs de (ky / kF) se trouvent dans la 
ligne centrale, avec des valeurs de l'erreur quadratique moyenne (MSE) de 0,000018 pour 
l'apprentissage et de 0,0000162 pour la validation. De plus, le modèle (k,/k;) fournit une bonne 
corrélation entre les valeurs prédites et expérimentales avec des valeurs de R° égales à 0,998 et 
0,998 pour l'apprentissage et la validation. De même, une excellente corrélation est observée entre 
la réponse prédite et observée pour (lhF/ur). Ainsi, des valeurs plus élevées de R2 = 0,999 et de 
faible valeurs de MSE = 0,000788 pour l'apprentissage ; et R2 = 0,99 et MSE = 0,000678 pour la 


validation, donnant ainsi un bon accord entre les valeurs du modèle et les données expérimentales. 
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De plus, l'approche RNA sélectionnée prévoit un apprentissage avec un AARD de 0.013 % pour 
(Unr/H) et de 0.0019 % pour (k,F/kr). 


Ces résultats montrent une excellente corrélation entre les réponses prédites et observées de 


(knr/Kr) et de (lnF/U), prouvant ainsi une bonne concordance entre les valeurs du modèle 


proposé et les valeurs des données expérimentales confirmant ainsi la validité et la faisabilité de la 


méthode adoptée. 


La comparaison avec les diverses corrélations empiriques confirme que le modèle proposé à 


des performances supérieures dans la prédiction du TCR et du DVR des nanofluides Alumine/Eau 


et peut être considéré comme un outil pratique pour la tâche considérée. 
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Conclusion 
générale 


Conclusion générale 


En l’absence de l’outil expérimental, il devient impératif d’utiliser l’outil informatique 
pour palier à cette carence. Les plans factoriels complets et les réseaux de neurones artificiels 
ont pu remplacer le travail expérimental en simulant tous les phénomènes observés au 


laboratoire avec une erreur très faible. 
Le travail entamé dans cette thèse se décline en deux parties distinctes : 


La première partie est consacrée à l’étude de l’influence de quatre facteurs sur la 
conductivité thermique des nanofluides, Ces facteurs sont : la fraction volumique, la 
température, le diamètre et le type des nanoparticules. Un modèle mathématique est 
développé pour prédire la conductivité thermique des nanofluides d’Oxydes métalliques/Eau 
avec des nanoparticules d'Al203, TiO: et CuO, en utilisant un plan d'expérience (Design of 
Experiments DOE), ensuit ce modèle est testé et analysé par la méthode ANOVA pour 
déterminer les effets des différents facteurs et leurs interactions. Une méthode factorielle 
complète à deux niveaux est considérée pour chaque facteur, puis les résultats sont analysés. 
Cette modélisation est basée sur des données expérimentales recueillies dans la littérature, et 
les résultats prédits ont été comparés à ceux d'autres modèles ainsi qu'aux résultats 
expérimentaux. Les résultats obtenus montrent que la fraction volumique des particules a 
l'effet le plus significatif sur l'amélioration de la conductivité thermique des nanofluides par 
rapport aux autres facteurs. Une augmentation de la fraction volumique et de la température 
des particules possède un effet favorable sur l'amélioration de la conductivité thermique des 
nanofluides, tandis qu’une taille de particule plus grande tend à la réduire. Ainsi, la 
concentration optimale de nanoparticules est essentielle à l’amélioration des performances des 
systèmes thermiques. Enfin, les résultats obtenus confirment la capacité du modèle proposé 
dans notre étude à prédire la conductivité thermique des nanofluides avec une grande 
précision. 

La seconde partie est dévouée à la modélisation de la conductivité thermique et de la 
viscosité de nanofluide AlO3/Eau se basant sur la méthode de réseau de neuronal artificiel. 
Prenant en compte la température, la fraction volumique, la taille et la forme de 
nanoparticules comme paramètres indépendants (en tant qu'entrées), le rapport de 
conductivité thermique (k,/k;) et le rapport de viscosité (u,F/U) sont considérés quant à 


eux comme des paramètres dépendants (en tant que sorties). Une approche MLP avec une 
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seule couche cachée est choisie comme la meilleure structure pour l'estimation du rapport de 


conductivité thermique et de viscosité dynamique du nanofluide Alumine-Eau, via de 
l'algorithme Levenberg—Marquardt (LM). Les nombre optimal de neurones dans les couches 
cachées sont 15 et 20 neurones pour le rapport de conductivité thermique et de viscosité 
dynamique du nanofluide considéré respectivement. La discussion des résultats obtenue, 
l’évaluation de la corrélation générée qui est effectuée par une comparaison avec d’autres 
corrélations appliquées au même régime, ainsi que la détermination des erreurs absolues et 


réelles, ont été représentés dans ce chapitre. 


Les principaux résultats qu’on peut ressortir de cette thèse sont les suivants : 

- Les valeurs R° et RMSE de 0,99 et 0,0074 respectivement, Ceci indique ainsi que le 
modèle de TCR fournit une bonne corrélation et que les valeurs prédites et 
expérimentales s’adaptent parfaitement. 

- La fraction volumique des particules a l'effet le plus significatif sur l'amélioration de 
TCR par rapport aux autres facteurs. 

- Une augmentation de la fraction volumique et de la température des particules possède 
un effet favorable sur l'amélioration de la conductivité thermique des nanofluides. 

- Plus la taille des particules est grande, plus le TCR est bas. 

- Une bonne concordance entre le modèle de TCR obtenu et les données 
expérimentales. 

- La meilleure architecture trouvée pour le taux de TCR et DVR est de 15 et 20 
neurones dans la couche cachée respectivement. 

- Le modèle (k,F/k;) fournit une bonne corrélation entre les valeurs prédites et 
expérimentales avec des valeurs de R° égales à 0,998 et 0,998 pour l'apprentissage et 
la validation. De plus, les valeurs de MSE de 0,000018 pour l'apprentissage et de 
0,0000162 pour la validation. 

- Une excellente corrélation est observée entre la réponse prédite et observée pour 
(Unr/Hy). Ainsi, des valeurs plus élevées de R2 = 0,999 et de faible MSE = 0,000788 
pour l'apprentissage, R2? = 0,99 et SME = 0,000678 pour la validation. 

- Les valeurs AARD de 0.013 % pour (lhF/u) et de 0.0019 % pour (k,F/kF) montrent 


une excellente corrélation entre les réponses prédites et observées de (k,F/k) et de 


(Unr/Ur). 


- La comparaison avec les diverses corrélations empiriques confirme que le modèle 


proposé à des performances supérieures dans la prédiction du TCR et du DVR des 
nanofluides Alumine-Eau et peut être considéré comme un outil pratique pour la tâche 


considérée. 


Une première perspective pour des améliorations futures est de continuée l’analyse du 
transfert de chaleur en mettent en exergue des nouveaux modèles prédictifs de la conductivité 
thermique et de la viscosité qui tiennent en compte des effets du mouvement Brownien, de la 
couche inter-faciale entre nanoparticules et fluide de base ainsi que le phénomène 
d’agrégation qui pourraient être des mécanismes à l’origine de l’augmentation de conductivité 
thermique du nanofluide. Il sera intéressant d’étudier l’impact de ces modèles sur le transfert 
thermique, mais il faudra de comparaisons avec de données expérimentales pour trouver le 
meilleur. 

Une autre voie est de simuler les nanofluides comme un fluide multi-composant. 
En d’autres mots, au lieu de considérer un seul milieu continu ayant de nouvelles 
propriétés et dans lesquels les particules sont dissoutes, on considère le fluide de base et les 
particules solides comme des éléments séparent. 

Aüïnsi que, l’avantage de l’emploi des nanofluides en tant que fluides caloporteurs 
dépend fortement d’un compromis entre l’augmentation de la conductivité thermique et 
l’augmentation de la viscosité. Les futures recherches devront donc être dirigées vers la 
sélection de matériaux, de forme et de taille de nanoparticules qui augmenteraient la 
conductivité thermique du mélange avec une augmentation modérée de la viscosité. 
Par exemple, l’usage de particules non-sphériques conférant un caractère rhéofluidifiant au 
mélange est une voie à explorer. En outre, les futures recherches devront également prendre 
en compte l’encrassement éventuel des conduites par les nanofluides et d’en étudier les 


causes. 
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Annexes 


I. Facteurs influençant les propriétés thermophysiques des 


nanofluides 


Les paramètres principaux qui influent sur la conductivité thermique du nanofluide sont les 
suivants: 


L.1 Effet de la concentration volumique de particules 


L'effet de la concentration volumique des particules sur l’amélioration de la 
conductivité thermique du nanofluide est montré dans les figures I-1.a et I-1.b, (travaux de 
sept groupes de chercheurs pour les nanoparticules d'A1:0; suspendus dans l’eau). La 
dimension de la particule et la température du nanofluide varient, cependant la tendance 
générale montre que la conductivité thermique augmente avec la fraction volumique des 
particules et que les concentrations de particules d’Oxyde métallique donne une amélioration 


allant jusqu’à 30% [1]. 


Masuda, at al. (1693 13 nm, 31.85C) 
Masuda, et al. (1893, 13 nm, 46.85C) 
—— Masude, et al. (1893, 13 nm, 66.85C) 


" 
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Thermal conductivity ratio (k /k ) 
Thermal conductivity ratio {k /k } 


Figure L.1 : Effet de la concentration volumique sur la conductivité thermique du nanofluide 
(ALO;/Eau) pour différents chercheurs [1]. 


1.2 Effet du type de nanoparticules 


Les figures I-2.a et I-2.b regroupent les résultats des études expérimentales de 
plusieurs chercheurs. Elles illustrent l’effet du type des nanoparticules sur la conductivité 
thermique pour deux types de nanofluides (un oxyde métallique et le second est un métal), 
tous deux sont dilués dans l’eau. Avec tous les paramètres (température et forme de la 
particule) sont supposés constants. Les résultats montrent que les particules métalliques 
produisent la même amélioration que les oxydes métalliques mais pour de faible 
concentrations volumiques. On observe aussi que le rapport de la conductivité thermique du 


métal augment plus rapidement que celui des oxydes métalliques [1]. 
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Figure L.2 : Effet des types de particules sur la conductivité thermique des nanofluides [1]. 


L.3 Effet de la température 

En général, la conductivité thermique du nanofluide est sensible à la température par 
rapport au fluide de base. Par conséquent, l’amélioration de la conductivité thermique du 
nanofluide est également sensible à la température et augmente en fonction d’elle. Les 
résultats obtenus par Das et al. [2] pour des nanoparticules d’Al,0:/Eau et CuO/Eau avec un 
intervalle de température de 28°C à 51°C, indiquent clairement que la conductivité thermique 
augmente lorsque la température augmente. Mais elle est plus élevée pour l’oxyde de cuivre 
CuO par rapport à l’oxyde d’aluminium Al:03. Où le taux de la conductivité thermique pour 
une concentration volumique p—4% et une température T=51°C, atteigne une valeur de 1,24 


pour l'ALO; et 1,36 pour le CuO. 
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—<— Das et al. [9] (38 nm, 51C) —— Das et al, [9] (29 nm, 51C) 
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Figure L.3 : Effet de la température sur la conductivité thermique des nanofluides à 
base d’eau : (a) pour ALO;, (b) pour CuO [1]. 


L.4 Effet de la taille des particules 

La conductivité thermique d’un nanofluide augmente avec la diminution de la taille 
des particules. Ce phénomène est dû aux effets du mouvement Brownien des particules, Teng 
et al. [3]; Le coefficient de diffusion Brownien augmente lorsqu'on réduit la taille des 
particules ; plus les particules sont petites plus la probabilité de collision entre elle augmente 
et meilleure est la conductivité thermique. 
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Figure L.4 : Influence du diamètre des particules sur la conductivité thermique. 


Le mécanisme de la couche inter-faciale solide-liquide renforce la théorie stipulant que 


la diminution de la taille des particules augmente la conductivité thermique des nanofluides. 


Yu et al. [1] prétendent le contraire ; plus les particules sont grandes plus la 
conductivité thermique augmente. Ozerinç [4] ainsi que Bazan [5] ont également trouvé de 
grandes disparités concernant l’évolution de la conductivité en fonction de la taille des 
particules. Pour un nanofluide Al-Ethylène-Glycol, Xie et al. [6] ont révélé une augmentation 
de la conductivité pour des particules de grandes tailles (environ 60 nm) et une diminution 
pour les petites particules. Ozerinç explique cela par l’effet du groupement (clustering) ; plus 


les particules sont petites plus cet effet est important. 


L5 Effet de la géométrie des particules 


SiC in Water (Particle Shape Effect) 
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Figure LS : Influence de la géométrie des particules sur la conductivité thermique [1]. 


Aucune amélioration du transfert de chaleur n’est observée dans le cas de particules en 
forme de disque. Le graphique de la figure LS montre que des particules à géométrie 
cylindrique (longueur indiquée) donnent une meilleure conductivité que des particules 
sphériques. Timofeeva et al. [7] ont étudié l’influence de la géométrie des particules sur les 
propriétés physiques (conductivité, viscosité et pH du fluide). Les nanofluides les plus utilisés 
ont des particules sphériques car leur production est moins coûteuse que celle des particules 


cylindriques, par exemple. 


L.6 Effet de pH du fluide 

Les études sur l’influence de l’acidité du fluide porteur d’un nanofluide sur la 
conductivité thermique de ce dernier sont assez limitées. Üzerinç, Yu et al. ainsi que Xie et al. 
[1, 4, 6] prétendent que la valeur du pH influence la conductivité thermique ; ils ont constaté 
un optimum pour un pH compris entre 2 et 3. 

Timofeeva et al. [7] ont déterminé l’influence du pH sur la conductivité thermique 
d’un nanofluide à base de particules cylindriques associées à de l’éthylène glycol ou de l’eau 
et avec des pH de 2,54 ; 3,33 et 4,10. Ils ont constaté que le pH influence peu la conductivité 


thermique et que la viscosité des fluides augmente en diminuant le pH. 


CuoO in Water (Acidity Effect) 
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Figure L.6 : Influence du pH sur la conductivité thermique d’un nanofluide CuO/Eau 
(APS=60nm) [1]. 


IL. Caractéristiques des nanofluides adoptées par différents chercheurs 
Les tableaux suivants présentent un résumé des résultats obtenu par différents groupes 
de recherche de taux de la conductivité thermique pour différentes variables principales et 
ceci pour des oxydes métalliques, des métaux et carbone silicone. 
a. Les oxydes métalliques 


Tableau IL.1 : Effet de la taille des particules sur la conductivité thermique [1]. 


ALO;/ Eau 
CuO/ Eau 
ALO;/ EG 
CuO/ EG 
ALO;/ Eau 
CuO/ Eau 


Wang et col (1999) AlLO;/ EG 
CuO/ EG 


ALO;/ Eau 


ALO;/ Eau 5.00 

Xie et col (2002c) ALO:/EG 5.00 

Wen et Ding (2004) Al:O;/ Eau 0.19-1.59 42 SE 
TiO;/ Eau 0.50-5.00 15 1.05-1.30 

MATSET ÉCOLOS, | ne 0.50-5.00 10 x 40 1.08-1.33 

Wen et Ding (2005) A1LO;/ Eau 0.31—0.72 [ 1.02-1.06 


CuO/ Eau 1.00 1.05 
Hwang et col (2006) SiO/ Eau 1.00 1.03 
CuO/ EG 1.00 1.09 


Kang et col (2006) | et col (2006) SO | D] Eau 1. EC 00 15-20 1.02-1.05 

En e Eau — 3)  _ 030.30 .04-1.12 

CES Cons) CuO/ Eau (pH=6)  |0.03-0.30 1.02-1.07 
29- T-02-1. 


_ ____ _ en 


Tableau II.2 : Effet de la température et de la taille des particules sur la conductivité 


Auteur (Année) 


Masuda et col (1993) 


Li et Peterson (2006) 


thermique [1]. 


Nanofluide 


AlLO;/ Eau (31.85°C) 
A1O;/ Eau (46.85°C) 
A1LO;/ Eau (66.85°C) 
Si0;/ Eau (31.85°C) 
SiO7/ Eau (46.85°C) 
SiO}/ Eau (66.85°C) 
TiO Eau (31.85°C) 
TiO7 Eau (46.85°C) 
TiO,/ Eau (86.85°C) 
ALO;/ Eau (21°C) 
A1LO;/ Eau (36°C) 
ALO3/ Eau (51°C) 
CuO/ Eau (21°C) 
CuO/ Eau (36°C) 
CuO/ Eau (51°C) 


ALOy Eau (27.5°0) 
ALO;/ Eau (32.5°C) 
ALO:/ Eau (34.7°C) 
CuO/ Eau (28.9°C) 
CuO/ Eau (31.3°C) 
CuO/ Eau (33.4°C) 


b. Les métaux et carbone silicone 


2.00-10.00 
2.00-10.00 
2.00-10.00 
2.00-6.00 
2.00-6.00 
2.00-6.00 


1.109-1.324 
1.100-1.296 
1.092-1.262 
1.010-1.011 
1.009-1.010 
1.005-1.007 
1.080-1.105 
1.084-1.080 
1.075-1.099 


Tableau IL.3 : Effet de la température et de la taille des particules sur la conductivité 


Auteur (Année) 


Wen and Ding (2004) 


thermique [1]. 


Nanofluide 


Ag/ Eau (30°C) 
Ag/ Eau (60°C) 
Au/ Eau (30°C) 
Au/ Eau (60°C) 
Au/ Eau (30°C) 
Au/ Eau (60°C) 


MWCNT/eau (20°C) 
MWCNT/eau (45°C) 


MWCNT/eau (20°C) 
MWCNT/eau (25°C) 
MWCNT/eau (30°C) 


0.04-0.84 20-60 (D) 1.04-1.24 
0.04-0.84 20-60 (D) 1.05-1.31 


Tableau II.4 : Effet de la taille des particules sur la conductivité thermique [1]. 


| : Taille de la ne : 
Auteur (Année) Nanofluide La concentration : L’amélioration 
particule 


Xuan et Li (2000) 250750 122-4175 


Eastman et Col (2001) .10-0. <10 1.016-1.100 
<10 1.031-1.140 
<10 1.002-1.140 


Xie, et Col (2002a) S1C/ Eau 0.78-4.18 26 
SiC/ Eau 1.00-4.00 600 
SiC/EG 0.89-3.50 
SiC/EG 1.00-4.00 


Xie et Col (2003) MWOCNT/ eau 0.40-1.00 15 nmx30pum 1.03-1.07 
so Hi Sn 
SE 
Assael et Col (2005) DWCNT/eau : : 

DWCNT/eau 

MWCNT/eau | 130 nmx10Oum 


MWCNT/eau : 130 nmx10Oum 


Hong et Col (2005) Fe/EG 0.20-0.55  : 13-1.18 
Liu et Col (2005) MWCNT/EG 0.20-1.00  _ (D) 1.02-1.12 
MWCNT/H M 1.00-2.00 20-50 (D) 1.09-1.30 
Hong et Col (2006) | et Col (2006) F/EG | EG 010-055 | 10-0.55 EC: 105-118 | 1.18 
a et Col (2006) TE — Eau a — 00 Co LRU 07 
MWCNT/ H Min 0.50 1.09 
Kang et Col (2006) Ag/ Eau 0.10-—0.39 8-15 1.03-1.11 
Diamant/EG 0.13-1.33 30-50 1.03-1.75 


EG : Ethylène glycol, H M : Huile de moteur, H Min : Huile minérale. 


_ _ _ _ _ 


III. Exemples des types de règles d'apprentissage 


Des exemples d'utilisations de différentes règles d'apprentissage, d’architectures et 
d’algorithme qui sont relié à des types supervisé et non supervisé sont illustrés dans les 


tableaux IIL.1 et IIL.2 


Tableau IIL.1 : Exemple de différents types d'apprentissages en mode supervisé [8]. 


Paradigme | Architecture Règle Algorithme Taches 
d'apprentissage 
Perceptron | Correction d'erreur | Perception Rétro- Classification 
Simple ou Propagation, Approximation de 
Multi- Adaline, Madaline fonctions, Prédiction, 
Couches Contrôle 
Récurrente Boltzmann Apprentissage de Classification 
Bolzmann 
Supervisé Multi- Hebb Analyse de Analyse de données 
Couches non discriminants Classification 
bouclés linéaires 
à Par compétition LVQ Catégorisation au sein 
compétition d'une classe Compression 
de données 
ART ARTMap Classification 
Catégorisation au sein 
d'une classe 


Tableau IIL.2 : Exemple de différents types d'apprentissages en mode non supervisé [8]. 


Paradigme | Architecture Règle Algorithme Taches 
d'apprentissage 
Multi- Correction d'erreur Projection de Analyse de données 
Couches non Sammon 
bouclés 
Non bouclés Hebb Analyse en Analyse de données 
ou à composantes 
Non compétition 
Supervisé à Par compétition VQ Catégorisation, 
compétition Compression de données 
Carte de SOM Catégorisation, 
Kohonen Analyse de données 
ART ART-1, ART-2 Catégorisation 
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1. Propriétés et prétraitement des données utilisées dans le chapitre IV 


Les propriétés thermophysiques utilisées dans le chapitre IV pour les nanoparticules 


d'ALO; et le fluide de base (Eau) ainsi que le nanofluide Al:03/Eau sont présentées dans les 


tableaux 1 et 2. 


Tableau 1: Propriétés de l'eau pure, de particule et de nanofluide Al1:0;/Eau [1]. 


Re pire Propriétés de Propriétés de 

particule Al:0; Nanofluid 

Densité, p (kg/m°) 987,17 3600 989,78283 
Chaleur spécifique, cp (J/Kkg.K) 4182 765 4169,571819 
Conductivité thermique, k (W/m.K) 0,645 36 0,646836333 
Viscosité dynamique, u (Ns/m°) 0,000528 0 0,122426401 
Coefficient de dilatation thermique, B (1/K) 0,0004712 00000058 0,000469507 
Conductivité thermique, k (W/m.K), Statique / / 0,646836333 
Conductivité thermique, k (W/m.K), Brwonain / , 0,122426401 


Tableau 2 : Variables de calcul pour Al1:0;/ Eau [1]. 


Variables calculassions valeurs 
0 0,001 
dp 2,00E-08 
Fraction de volume de liquide, f pour ALO; 99,86618619 
Température (K ) 325 
CT, ) 0,000752209 
constante de Boltzmann J/K 1,38E-23 
Poids moléculaire du fluide de base, M 0,01801528 
(kg/mol) 
Masse volumique de l'eau at 293k 9,87E+02 
Numéro Avogadro, N (mol”(-1)) 6,02E+23 
Diamètre moléculaire de l'eau, df (m) 3,86815E-10 
Diamètre hydraulique 0,0203 


On utilise un tableau nommé matrice d'expériences pour récapituler l'ensemble des 
essais. Les données expérimentales de conductivités thermiques effectives et de viscosités du 


nanofluide Al,03/Eau obtenues à partir de la littérature, sont présentées dans les tableaux 3 et 


4 respectivement [1]. 


Tableau 3: Matrice de k,/k; 


X1 X2 X3 X4 Yi: Kaç/ky 
1 300 20 sphérique 1,11 
1 300 40 sphérique 1,09 
1 300 60 sphérique 1,077 
1 300 80 sphérique 1,075 
1 300 100 sphérique 1,069 
2 300 20 sphérique 1,14 
2 300 40 sphérique 1,11 
2 300 60 sphérique 1,109 
2 300 80 sphérique 1,1 
2 300 100 sphérique 1,096 
3 300 20 sphérique 1,169 
3 300 40 sphérique 1,14 
3 300 60 sphérique 1,136 
3 300 80 sphérique 1,13 
3 300 100 sphérique 1,125 
4 300 40 sphérique 1,174 
4 300 60 sphérique 1,167 
4 300 80 sphérique 1,159 
4 300 100 sphérique 1,155 
5 300 20 sphérique 1,22 
5 300 40 sphérique 1,2 
5 300 60 sphérique 1,197 
5 300 100 sphérique 1,186 
1 300 20 lame 1,027 
2 300 20 lame 1,059 
3 300 20 lame 1,082 
4 300 20 lame 1,11 
5 300 20 lame 1,138 
1 300 20 brique 1,03 
2 300 20 brique 1,06 
3 300 20 brique 1,1 
4 300 20 brique 1,13 
5 300 20 brique 1,17 
1 300 20 cylindrique 1,04 
2 300 20 cylindrique 1,08 
3 300 20 cylindrique 1,12 
4 300 20 cylindrique 1,159 


5 300 20 cylindrique 1,2 

1 300 20 plaquette 1,27 
4 300 20 plaquette 1,11 
1 300 20 sphérique 1,11 
2 300 20 sphérique 1,14 
3 300 20 sphérique 1,17 
4 300 20 sphérique 1,2 
5 300 20 sphérique 1,22 
1 300 20 sphérique 1,13 
2 300 20 sphérique 1,14 
3 300 20 sphérique 1,17 
4 300 20 sphérique 1,19 
5 300 20 sphérique 1,23 
1 305 20 sphérique 1,131 
2 305 20 sphérique 1,15 
3 305 20 sphérique 1,18 
4 305 20 sphérique 1,22 
5 305 20 sphérique 1,24 
1 310 20 sphérique 1,15 
2 310 20 sphérique 1,17 
3 310 20 sphérique 1,2 
4 310 20 sphérique 1,23 
5 310 20 sphérique 1,26 
1 315 20 sphérique 1,165 
2 315 20 sphérique 1,192 
3 315 20 sphérique 1,21 
4 315 20 sphérique 1,25 
5 315 20 sphérique 1,28 
1 320 20 sphérique 1,18 
2 320 20 sphérique 1,2 
3 320 20 sphérique 1,23 
4 320 20 sphérique 1,26 
5 320 20 sphérique 1,3 

1 325 20 sphérique 1,19 
2 325 20 sphérique 1,22 
3 325 20 sphérique 1,25 
4 325 20 sphérique 1,28 
5 325 20 sphérique 1,31 
1 300 20 sphérique 1,117 
1 300 40 sphérique 1,09 
1 300 60 sphérique 1,08 
1 300 80 sphérique 1,07 
1 300 100 sphérique 1,07 
2 300 20 sphérique 1,14 
2 300 40 sphérique 1,117 
2 300 60 sphérique 1,1 


2 300 80 sphérique 1,1 
2 300 100 sphérique 1,097 
3 300 20 sphérique 1,16 
3 300 40 sphérique 1,14 
3 300 60 sphérique 1,138 
3 300 80 sphérique 1,34 
3 300 100 sphérique 1,126 
4 300 20 sphérique 1,196 
4 300 40 sphérique 1,17 
4 300 60 sphérique 1,167 
4 300 80 sphérique 1,163 
4 300 100 sphérique 1,15 
5 300 20 sphérique 1,22 
5 300 40 sphérique 1,2 
5 300 60 sphérique 1,2 
5 300 80 sphérique 1,19 
5 300 100 sphérique 1,186 
Tableau 4: Matrice de u,/u 
X1 X2 X3 X4 Y2: Unf/Uf 
1 300 20 sphérique 1,1 
1 300 40 sphérique 1,092 
1 300 60 sphérique 1,084 
1 300 80 sphérique 1,075 
1 300 100 sphérique 1,058 
2 300 20 sphérique 1,23 
2 300 40 sphérique 1,18 
2 300 60 sphérique 1,16 
2 300 80 sphérique 1,15 
2 300 100 sphérique 1,13 
3 300 20 sphérique 1,39 
3 300 40 sphérique 1,3 
3 300 60 sphérique 1,26 
3 300 80 sphérique 1,24 
3 300 100 sphérique 1,22 
4 300 20 sphérique 1,64 
4 300 40 sphérique 1,46 
4 300 60 sphérique 1,38 
4 300 80 sphérique 1,33 
4 300 100 sphérique 1,32 
5 300 20 sphérique 1,96 
5 300 40 sphérique 1,66 


EE 


5 300 60 sphérique 1,54 
5 300 80 sphérique 1,47 
5 300 100 sphérique 1,42 
1 300 20 sphérique 1,1 
2 300 20 sphérique 1,23 
3 300 20 sphérique 1,4 
4 300 20 sphérique 1,64 
5 300 20 sphérique 1,95 
1 300 40 sphérique 1,07 
2 300 40 sphérique 1,18 
3 300 40 sphérique 1,3 
4 300 40 sphérique 1,46 
5 300 40 sphérique 1,66 
1 300 60 sphérique 1,07 
2 300 60 sphérique 1,16 
3 300 60 sphérique 1,26 
4 300 60 sphérique 1,39 
5 300 60 sphérique 1,53 
1 300 80 sphérique 1,07 
2 300 80 sphérique 1,14 
3 300 80 sphérique 1,24 
4 300 80 sphérique 1,34 
5 300 80 sphérique 1,47 
1 300 100 sphérique 1,057 
2 300 100 sphérique 1,13 
3 300 100 sphérique 1,21 
4 300 100 sphérique 1,31 
5 300 100 sphérique 1,42 
1 300 20 lame 1,18 
2 300 20 lame 1,53 
3 300 20 lame 1,56 
4 300 20 lame 1,77 
5 300 20 lame 2,04 
1 300 20 brique 1,16 
2 300 20 brique 1,29 
3 300 20 brique 1,45 
4 300 20 brique 1,81 
5 300 20 brique 2,31 
1 300 20 cylindrique 1,22 
2 300 20 cylindrique 1,66 
3 300 20 cylindrique 2,22 
4 300 20 cylindrique 2,97 


EE 


5 300 20 cylindrique 3,93 
1 300 20 plaquette 1,43 
2 300 20 plaquette 1,97 
3 300 20 plaquette 2,66 
4 300 20 plaquette 3,45 
5 300 20 plaquette 4,37 
1 300 20 sphérique 1,08 
2 300 20 sphérique 1,25 
3 300 20 sphérique 1,41 
4 300 20 sphérique 1,64 
5 300 20 sphérique 1,95 
1 300 20 sphérique 1,11 
2 300 20 sphérique 1,23 
3 300 20 sphérique 1,41 
4 300 20 sphérique 1,64 
5 300 20 sphérique 1,95 
1 305 20 sphérique 0,84 
2 305 20 sphérique 94 

3 305 20 sphérique 1,07 
4 305 20 sphérique 1,25 
5 305 20 sphérique 1,5 

1 310 20 sphérique 0,75 
2 310 20 sphérique 0,83 
3 310 20 sphérique 0,97 
4 310 20 sphérique 1,12 
5 310 20 sphérique 1,35 
1 315 20 sphérique 0,68 
2 315 20 sphérique 0,77 
3 315 20 sphérique 0,89 
4 315 20 sphérique 1,02 
5 315 20 sphérique 1,21 
1 320 20 sphérique 0,64 
2 320 20 sphérique 0,69 
3 320 20 sphérique 0,8 

4 320 20 sphérique 0,93 
5 320 20 sphérique 1,12 
1 325 20 sphérique 0,56 
2 325 20 sphérique 0,65 
3 325 20 sphérique 0,74 
4 325 20 sphérique 0,86 
5 325 20 sphérique 1,03 


Chaque forme est codifiée sous une forme numérique selon le tableau ci-dessous : 


Tableau 5 : Codification de facteur x4 


Forme de Nh | Codification 
sphérique 0 
lame 0,25 
brique 0,5 
cylindrique 0,75 
plaquette 1 
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Abstract 

In this paper, a full factorial design analysis is proposed 
for predicting nanofluid thermal conductivity ratio 
(TCR) as well as determining the effects of critical 
factors and their interactions. A statistical design of 
experiment approach with three variables (volume 
fraction, temperature, and nanoparticle diameter) at 
two levels is carried out. Three types of oxide-water 
nanofluids (ALO;-water, CuO-water, and TiO,-water) 
are used to evaluate the effectiveness of the proposed 
mathematical model. The significance and adequacy of 
the regression model were evaluated by the analysis of 
variance. The predicted model has a root mean square 
error equals to 0.0074, R?=0.99, and P<.0013, thus 
showing good results compared to a set of experimental 
data as well as other mathematical model results. The 
results illustrate that the TCR of metallic oxide 
nanofluids increases with temperature and nanoparti- 
cles volume fraction but decreases when nanoparticle 
size intensifies. Furthermore, it is found that the 
nanoparticles volume fraction has a great impact on 
the nanofluids thermophysical properties. Finally, the 
obtained results confirm that the proposed model is 
considerably accurate and capable of predicting na- 
nofluids thermal conductivity and that it can be used 
with ease as an alternative to many other models. 
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1 | INTRODUCTION 


Over the last few decades, many researchers have been trying to find the appropriate technique 
to improve the heat transfer rate and to increase the thermal conductivity of fluids. One of the 
recent technologies of improving heat transfer is suspending nano-sized solid particles of 
metallic or oxide-metallic with high intrinsic thermal conductivity in a base fluid. This mixture 
is called a nanofluid.!* 

Nanofluids thermophysical properties can be controlled by chemical and physical 
specifications of the base fluid and nanoparticles. Nanofluids thermal conductivity is one of 
the factors that play an important role in heat transfer performance. Thus, using a suitable type 
of nanofluids is highly important in improving the heat transfer rate.®” However, research 
results illustrate wide discrepancies and inconsistencies in nanofluids thermal conductivity due 
to different parameters, such as the difference in nanofluids production and stabilization 
techniques, nanofluids characterization as well as thermal conductivity measuring methods.f 

A strong relationship exists between nanofluids thermal conductivity and various 
parameters, namely: nanoparticle volume fraction, particle material, size and shape, base fluid 
material and temperature, specific heat capacity, type of surfactant, and density.” !? 

The overall aim of using nanofluids is to have optimum thermophysical properties at 
minimum particle volume fractions. Moreover, it has been reported that the thermal 
conductivity of nanofluids increases with the increase of nanoparticle concentration.”""#?° It 
has also been revealed that the thermal conductivity enhancement ratio declines with that of 
the base fluid, while it is directly proportional to the nanoparticle thermal conductivity.?!?? 
Nevertheless, for a certain nanoparticles concentration, it is shown that for a lower thermal 
conductivity the heat transfer coefficient will be enhanced. This fact agrees with conventional 
models such as the Maxwell model. On the other hand, it has been shown by Alrashed et al 
that nanofluids thermal conductivity enhancement is independent of the base fluid. 

Another important factor that affects the nanofluids thermal conductivity is the 
nanoparticles’ size, since it has been confirmed that the highest thermal conductivity is 
achieved by a uniform distribution of small particles! It is also assumed that an increase in the 
nanoparticles’ specific surface area could lead to an increase in the effective nanofluids thermal 
conductivity because heat transfer between particles and the base fluid occur at particle/fluid 
interface. Moreover, the Brownian motion of nanoparticles could be intensified by reducing the 
particle's size. Hence, the latter is claimed to be one of the most important factors in enhancing 
nanofluids’ heat transfer. /"# The shape of nanoparticles is another parameter affecting the 
enhancement degree of nanofluids thermal conductivity. The enhancement ratio of the 
elongated particles is greater than that of the spherical ones.” In general, a decrease in 
nanoparticle sphericity leads to an increase in nanofluids thermal conductivity. 

Generally, nanofluids thermal conductivity is more sensitive to temperature and the increase 
in temperature ensued an increase of the thermal conductivity enhancement.**®?% Many 
studies showed the existence of an optimum temperature at which nanofluids thermal 
conductivity is maximum.” Also the pH has an effect on the nanofluids thermal conductivity as 
its enhancement ratio increases with acidity.!” Furthermore, the thermal conductivity 
enhancement ratio can be affected by surfactants and additives used to prevent particle 
agglomeration.!*?7 A summary of experimental studies on nanofluids thermal conductivity can 
be found in Rashmi et al.?* 

From the above review, it can be noticed that all the different parameters and mechanisms 
affecting the thermal conductivity of nanofluids cannot be accounted by the traditional models 
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TABLE 1 Summary of representative studies on thermal conductivity of water-based nanofluids adopted in 
this study 


Particle volume Particle 
Nanoparticle fraction (%) size, nm Temperature, K Researchers 
CuO 0-18 29, 28.6 Room temperature Aminian** 
ALO; 1-4 20, 25, 28 
CuO 0-5 28.5 Differenttemperature Kazemi-Beydokhti 
et al*° 
ALO; 0-5 20, 25 
TiO; 0.2-2 21 288-308 Duangthongsuk et al* 
TiO; 0-8 21 288-338 Yiamsawasd et al”? 
ALO; 1-4 30 +10 274-313 Longo et al 
TiO; 1-6 30-50 + 10 


of particle suspensions.” Thus, several models have been proposed, such as the work of 
Alrashed et al”?? using an enhanced artificial neural network (EANN) and an adaptive neuro- 
fuzzy inference system approach and the results showed a suitable accuracy with experimental 
data. Hosseini et al” presented a nanofluid thermal conductivity model as a function of 
nanoparticle thermal conductivity, taking into consideration the working fluid, temperature, 
nanoparticle volume fraction, and size as well as the properties of the interfacial shell. Their 
results demonstrated good compatibility with experimental data and other models. Shahsavar 
et al*° explored several experiments with the artificial neural network method to measure and 
model nanofluid thermal conductivity and viscosity. The results of their model can be employed 
to provide the calculation of properties. Bagherzadeh et al°! developed an EANN approach of 
nanofluid thermal conductivity. Their proposed model gave accurate results with less 
computational time. À complete review of these models can be found in references.*?** 

Due to wide discrepancies and inconsistencies results shown in literature, additional 
investigations to elucidate nanofluids thermal conductivity effects must be carried out. Thus, the 
present research work aims to find a new correlation that takes into account all the parameters 
influencing the thermal conductivity by using expanded and extensive experimental data. 

Therefore, in this study, a mathematical model is developed to predict the thermal 
conductivity of oxide-water nanofluids. Four major influencing factors on the nanofluid 
thermal conductivity have been selected, and design of experiment (DOE) with the analysis of 
variance (ANOVA) method are used to determine the effects of these factors and their 
interactions. A full factorial method is considered at two levels for each factor.? and the results 


TABLE 2 Thermophysical properties of water and selected nanoparticles 


k, kg/m:s P, kg/m° Cp, J/Kkg-K K, W/m-K 


Base fluid Water (H20) 0.000912 997.047 4184 0.6063 
Nanoparticles Alumina (ALO:;) . 3970 765 40 
Titanium dioxide (TiO;) a. 4250 686.2 8.9538 


Copper oxide (CuO) 6500 535.6 20 
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TABLE 3 Factors and levels used in the factorial design 


Factors Low level (—1) High level (+1) 
(X1) Particle volume fraction (%) 1 4 

(2) Temperature, K 295 298 

(33) Particle size, nm 21 29 


are analyzed and plotted by DOESs software. The present modeling is based on experimental 
data collected from the literature and the predicted results were compared to those of other 
models as well as with experimental results. 

The novelty of the present work is that it only uses significant factors, thus needing less 
experimental data and computational time. Furthermore, operating conditions of the present 
research work are at low temperature and medium volume fraction, while the majority of the 
empirical and statistical models in use, are applicable at higher temperature ranges. Hence, the 
purpose of this study is to formulate a new model for predicting the thermal conductivity of 
three types of nanofluids at low temperatures and medium volume fractions. 


2 | EXPERIMENTAL DATA 

The experimental data used in this study are extracted from literature***; the details are 
presented in Table 1. Metal oxide nanoparticles are used in the present work because they are 
the most employed in heat transfer applications due to their low cost. They are produced by 
conventional methods and are more stable than metallic nanoparticles.*? Three types of oxide 
metal nanoparticles are studied: Aluminum oxides (Alumina), ALO;“*; Titanium oxides 
(Titania), TiO:; and Copper oxides, CuO.*!*? The corresponding thermophysical properties of 
the base fluid (water)** and solid phases are given in Table 2. 


3 | FULL FACTORIAL DESIGNS 


A full factorial design is a powerful technique, quick and cost-effective, employed to achieve the 
best optimization. This design allows us to determine the effects on the response of each of the 
factor levels. Furthermore, it takes into consideration all the combinations between factors’ 
levels as well as the effects of interactions between different factors. This technique is used to 
reduce the number of experiments, time and overall process cost with a better response. The 
importance of this method is shown by different studies using experimental designs.*** 


4 | MATHEMATICAL MODEL 


In this study, an experimental design method is adopted to optimize nanofluids thermal 
conductivity. The DOE technique is performed to identify the key variables that affect the thermal 
conductivity ratio (TCR; K,#/Ks) and to determine at what levels these factors must be kept to 
reduce measurement errors. The full factorial design method is used to study the effects caused by 
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independent factors as well as interactions between those self-governing factors.“*"7 In the 
present work, a full factorial design at two levels for three independent factors is considered: 
volume fraction (X,), temperature (X,), and nanoparticle diameter (X3), with two factorial levels 
coded as —1, +1, for low and high levels, respectively. The types of the nanoparticles (X,) are 
coded as L; for (ALO:), L; for (TiO:), and L, for (CuO). A full factorial design 2° (ie, three factors 
are varied over two levels) for three types of nanoparticles is utilized. Therefore, a total of 24 
experimental runs are suggested for analyzing the interaction of each level (Table 4). The ratio 
(Knt/K5) is considered as the response factor (dependent factor). Results are analyzed and plotted 
with 95% confidence intervals using experimental design software. 

The polynomial equation for the experimental design with three factors is given by the 
following Equation (1): 


Y = do + >, dix; + >» (EE D aux? + Aij.zXiXj…. Xy + € (1) 


TABLE 4 Design matrix and the experimental responses 


Independent factors 


Dependent factors 
X X X, X 
Nanoparticles Experiment ? T dp Np Y= Kar/Kf 
ALO: 01 ue “4j L1 1.017 
02 1 —1 —1 L1 1.09 
03 —1 1 —1 L1 1.02895 
04 1 1 —1 L1 1.100320 
05 —1 —1 1 L1 1.0257 
06 1 —1 1 L1 1.093 
07 —1 1 1 L1 1.025400 
08 1 1 1 L1 1.105400 
TiO;, 09 —1 —1 —1 L2 1.052300 
10 1 —1 —1 L2 1.096330 
11 —1 1 —1 L2 1.047505 
12 1 1 —1 L2 1.130000 
13 —1 —1 1 L2 1.047070 
14 1 —1 1 L2 1.111920 
15 —1 1 1 L2 1.018215 
16 1 1 1 L2 1.115060 
CuO 17 —1 —1 —1 L3 1.057 
18 1 —1 —1 L3 1.154 
19 —1 1 —1 L3 1.136 
20 1 1 —1 L3 1.16 
21 —1 —1 1 L3 1.05835 
22 1 —1 1 L3 1.088200 
23 —1 1 1 L3 1.0609 
24 1 1 1 L3 1.110300 


Independent factors: is the factor that is changed or controlled in the experiment, in this study three independent factors is 
considered: volume fraction (@), temperature (T) and particle size (dp) which are coded X1, X2 and X3 respectively. Np: is the 
types of the nanoparticles (X4) are coded as L1 for (A1203), L2 for (TiO2) and L3 for (CuO). Dependent factors: is the variable 
that is being measured or tested in an experiment, in this study the thermal conductivity ratio (Knf/Kf) is considered as the 
response factor (dependent factor) which coded by Y. 


GRINE ap BENHAMZA HEAT TRANSFER WI LEY 2935 


Asian Research 


TABLE 5 Constants values 


Constants 


Nanoparticles C1 X10 * C:X10 * C;X10 * C,X10 * CsX10 * C,X10 * C,x107* 


ALO; 39.52 —18.36 103.99 2.13 5.18 —23.92 —196.5 
TiO; 35.21 57,597 45.16 76.37 45.48 —49.92 —30.71 
CuO —74.74 75.73 —149.1 —78.51 —50.66 73.85 227.22 


Where y is the response of a domain determined experimentally with a definite accuracy, x is 
the ith variable or a reduced centered coordinate; a,, a; a;; a;; are model coefficients determined 
experimentally and e is the total standard deviation including the experimental and model 
selection errors. Tables 3 and 4 show the selected four factors and different factor level settings 
as well as the obtained response matrix at each run. A convenient selection of the experimental 
domain allows to cover the required optimal zone. 


5 | RESULTS AND DISCUSSION 


The mathematical model generated by the experimental design method, obtained on the basis 
of experimental data has the following polynomial form: 


Y = 1.083716666 + (0.032505833 + C1)X1 + (0.0061325 + C)X 
+ (-0.008745416 + C;)X3 + (-0.005212916)X2X; (2) 


+ (C)X2X + (C)XX3 + (0.00418375 + Ci)XX2X3 + Cz. 
— 0.02 -296. dy - 0.000000025 
Where X — Te Ë. 2 — F — 3, 3 — Ro à and Cï, C, GC Ca Cs, Cé, and C; are 
constants depending on the type of nanoparticle determined by its physical properties. The 
values of these constants for each nanoparticle are given in Table 5. 
The partial models of each of the three types of nanofluids are as follows: 


AbO3: Y = 1.0640665 + 0.036458749X; + 0.00429625X; + 0.00165375X3-0.005212916X; X; 
+ 0.00021375X,X> + 0.000518333X,X3 + 0.00179125X, X X3. (3) 


TiO;: Y = 1.080645 + 0.036027499X, + 0.000395X,-0.0039375X;-0.005212916X;X; 
+ 0.0076375X,X; + 0.004548333X, X;-0.00080875X X X3. (4) 


CuO: Y = 1.106438749 + 0.02503125X, + 0.01370625X2-0.023656249X;-0.005212916X;X; 
-0.00785125X,X2-0.005066666XX; + 0.0115675X XX. (5) 


5.1 | Testing of the proposed model 


The proposed model is tested by comparing the TCR (K,;#/K5) of the predicted values with the 
actual values, taking into account the standard error (Figure 1). 
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1 1.05 11 1.15 
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FIGURE 1  Actual and predicted values, for TCR models. (P = .0013; R? = 0.99: RMSE = 0.0074). RMSE, 
root mean square error [Color figure can be viewed at wileyonlinelibrary.com|] 


It is clearly shown that all the values of TCR are within the central line, with a value of root 
mean square error of (0.0074). Furthermore, the TCR model provides a good correlation 
between predicted and experimental values with R° equals to 0.99. 

Table 6 gives the estimated parameters and their standard errors and the corresponding f test 
to each value of the independent variable. The obtained results show that all of the standard 
errors are very small (<0.003), which indicates that the results are statistically significant. The 
ANOVA for the full factorial design of the polynomial model is presented in Table 7. It can be 
noticed that the proposed model is significant with a P value of .0013. Besides, the confidence 
level for linear terms X, (P < .0001), X3 (P < .0045), and X4 (P < .0008) are most significant on 
TCR, while X, (P < .0157) is significant. Most of the square terms and their interactions are also 
highly significant. Likewise, it can be seen from Table 7 that the particle volume fraction is the 
most significant term since it has a very high F value. Hence, the particle volume fraction has 
the most significant effect factor on the thermal conductivity enhancement ratio. The particle 
diameter is also an important factor (with high F value). Although the temperature term has a 
lower F value compared to other main parameters, an F value of 16.26 suggests that the 
temperature has also a significant effect on the response. In addition, R° and adjusted R? values 
of 0.99 and 0.97, respectively, are very close to 1, suggesting a strong relationship between 
variables and the higher correlation between experimental and fitted values. This confirms that 
the proposed model has excellent predictive capability. 

Figure 2 shows the profiles of the predicted response TCR vs the four factors: (#, T, dp, and 
Np), coded X;, X2, X3, and X4, respectively. They indicate that the volume fraction (g) and the 
temperature (T) have a positive effect on the response. These results signify that TCR increases 
when X, and X, change from a lower level to a higher level. On the other hand, the particle size 
(dp) has a negative effect on the response. This means that the response reduces when this 
factor passes from low levels to high levels. Thus, it can be confirmed that CuO is the 
preeminent nanoparticle used to increase the TCR response. 


GRINE ap BENHAMZA 


HEAT TRANSFER 


TABLE 6 Prediction parameter estimates 


Term 


Intercept value à 
X; (0.01, 0.04) 

X, (295, 298) 

X3 (2.1e—8, 2.9e—8) 
Xà [L1] 

X3 [L2] 

X3 [L3] 

X1 XXa [L1] 

Xi X X4 [L2] 

X1 X X4 [L3] 

XX X3 

XX X4 [L1] 

X3 X X4 [L2] 

X2 X X4 [L3] 

X3X Xa [L1] 

FO (LA 

X3X X4 [L3] 


X1 X X2 X X3 

XX X2XX3 XX [L1] 
XX XX X3 X X4 [L2] 
X XX X X3 X X4 [L3] 


“Significant at P <.05. 


Estimation 


1.08037167 
0.03250583 
0.0061325 
—0.00874542 
—0.01965042 
—0.00307167 
0.02272208 
0.00395292 
0.00352167 
—0.00747458 
—0.00521292 
—0.00183625 
—0.0057375 
0.00757375 
0.01039917 
0.00451167 
—0.01491083 
0.00021375 
0.0076375 
—0.00785125 
0.00051833 
0.00454833 
—0.00506667 
0.00418375 
—0.0023925 
—0.0049925 
0.007385 


Asian Research 


Std error 


0.00152043 
0.00152043 
0.00152043 
0.00152043 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00152043 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00215021 
0.00152043 
0.00215021 
0.00215021 
0.00215021 


Wiley? 


t Ratio P>|t| 
710.57 <.0001* 
21.38 <.0001* 
4.03 .0157* 
75 .0045* 
—9.14 .0008* 
143 .2263 
10.57 .0005* 
1.84 .1399 
1.64 .1768 
—3.48 .0254* 
—3.43 .0266* 
—0.85 .4412 
—2.67 .0559 
3.52 .0244* 
4.84 .0084* 
21 .1039 
—6.93 .0023* 
0.1 .9256 
3.55 .0238* 
—3.65 .0217* 
0.24 .8214 
212 .1019 
—2.36 .078 
2.75 .0513 
—1.11 .3282 
2,92 .081 
3.43 .0264* 


Maximizing desirability is performed to obtain optimal factor settings for a maximum TCR. 
The profiler in Figure 3 displays the following optimal settings: 0.042 for volume fraction, 
294.73K for temperature, 21.27 nm for particle diameter, and CuO for the nanoparticle type, 
thus giving an estimated response equals to 1.162. 


5.2 | Validation of the proposed model 


To validate the proposed model of nanofluids conductivity ratio, a comparison with 
experimental and theoretical models results is performed. Figures 4-9 show the effects of 
different parameters on the nanofluids TCR of experimental data set compared with the 


proposed model. 
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TABLE 7 Analysis of variance ANOVA results for experimental design 


source SS DF MS F value P>F 
Model 0.04270618 19 0.0022477 40.5130 -0013* 
X1 (0.01, 0.04) 0.02535910 1 0.0253591 457.0790 < .0001* 
X2 (295, 298) 0.00090258 1 0.0009026 16.2684 .0157* 
X3 (2.1e—8, 2.9e—8) 0.00183558 1 0.0018356 33.0849 -0045* 
x: 0.00729494 2 0.0036475 65.7429 -0009* 
REX 0.00067118 2 0.0003356 6.0487 -0617 
XX; 0.00065219 1 0.0006522 11.7552 .0266* 
XX 0.00074922 2 0.0003746 6.7521 -0522 
X3XX4 0.00280665 2 0.0014033 25.2939 .0054* 
XXI XX 0.00096015 2 0.0004801 8.6530 .0352* 
XX X3X X4 0.00037302 2 0.0001865 3.3617 1391 
ROGUE 0.00042009 1 0.0004201 7.5718 .0513 
XX XX3X X4 0.000681 50 2 0.0003407 6.1418 -0603 
Residual 0.00022192 4 0.0000555 

Total 0.04292811 23 

R° 0.99483035 

Adj R° 0.97027453 


Abbreviations: DF, degree of freedom; MS, mean square; SS, sum of square. 
“Significant at P <.05. 


The effect of volume fractions on the predicted TCR with different diameters of 
nanoparticles of CuO-water nanofluids and at temperatures of 295K and 298K are given in 
Figure 4A,B, respectively. As it can be seen in these figures, the proposed model is in good 
agreement with the experimental data which are extracted from Aminia,** Zerradi et al,“ and 
Khedkar et al.“ 

Figure 5A shows the proposed model's ability to predict the TCR of ALO;-water nanofluid at 
298K, compared with the experimental data.** Thus, the proposed model has an accurate 
prediction of the TCR for different particle volume fractions. 


1.071885 
Ü [1.06447, 
F 1.0793] 


LUS Ge . CESR = a a 
SISSSS ALALSAS à 3 8 8 7 > 3% 
= = = à & & à à à 
0.02272 296.381 2.5e-8 L2 
TK) dp (m) Np 


FIGURE 2  Predict TCR vs X,(9), X(T), X:(dp), and X, (Np) (Red lines, predicted TCR; blue lines, 95% 
confidence interval of the mean TCR). TCR, thermal conductivity ratio [Color figure can be viewed at 
wileyonlinelibrary.com| 
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FIGURE 3 Prediction profiles at the optimum settings. TCR, thermal conductivity ratio [Color figure can be 
viewed at wileyonlinelibrary.com] 


Similarly, Figure 5B illustrates a good agreement of the TCR of TiO,-water nanofluid of the 
proposed model compared with the experimental data of Duangthongsuk and Wongwises, * 
Yiamsawasd et al,*” and Zerradi et al** at the same conditions of temperature and volume 
fraction for different nanoparticles diameters. 

The effects of temperature variation (T) on the nanofluids TCR of CuO-water and Al,03- 
water for nanoparticles diameter of 29 nm are presented in Figure 6A.,B, respectively. The 
predictions of the TCR of the proposed model are within an acceptable accuracy for an available 
experimental data.** 

Results are shown in Figure 7 of CuO-water nanofluids TCR at an ambient temperature 
of our proposed model are compared to experimental results, that are extracted from 
Aminia** as well as to correlations developed by Chon et al,*° Koo and Kleinstreuer,”! and 
Li and Peterson, *? Hamilton-Crosser classical model** and to a neural model.** As it can 
be seen, the proposed model provides a better matching than the results obtained by Li 
and Peterson,”* Koo and Kleinstreuer,*! and neural models**. Similarly, the correlation 
developed by Chon et al*° fits well the experimental data. Furthermore, one can notice 
that the Hamilton-Crosser relationship provides good results for small values of the 
particle volume fraction. 

In our proposed model, the variation of nanofluids TCR as a function of temperature is 
compared to the results of various correlations as well as to data set of Aminia,** and that is for 
CuO in the water of a 29 nm diameter and particle volume fraction of 4% (Figure 8). As it can be 
clearly seen, the proposed model generates a better fitting with experimental data. Nevertheless, 
both relationships developed by Koo and Kleinstreuer”' and Chon et al”° present relatively fine 
results. In contrary, Li and Peterson®* model cannot predict well the thermal conductivity 
enhancement trend of the experimental results of Aminia** It is also interesting to note that the 
Hamilton-Crosser model”* remains an acceptable prediction of nanofluids effective thermal 
conductivity at least at ambient temperatures. 

Figure 9 depicts a comparison between the thermal conductivity obtained from experiments 
and that obtained from the present model. The experimental data sets used in this figure are 
from Kazemi-Beydokhti et al,” Yiamsawasd et al,*” and Khedkar et al*° and the predictions’ 
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FIGURE 4 Effects of nanoparticles volume fraction on the thermal conductivity ratio of CuO-water 
nanofluid for different diameters at (A) 295K and (B) 298K, compared with the proposed model [Color figure 
can be viewed at wileyonlinelibrary.com| 
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FIGURE 5 Effects of nanoparticles volume fraction on the nanofluids thermal conductivity ratio of (A) 
ALO;-water and (B) TiO,-water for different diameters, compared the proposed model [Color figure can be 


viewed at wileyonlinelibrary.com] 
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FIGURE 6 Effects of temperatures on the nanofluids thermal conductivity ratio of (A) CuO-water and (B) 
ALO;-water compared with the proposed model [Color figure can be viewed at wileyonlinelibrary.com| 
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FIGURE 9 Effects of particle volume fraction on the nanofluids thermal conductivity ratio for different 
nanoparticles compared to the proposed model [Color figure can be viewed at wileyonlinelibrary.com] 


results of the present model are within a good accurateness compared with experimental 
results. 


6 | CONCLUSION 


In this study, a full factorial design approach was performed to estimate the thermal conductivity 
enhancement of metallic oxide-water nanofluids of ALO3, TiO;, and CuO nanoparticles. A new 
model of TCR of oxide-water nanofluids is developed using expanded and extensive experimental 
data. This model is applicable to the particle volume fraction range of 1%-4%, at low temperatures 
of 295K to 298K and particle diameters of 21-29nm. The factors that affect the thermal 
conductivity enhancement ratio are identified and analyzed by the ANOVA method. These 
factors are volume fraction, temperature, nanoparticle average diameter, and nanoparticle type. 
The results demonstrate that the particle volume fraction has the most significant effect on 
nanofluids thermal conductivity enhancement compared with other factors. An increase in the 
particle volume fraction and temperature has a favorable effect on the nanofluids thermal 
conductivity enhancement, while the higher particle size tends to reduce it. Furthermore, it is 
found that the nanoparticles volume fraction has a great impact on the thermophysical properties 
of nanofluids. Thus, the optimum nanoparticle concentration is essential to be determined in 
forming nanofluids that can enhance thermal systems performance. Finally, the evaluation of the 
results shows an excellent correlation between predicted and observed TCR responses, indicating 
good accordance between the suggested model and the experimental data. In addition, the 
obtained results confirm the effectiveness and capability of the proposed model in predicting 
nanofluids thermal conductivity with high accuracy. 
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